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Sommaire 


L’enyoi dans les matériaux d’impulsions élastiques ultrasoniques de forme connue, au 
moyen de céramiques piézoélectriques, et leur réception aprés propagation, pose le pro- 
bléme de la détermination des fonctions de transfert 4 l’émission et a la réception. Ces 
fonctions sont calculées dans le cas de pastilles piézoélectriques cylindriques vibrant 
parallélement 4 leur axe. On peut ainsi définir les impulsions électriques d’attaque et 
déterminer les impédances acoustiques des matériaux a appliquer sur les faces des pastilles 


pour obtenir des impulsions de forme désirée. 


Zusammentassung 


Beim Senden elastischer Ultraschallimpulse bekannter Form mittels piezoelektrischer 
Keramiken in ein zu untersuchendes Material und beim Empfang dieser Impulse nach 
Durchlaufen der Probe tritt die Frage nach der Bestimmung der Ubertragungsfunktionen 
auf. Diese Funktionen werden fiir zylindrische piezoelektrische Scheibchen, die parallel zu 
ihrer Achse schwingen, berechnet. Man kann so die elektrischen Sendeimpulse festlegen 
und die akustischen Impedanzen der Stoffe bestimmen, die man an den Flachen der Scheib- 
chen anbringen mu}, um Impulse der gewtinschten Form zu erhalten. 


Summary 


The transmission into materials of ultrasonic elastic impulses of known form, by means 
of piezoelectric ceramic materials, and their reception after propagation, presents the 
problem of knowing the transfer coefficient on sending and on receiving. These coefficients 
have been calculated for a piezoelectric cylindrical disc vibrating parallel to its axis. It is 
thus possible to define the sender electrical impulse and to determine the acoustic impe- 
dance of the material to be applied on the opposite faces of the model being studied, in 


Patines Hl DE CAPTEURS PIEZOELECTRIQUES 


order to obtain transmitted impulses of the desired form. 


1. Notations Utilisées 


Capacité propre des pastilles 
Constante élastique longitudinale 


~~. 


nom 


Courant électrique aux bornes 
Vet 


Charge électrique aux bornes 


: : Surface 
onstante élastique longitudinale du Ti 
Constante élastique longitudinale du TiO , 
Constante piézoélectrique Cite oV 
D 1015 dynes le Ti0.B * ; 
Pecan lsmb Pe eae a T Fonction de transfert de |’émetteur 
,3 coulombs A ; 
Epai Fee criles onction de transfert du capteur 
aisseur des ae F ; 
bE P U_ Différence de potentiel aux bornes de l’émetteur 
Fréquence Decne ‘al 
Fré See aly Rag ee U’_Différence de potentiel aux bornes du capteur 
réquence de résonan : ae : 
F aS : ee eee ae V, Vitesse longitudinale de propagation dans une 
orce exercée sur la face | de l’émetteur F . 
: : tige de TiO.Ba 
(sens positif = sens de l’axe électrique) . a \ 
Force exercée sur la face 2 de |’émetteur Uys tesse Gee uarticiles ee Be lee a 
(sens positif = sens de l’axe électrique) enor eae ech ace axe Clecirique) 
Force exercée sur la face 1 du capteur v2 Vitesse des particules sur la face 2 ; 
(sens positif = sens de I’axe électrique) (sens positif = sens de l’axe électrique) 
Force exercée sur le face 2 du capteur Z  Impédance acoustique 
(sens positif = sens de l’axe électrique) Z, Umpédance acoustique sur la face | 
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Z, Impédance acoustique sur la face 2 
® Rapport du transformateur du schéma équi- 


valent (®©=DC)) 
2) Densité du TiO;Ba 
o Densité 
@® Pulsation 
@p Pulsation de résonance en faces libres. 


2. Introduction 


Dans les techniques de modéles sismiques on 
cherche a obtenir des impulsions ultrasoniques de 
forme appropriée, et pratiquement aussi simple que 


3. Emetteur 


La fonction de transfert de |’émetteur est le rap- 
port de la force transmise aux matériaux au contact, 
a la tension électrique d’attaque aux bornes de la 


pastille. 
Du cété de la face 2, elle s’écrit: 
Taya ea 
DU 


ot ®=DC, est le rapport du transformateur repré- 
senté dans schéma équivalent. 

La résolution des équations du réseau donne la 
fonction de transfert: 


Zz (Z,;+j Zot) 


LEG) == 


5 : 
a eae 
joC,Zo 


possible, pour permettre une bonne analyse des 
phénoménes de propagation. Ceci pose le probléme 
de la détermination des fonctions de transfert des 
émetteurs et capteurs piézoélectriques. 

Ces fonctions de transfert sont déterminées dans 
la présente étude pour des pastilles piézoélectriques 
cylindriques vibrant longitudinalement, dont les 
axes électriques et mécaniques sont confondus, dont 
le diamétre est faible devant la longueur d’onde des 
impulsions qui prennent naissance, et dont les faces 
sont au contact de matériaux homogénes d’impé- 
dances acoustiques connues. 

Pour les pastilles ainsi définies, le systeme d’équa- 
tions reliant le courant, la tension électrique aux 
bornes des pastilles, la vitesse des particules et la 
force appliquée sur chacune des faces s’écrit [5] 

v9 =v; cos © +j re N) ain 


Vo Zy Vo 


° 


Ps +DO SF) DO) toa 0 Ze sith 
Vo Vo 


joc U =i1+ DC) (v.—2) - (1) 


On déduit du systéme d’équations (1) le circuit 
électromécanique équivalent a la pastille piézoélec- 
trique (Fig. 1) dans l’analogie qui fait correspondre 
respectivement un courant et une tension électri- 
ques a la vitesse des particules et a la force appli- 
quée [7]. Zy= S V9 co est ’impédance mécanique 
caractéristique des pastilles. Z, et Z, , quotients pour 
chacune des faces de la force appliquée a la vitesse 
des particules, sont les impédances des matériaux au 
contact. Ces impédances sont réelles dans le cas de 
barres appliquées contre les faces de la pastille. 
Dans le cas de plaques ou de blocs, elles sont com- 
plexes [6]. 


sin eV 


; (2) 


) Zit Za+2j 20) + (Zr +i Zo Cari Zao) 


Fig. 1. Circuit électromécanique équivalent a la pastille 
piézoélectrique. 


3.1. Impédances Z, et Zz réelles 


Les variations de 7'(@/@g) sont représentées sur 
la Fig. 2a en module et en phase pour différentes 
valeurs réelles de l’impédance Z,, Z, étant choisie 
égale a Z). Pour une capacité propre Cy = 12 pico- 
farads, l’ordonnée T (w/z) = 1 correspond a 2-10? 
dynes/volt. Remarquons que pour les faibles valeurs 
de w/z la force transmise au matériau croit avec Z, 
et que pour Z, voisin de Z), T(@/wa) est sensible- 
ment linéaire jusqu’a @/Wp = 0,2. Si l’on désire une 
fonction de transfert linéaire jusqu’a des fréquences 
relativement hautes, il y a donc intérét a utiliser des 
pastilles émettrices minces, a fréquence de résonance 
élevée. 


3.2. Impédance Z, réelle, Z, complexe 


Z,; étant choisie égale a Zy pour les raisons énon- 
cées ci-dessus, il est intéressant de voir dans quelle 
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Module 


Z, parametre réel. 


_ mesure la fonction de transfert est déformée lorsque 
i la face 2 est au contact de corps autres que des 
_ barres, c.a.d. lorsque Z est complexe. Les variations 
_ de T(@/@g) sont représentées sur la Fig. 2 b, en mo- 
: dule et en phase, pour différentes valeurs de |’impé- 
_ dance complexe Z, , Z; étant choisie égale a Z, . 

Le terme imaginaire n’affecte que les amplitudes, 
et de facon relativement limitée. 


4. Capteur 


La fonction de transfert du capteur est le rapport 
de la tension électrique recueillie aux bornes de la 
pastille 4 la force exercée sur ses faces. 


t 
i 
if (eee 2=(142)) Zo 


Phase ae Ze (142))Z5 
SEZ ez 


nia 


nid 


aa 


Fig. 2 b: Fonction de transfert de l’émetteur: Z,=Z, ; 
Z» parameétre complexe. 


Pour la face 2 elle s’écrit: 
DU’ 


ou @ est défini comme plus haut. 


T’(w) = 


b) 


En supposant, pour simplifier, que la tension U’ 
est mesurée aux bornes d’une forte impédance, ce 
qui généralement est le cas dans les expériences sur 
modeles ultrasoniques, la résolution des équations 
du réseau donne la fonction de transfert 
fi (@) tS Se! 5 ef gis uk : 

j@CyZ, oe TS ie 


sino e/Vy (j Zyt+Z;) 
(3) 


ad! jes eS ee 
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Fig.3 a. Fonction de transfert du capteur: Z,=Z,) ; 


Z, parametre réel. 


vi tg we . 
2V, 


2.1. Impédances Z, et Zz réelles 


Les variations de 7’ (~/mg) sont représentées sur 
la Fig. 3 a en module et en phase pour différentes va- 
leurs réelles de l’impédance acoustique Z,, Z, étant 
choisie égale a Zy. Pour une capacité propre 
Cy = 12 picofarads, l’ordonnée 7” (@/mg) = 1 corres- 
pont 45°10 4 volts/dyne. 

Remarquons que pour les faibles valeurs de 
w/z la tension électrique délivrée par le capteur 
croit avec-Z,-et que pour Z,;=Z), T’ (w/@p) admet 
une pente sensiblement nulle tant que @/@g reste in- 
férieur a 0,3. La condition Z,=Z) permet done d’ 
obtenir une réception sans déformation. 


On remarque également que |’amplitude est in- 
versement proportionelle a la capacité Cy. 


0,5 


Z,=(1*33) Z, 
Z=(1*5 ) Zo 
Z=(*th ). Ze 
Z=(* Kw) % 


Fig.3 b. Fonction de transfert du capteur: Z,;=Zy ; 
Zz parametre complexe. 
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On utilisera donc des pastilles réceptrices relative- 
ment épaisses a condition que la fréquence de réso- 
nance en oscillations libres ne soit pas trop basse, ce 
qui aurait pour effet de limiter la bande passante. 


Bee Impédance Z, réelle, Z, complexe 


Lorsque la face 2 est au contant de corps 
autres que des barres, Z, est complexe. Z, étant 
choisie égale a Zy pour obtenir une réception sans 
déformation, les variations de 7’ (w/w) sont repré- 


| sentées sur la Fig. 3b en module et en phase, pour 


différentes valeurs de |’impédance Z, complexe. Ici 
encore le terme imaginaire affecte surtout les ampli- 
tudes. 


5. Exemples 


Des impulsions élastiques de formes particuliére- 
ment simples, comprises dans la bande de 10 a 


50 ke/s, sont utilisées pour la réalisation de sismo- 


grammes synthétiques par le procédé de la barre 
analogique [1]: une pastille émettrice et une pastille 
réceptrice accolées permettent d’introduire des impul- 
sions a une extrémité de la barre et de capter les 
arrivées d’énergie réfléchie aprés propagation. La 
barre analogique est appliquée d’un cété du dis- 
positif émetteur-récepteur; une barre de diametre 
constant et d’impédance connue est appliquée de 
Yautre cote. En choisissant les impédances mécani- 
ques caractéristiques des matériaux au contact sen- 
siblement égales a celles des pastilles, les écarts sur 
les valeurs du produit @ V étant au besoin compensées 
par des sauts de diamétre appropriés [3], et en 
prenant une fréquence de résonance en face libre de 
250 ke/s pour les pastilles, les fonctions de transfert 
conservent a |’émission une pente | et a la réception 
une pente nulle, jusqu’a des fréquences voisines de 
50 ke/s; on obtient ainsi a |’émetteur une force sur 
la face 2 sensiblement proportionnelle a la dérivée du 
signal électrique d’attaque, et au capteur un signal 
électrique sensiblement proportionnel a la force 
exercée sur la face 2. La Fig. 4 est la photographie sur 


- 


(a) 


(b) 


Fig.4. Réponse (b) du capteur a un signal d’attaque 
(a) de 20 microsecondes de durée. - 


lécran de loscilloscope de la réponse du capteur a 
un signal d’attaque particuliérement simple. La Fig.5 


- représente la fonction de transfert expérimentale de 


lensemble émetteur-capteur accolés, comparée a la 
fonction de transfert théorique. La bonne concor- 
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{T(e) T'(o) 
MILLIVOLTS /VOLTS 


réponse theorique 
atténuée 50 fois 
== - reponse expérimentale 


Fig.5. Fonction de transfert expérimentale de len- 
semble émetteur-capteur, pastilles accolées, com- 
parée a la fonction de transfert théorique. 


ZL; —Ly . 


dance des résultats théoriques et expérimentaux 
semble confirmer les hypothéses faites sur la nature 
des impédances acoustiques. 

Les Figs. 2b et 3b laissent prévoir que méme 
lorsque les impédances chargeant les faces des pastil- 
les n’ont plus des valeurs réelles, mais des valeurs 


complexes, l’impulsion doit s’écarter assez peu de 
la forme simple représentée sur la Fig. 4. C’est ce 
qu'on voit sur la Fig. 6, tirée d’une étude expérimen- 


ty 


WEEE 


i 
VS. 


Fig. 6. Exemple d’utilisation de Vimpulsion de la Fig. 4 
a Vétude des diffractions créées par une fente. 
Echelle verticale: 10 microsecondes entre lignes 
de temps. 


tale sur les diffractions causées dans un solide par 
un écran semi-infini parfaitement réfléchissant [4]. 
Un émetteur E et des capteurs Ry a Ryy ont été dis- 
posés sur les tranches opposées d’une plaque de 
plexiglas entamée par une fente d’air. On a repré- 
senté la trace synthétique pour chacune des positions 
du capteur. L’extrémité I de la fente devient 
de au Vonde 
longitudinale incidente issue de E parvient en ce 


source diffraction moment ou 


point, et émet une onde de dilatation et une onde 
de distorsion diffractées. La bonne qualité des 
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impulsions permet de séparer nettement les arrivées 
directes (1), les arrivées de dilatation diffractées 
(2), les arrivées de distorsion diffractées (3), et les 
réflexions multiples (4) et (5). 


6. Conclusion 


Dans les études sur modéles sismiques, les im- 
pulsions de forme simple nécessaires 4 une bonne 
analyse des phénomenes de propagation sont obtenues 
en attaquant |’émetteur par une impulsion électrique 
de forme appropriée, et en disposant sur les faces 
des pastilles émettrices et réceptrices des matériaux 
d’impédance acoustique convenable. La connaissance 
des fonctions de transfert de l’émetteur et du capteur 
permet de déterminer parmi les matériaux dispo- 
nibles et les signaux d’attaque pratiquement réali- 
sables, ceux qui permettent d’obtenir au capteur 
Vimpulsion désirée. 

Les auteurs remercient Mesdemoiselles F. Bar- 
rier et H. Onacoiry, de |’Institut Francais du Peé- 


. trole, qui ont effectué les calculs numériques. 
(Regu le 4 Décembre 1960.) 
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ABSORPTION LONGITUDINALER ULTRASCHALLWELLEN 
IN ALUMINIUM BEI HOHEN TEMPERATUREN 


von K. TsapEen 


III. Physikalisches Institut der Universitat Gottingen 


Zusammenfassung 


Die Absorption longitudinaler Ultraschallwellen wurde in zwei polykristallinen Alu- 
miniumlegierungen und in einem Aluminiumeinkristall mit einem Echo-Impuls-Verfahren 
ermittelt. Die Messungen wurden im Frequenzbereich von 5 bis 72 MHz und bei Tempe- 
raturen zwischen Zimmertemperatur und 400 °C ausgefiihrt. Fiir alle Proben ergab sich 
ein quadratischer Anstieg der Absorption mit der Frequenz. Wahrend im Einkristall die 
Absorption mit der Temperatur kontinuierlich zunimmt, erhalt man in den polykristallinen 
Proben einen plotzlichen Anstieg der Absorption bei Temperaturen oberhalb 300 °C. 

Die bei dem angewandten Verfahren auftretende geometrische Ausbreitungsdampfung 
wurde unter speziellen Voraussetzungen berechnet und mit Messungen in einem Quarz- 
glaszylinder verglichen. Dabei ergab sich, da8 die Ausbreitungsdampfung kleiner als er- 

_ wartet und fiir die Absorptionsmessungen in Aluminium im betrachteten Frequenzbereich 
praktisch zu vernachlassigen ist. 


Sommaire 


A Vaide de la méthode utilisant les échos des impulsions sonores on a déterminé |’ab- 
sorption des ondes ultrasonores longitudinales dans deux alliages d’aluminium, puis dans 
un cristal unique d’aluminium. Les mesures ont été faites dans une gamme de fréquences 
de 5 a 72 MHz et pour des températures variant depuis la température normale des appar- é 
tements jusqu’a 400 °C. Pour tous les échantillons on obtint un accroissement quadratique i 
de l’absorption avec la fréquence. Mais tandis que |’absorption augmente continfiment avec F 
la température dans le cas du cristal unique, elle présente dans le cas des échantillons 
polycristallins un accroissement soudain pour des températures dépassant 300 °C. 

En utilisant cette méthode de mesure on a calculé |’atténuation géométrique de la pro- 
pagation sous des hypothéses particulieres et on a comparé le résultat avec celui des 
mesures faites sur un cylindre de quartz fondu. On a trouvé que |’atténuation de propa- 
gation est plus faible qu’on ne l’aurait attendue et que, dans les mesures d’absorption 
dans laluminium pour le domaine des fréquences considéré, elle est pratiquement négli- 
geable. 


Summary 


The sound-echo impulse method has been used to measure the absorption of longitudinal 
ultrasonic waves in two alloys of aluminium and in a single crystal of aluminium. The meas- ee 
urements haye been made in a frequency range from 5 to 72 Mc/s and for temperatures 
varying from normal room temperature up to 400 °C. For all the specimens a quadratic 
inerease of absorption with frequency was obtained. However, whereas there was a conti- 
nual increase of absorption with temperature in the case of the single crystal, there was a 
sudden rise for temperature above 300 °C in the polycrystalline specimens. ie 

In using this method of measurement the geometrical attenuation of the transmitted 
waves, making particular assumptions, has been calculated and the result compared with 
measurements made on a quartzglass cylinder. It has been found that the attenuation of the 
transmitted waves is less than to be expected and that, in the absorption measurements in 
the aluminium for the frequency range considered, it is practically negligible. 


1. Einfiihrung 


Messungen der Ultraschallabsorption in Metallen 
bei Frequenzen oberhalb 1 MHz wurden bislang 
iiberwiegend bei Zimmertemperatur ([1] bis [3]) 
und bei tieferen Temperaturen [4], [5] durchge- 
fiihrt. Bei hoheren Temperaturen liegen dagegen nur 
vereinzelte Ergebnisse vor. Mit bei tieferen Fre- 
quenzen anwendbaren Verfahren wurde jedoch schon 
friih die Dampfung in Metallen bei hoher Tempe- 


ratur untersucht [6], [7]. Auf Grund der hierbei 
erzielten Ergebnisse und einer Messung in einem 
Bleieinkristall bei 26,5 MHz entwickelte Mason [8] 
eine Theorie, in der er die mit zunehmender Tempe- 
ratur wachsende Dampfung auf den Einflu8 der 
Versetzungen zuriickfiihrt. Er macht fiir die Damp- 
fung bei hohen Temperaturen diejenigen Verset- 
zungsschleifen verantwortlich, die durch die Schall- 
welle von den als Haftpunkten dienenden Fremd- 
atomen momentan hinweggerissen werden k6nnen. 


128 


Nicht haftende Versetzungsschleifen schwingen nicht 
koharent mit dem Gitter und entziehen der Schall- 
welle Energie. Die Wahrscheinlichkeit des Losbre- 
chens der Versetzungsschleifen wachst mit steigender 
Temperatur. Als Maf fiir die Dampfung verwendet 
Mason den reziproken Wert der Giite Q des Mate- 
rials, fiir den er bei hohen Temperaturen folgende 
Beziehung angibt: 


<f- = Ce —H/RT 


Q 
(C Konstante, H Aktivierungsenergie, R Gaskon- 
stante, T absolute Temperatur). Der Ausdruck 1/Q 
hangt mit dem Absorptionskoeffizienten « einer fort- 
schreitenden ebenen Welle mit dem Druck 


p= poe (et-*2) 


(1) 


(po Druck an der Stelle x=0O zur Zeit t=O, 
ik=a+if Ausbreitungskonstante, « Absorptions- 
koeffizient, 6 =2x/A WinkelmaB) auf die einfache 
Weise 


ete pices 


(2) 


Q <a f 
(v Schallgeschwindigkeit, f Frequenz) zusammen. 
Auf Grund der Beziehungen Gl. (1) und Gl. (2) 
mute die durch Losbrechen der Versetzungsschlei- 
fen verursachte Absorption linear mit der Frequenz 
anwachsen. Bei den bisher vorliegenden Absorptions- 
messungen in Metalleinkristallen [3], [9] wurde 
jedoch eher ein quadratischer Anstieg der Absorp- 
tion mit der Frequenz gefunden. Lams, Repwoop 
und SurernsHierer fanden ferner in Germanium 
und Silizium zwischen Temperaturen von — 154 °C 
bis +197 °C nur eine geringe Temperaturabhingig- 
keit. 

Als einzige bekannte Messung bei noch héheren 
Temperaturen sei die von Zucker [10] bei 5 MHz 
in einer polykristallinen Aluminiumlegierung bis 
375 °C erwahnt. Er findet fiir Longitudinalwellen 
einen sehr starken Anstieg der Absorption oberhalb 


300 °C. 


2. Das MeBverfahren 


Fiir Messungen der Ultraschallabsorption im Fre- 
quenzbereich von 3 bis 72 MHz wurden zwei Appa- 
raturen entwickelt, die nach dem Echo-Impuls-Ver- 
fahren arbeiten. Die mit diesen Apparaturen zu 
messenden Proben hatten zylindrische Form, ihre 
Stirnflachen waren auf 1 parallel und auf einen 
Interferenzring plan geschliffen und poliert. Auf 
eine Stirnflache der Probe wurde als elektromecha- 
nischer Wandler eine in X-Richtung geschnittene 
Quarzscheibe aufgedriickt, wobei als Kontaktsub- 
stanz ein dinner Film aus temperaturbestaéndigem 
Silicono! (Dow Corning Corporation DC 710) ver- 
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wendet wurde. Die Probe diente gleichzeitig als die 
eine Elektrode zur Anregung des Quarzes, die Ge- 
genelektrode aus Messing wurde mit einer Feder 
gegen den Quarz gedriickt und fettfrei gehalten, so 
daf§ der Quarz nach hinten akustisch praktisch an 
Luft grenzte und in diese Richtung keinen Schall 
abstrahlen konnte. 

Der Quarz wurde mit Schwingungsimpulsen von 
etwa 6us Lange in seiner Grundschwingung oder 
ungeraden Oberschwingung angeregt. Die durch 
mehrfache Reflexion an beiden Stirnflachen der 
Probe entstehende Echofolge wurde durch den glei- 
chen Quarz registriert, in einem Uberlagerungs- 
empfanger verstarkt, demoduliert und auf dem Oszil- 
lographenschirm sichtbar gemacht. 

Fur die genaue Ablesung der Hohe der Echo- 
impulse wurde zu einer gegen den Ausléseimpuls 
veranderbaren Zeit durch einen weiteren Oszillator 
ein Vergleichsimpuls erzeugt, der tiber eine Eich- 
leitung dem Empfanger zugeftihrt wurde. Die Hohe 
des Vergleichsimpulses wurde jeweils gleich der 
Hohe des zu messenden Echoimpulses gemacht, dann 
konnte die Amplitudenabnahme der Echos auf der 
Kichleitung in dB abgelesen werden. 

Die Quarzhalterung und die Zufiihrung des 
Kabels wurde aus temperaturbestandigem Material 
ausgeftihrt, so daf} die Probe in einem elektrischen 
Rohrofen auf die erforderlichen Mefitemperaturen 
bis zu 400 °C gebracht werden konnte. Die Tempe- 
raturregelung erfolgte mit einem Kontaktthermo- 
meter, die Temperaturmessung mit einem Wider- 
standsthermometer. 


3. Die geometrische Ausbreitungsdimpfung 


Fur die Erzeugung von ebenen Wellen ist erfor- 
derlich, dafS der Durchmesser des Schwingungsquar- 
zes und der MeBprobe jeweils méglichst viele Wel- 
lenlangen betragt. Bei tiefen Frequenzen wird man 
jedoch diese Forderung nur ungeniigend befriedigen 
konnen und Einflisse der seitlichen Begrenzungen 
des Quarzes oder der Probe beriicksichtigen miissen. 

Wahlt man den Durchmesser der MefSprobe gleich 
dem Durchmesser des Quarzes, so betrachtet man 
die Schallausbreitung im stabformigen Medium. Fiir 
dieses bereits von Mason und McScrimin [1] benutzte 
Verfahren haben Repwoop und Lams ([11] bis 
[13]) die Ausbreitung der verschiedenen Moden 
diskutiert. Jeder Mode entspricht eine Amplituden- 
verteilung tiber den Querschnitt des Stabes. Im all- 
gemeinen gelingt es nicht, die Amplitudenverteilung 
so herzustellen, da nur eine der betrachteten Moden 
angeregt wird. Das hat zur Folge, daB auf Grund 
einer Interferenz der mit verschiedener Schallge- 
schwindigkeit sich ausbreitenden Moden dem expo- 
nentiellen Abfall der Echofolgen periodische Schwan- 
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kungen tiberlagert werden. Aufierdem macht sich die 
zylindrische Begrenzung der Probe auch gerade fir 
die Ausbreitung von Impulsen storend bemerkbar. 
Bei schragem Einfall der Longitudinalwellen auf die 
zylindrische Manteloberflache findet eine Aufspal- 
tung in longitudinale und transversale Wellen statt, 
die unter verschiedenen Winkeln reflektiert werden. 
Abgesehen davon, dafi der Longitudinalwelle durch 
diesen Vorgang Energie entzogen wird, treten bei 


tiefen Frequenzen (<10MHz) Storimpulse auf, 


welche dadurch entstehen, dai die Transversalwellen 
bei erneuter Reflexion wieder in Longitudinalwellen 
verwandelt werden, die wie die urspriingliche Welle 
laufen und gegeniiber dieser um gewisse Zeiten ver- 
zogert sind. 

Wahlt man den Durchmesser der Probe grofer als 
den des Quarzes, und zwar so gro}, daB die zylin- 
drische Begrenzung der Probe keinen Einflu8 auf 


‘die Schallausbreitung hat, so betrachtet man die 


Schallabstrahlung eines kreisformigen Senders in 
einen unendlichen Halbraum. Dieses Verfahren 
wurde von Sexi, Granato und Truett [14] disku- 
tiert und bei Messungen in Germaniumkristallen [3] 


- angewandt. Man muf in diesem Fall eine durch Beu- 
.gung an der Begrenzung des Quarzes auftretende 


geometrische Ausbreitungsdampfung _ beriicksichti- 
gen, die sich unter speziellen Bedingungen rechne- 
risch ermitteln lat. 

Da ein in X-Richtung geschnittener Quarz aus- 
schlieBlich Longitudinalwellen abstrahlt und weder 
durch die Reflexion bei senkrechtem Einfall noch auf 
irgendeine andere Weise Transversalwellen erzeugt 
werden, lat sich das Problem mathematisch so be- 
handeln, als ob das Medium eine Fliissigkeit ware. 
Ferner wird angenommen, dafi der Quarz gleich- 
formig als Kolbenmembran schwingt, obwohl die 
Schwingungsverteilung einer Quarzplatte im allge- 
meinen komplizierter ist. Die an dem als Empfanger 
dienenden Quarz entstehende elektrische Spannung 
soll dem Integral des Druckes tiber die Empfanger- 


_ flache proportional sein. Das Medium wird als ver- 
lustfrei angenommen, das heift, man setzt a=0; 


k=22n/2. SchlieBlich sollen alle Reflexionen ohne 


Energieverlust vor sich gehen, so da} man sie nicht 
zu beriicksichtigen braucht und die Rechnung so 
durchfiihren kann, als ob getrennte Sender und Emp- 
finger vom Radius a einander im Abstand X dia- 
metral gegeniiberstiinden. 

Ausgehend von der bereits von Lord Ray eicu 
angegebenen Beziehung fiir das Schallfeld einer 
Kolbenmembran 


i(wt—kr) 
p= pe fs ay (3) 
20 r 
Fy 
_(®  Geschwindigkeitspotential im  Aufpunkt, 


_V,=9@/9n vorgegebene Normalschnelle der Kol- 
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benmembran, k = w/c = 2x/A Ausbreitungskonstante, 
r Abstand des Aufpunktes vom Flachenelement der 
Senderflache F,) erhalt man fir das Integral des 
Druckes tiber die Empfangerflache, das laut Voraus- 
setzung der am Empfanger auftretenden Spannung 
proportional sein soll: 


[ram--e [SF 


Fs F, 
i(wt—kr) 
Eine ber Soir dk. 
20 ir 


Pa oF 


O®@ 


dF, (4) 


SEKI, Granato und Truett [14] fihrten konse- 
quenterweise zunachst die Integration uber F, aus, 
berechneten also das Schallfeld der Kolbenmembran, 
muften jedoch dann bei der Integration tiber Fy auf 
graphische Methoden zuriickgreifen. Bei geeigneterer 
Wahl der Koordinaten auf F, und F, und durch eine 
andere Reihenfolge der Integration ist es jedoch 
moglich, das Problem auf rein analytischem Wege 


befriedigend zu losen. Man beschreibt die Sender- 
flache F, durch 
rears Os Se, 
Y,=SCOSP, OXpvS2n 
Zz, =ssing 


und die Empfangerflache Fy durch 
Om ho 2ia,, 
OSyS P (s,h) 


Cie ca 
Yo=scosp+hcosy, 
Z =ssing+hsiny, 

(a Radius des Senders und des Empfangers, x Ab- 


stand zwischen Sender und Empfanger). Dabei ist 
rahe and 


2% fir Ors: a he, 

Shep eestp 
W(s,h) = \ 2 are cos eS fir BMS Sr 
0 fir OSs <h-a. 


Durch unbestimmte Integration tiber 7p, y und s er- 
halt man nach Ejinsetzen der Grenzen 


V2+4e 


fo dF, = eS Tone: x 
2 
F L 
x (4 a> arc cos Voie oe — 
2a 
— V(r? —2?) (4@—7? +22) dr. (5) 
Hieraus erhalt man durch die Substitution 


r= (x? + 4a cos? 6)" und durch einmalige partielle 
Integration den von Wiiurams [15] angegebenen 
Ausdruck, aus dem Bass [16] eine Naherungsformel 
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entwickelte, 
2 
[rar- REDO We (6) 
k 
xl2 
- (x Ly | enik Vt + datcos* On? 0 dO ; 
0 


Fur den speziellen Fall von x =0 wird 


z V I, (2ka) 
dre 2¥ OV ¢ I, ) 7 
| p ; ERE te 
F,=F; 
2S Eee | 
aR 2ka i 


Dabei ist I, die Bresseusche und S, die Srruvesche 
Funktion erster Ordnung [17]. Fir groBe Werte 
von ka geht der Ausdruck Gl. (7) in aa? ocVy 
uber, das hei®t der iiber die Senderflache gemittelte 
Druck wird 


ig | Pah=2eVo, 
a 


F 


also gleich dem bekannten Wert fiir den Druck im 
Nahfeld einer Kolbenmembran. Es ist sinnvoll, den 
Grenzwert 1a”,@cV , der sich fir groBe Werte von 
ka ergibt, als Bezugswert fiir die zu berechnende 
Ausbreitungsdampfung 


hes | [pak | 
a \F. | 
—— = — 20 log —*—___! 8 
dB we aoc, (8) 


zu wahlen. 


Zur Berechnung des auf ma?o0cVy bezogenen 
Integrals Gl. (5) fiir den allgemeinen Fall x >0 sub- 
stituiert man zweckmabig 


r=2#+0 (u+1) 
mit 26= (a2 +4 a?)"". (9) 
Dann wird 


te Te fo dF, = ikd e ik @t 9) 4 
sf 


Qcro 
F; 


(10) 


+1 
x | ei“ (arecosv— V1 —v?) du 


-1 
mit 


v= — oe (ilies) 
a 


(11) 


Bei Beschrankung auf groBe Werte von x/a(x/a> 10) 
wird 0/a~a/x und (6/a)? so klein, daB man in 
Gl. (11) v=w setzen kann. Man kommt dann nach 
erneuter Substitution w= —cosa und partieller In- 
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tegration auf die Integraldarstellungen der Bessel- 
funktionen Ip und I, und erhalt 


My te os fe aR, =e**{1 —e ikrale x (12) 


naocV, 
2 2 
w(t) on (t2 
x x 


Py 
woraus sich die Ausbreitungsdampfung, definiert 
durch Gl. (8), in folgender Weise ergibt: 


se (2 \/"+ (13) 
x 


ues 10 log ||eos *% 
z 
o: s ka eae 
+. |sin ——= | See : 
x Be 


dB 

Den Verlauf dieser Ausbreitungsdampfung zeigt 
Bild 1. Im Nahfeld (x<15a?/A) hat sie den von 
Sexi, Granato und Truett angegebenen durch- 
schnittlichen Betrag von 1 dB pro a?/A. In Bild 2 ist 
die Abszisse im logarithmischen Mafstab aufgetra- 
gen. Man sieht den Ubergang zum Fernfeld 
(w>15a?/A), in dem die Ausbreitungsdémpfung 


ag a 


x 


ax 
5 


& lke 
0 5 10 5 
XA /a? 


Bild 1. Berechnete geometrische Ausbreitungsdimpfung 
nach Gl. (13). Die gestrichelte Gerade entspricht 
dem Durchschnittswert von a*=1 dB pro a?/A. 


° 03 1 3 0 


XA fa? ——» 


Bild 2. Berechnete geometrische Ausbreitungsdiémpfung 
nach Gl. (13). Die gestrichelte Gerade entspricht 
der fiir das Fernfeld (x >15a?/A) giiltigen 
Ausbreitungsdampfung nach Gl. (14). 


cae ee ole eS > 5 alain td coe 
“, a “Os . 
, > x bi 
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durch die Asymptote 
St 
dB 


dargestellt wird }. 

Zur experimentellen Ermittlung der Ausbrei- 
tungsdampfung wurden Messungen in einem Quarz- 
glaszylinder ausgefuhrt, der in Bild 3 dargestellt ist. 
Sein Durchmesser betrug 6 cm, seine Hohe 5 cm. 


2 
Se -20loz aa (14) 
x 


_ Bild3. Zwei MeBproben. Links der verwendete Quarz- 
glaszylinder (Durchmesser 6cm, Héhe 5 cm), 
rechts der Aluminiumeinkristall (Durchmesser 
3,6cem, Hohe 2,8cm). Die Schwingquarze der 
Grundfrequenz 4,8 MHz haben einen Durchmes- 
ser von 2cm. Als Unterlage dient eine Richt- 


platte aus Hartgestein, auf der das Parallel- 


schleifen der Metallproben erfolgte. 
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Bild 4. Gemessene Gesamtextinktion in geschmolzenem 
Quarz. Die gestrichelte Linie entspricht dem 
Durchschnittswert der berechneten geometrischen 
Arbeitsdimpfung im Nahfeld von 1 dB pro a?/A. 
X Grundfrequenz des Quarzes: 4,80 MHz, 

@ Grundfrequenz des Quarzes: 1,126 MHz, 
© Messungen von Repwoop, Lams. und Surein- 
SHLEIFER [9]. 
1 Die in Gl. (12) enthaltene Phasenbeziehung kann 
_ zur Berechnung des bei Schallgeschwindigkeitsmessun- 
_ gen mit dem Echo-Impuls-Verfahren durch Beugung an 
der kreisférmigen Begrenzung des Senders auftreten- 
_ den Fehlers dienen (vergleiche [20]). 


f 


7 


q 
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Der Durchmesser der verwendeten Schwingquarze 
betrug 2 cm, ihre Grundfrequenzen f=1,126 MHz 
und f=4,80 MHz. In die Mantelflache des Quarz- 
glaszylinders wurden ein Rechts- und ein Links- 
gewinde unterschiedlicher Steigung und Tiefe ein- 
geschliffen. Die hierdurch entstehende Struktur der 
Flache sollte bewirken, da der an die AuBenwand 
gelangende Schall zerstreut wird. 

Wie aus Bild 4 ersichtlich schlieBen die Messun- 
gen bei hohen Frequenzen gut an die von Lams, 
Repwoop und SuTeinsHueirer [9] in geschmolzenem 
Quarz erzielten Ergebnisse an. Die Materialdamp- 
fung in geschmolzenem Quarz ist so klein, dai ihr 
Anteil an der Gesamtextinktion bei Frequenzen 
unterhalb 10 MHz gegentiber dem der geometrischen 
Ausbreitungsdampfung vollig in den Hintergrund 
tritt. Man sieht, daf$ in diesem Frequenzbereich die 
Extinktion sehr stark von der Qualitat der verwen- 
deten Quarze abhangig ist, da die mit verschiedenen 
Quarzen erzielten Ergebnisse sehr unterschiedlich 
sind. Ferner zeigt ein Vergleich der MeSpunkte mit 
der gestrichelt eingezeichneten Geraden, die dem 
Durchschnittswert der geometrischen Ausbreitungs- 
dampfung im Nahfeld von 1 dB pro a?/A entspricht, 
daf die experimentell gefundene Extinktion im all- 
gemeinen viel kleiner ist als man allein auf Grund 
der unter den angegebenen Bedingungen berechne- 
ten geometrischen Ausbreitungsdampfung erwarten 
sollte. So ergab zum Beispiel eine Messung bei 
4,8 MHz, der Grundfrequenz des verwendeten Quar- 
zes, einen Extinktionswert von 6,0: 107? dB/cm, ob- 
wohl theoretisch allein die geometrische Ausbrei- 
tungsdimpfung im Nahfeld 0,115 dB/cm betragen 
sollte. Bild 5 zeigt die fiir diese Messung erhaltene 


Bild5. Im Quarzglaszylinder bei 4,8 MHz erhaltene 
Echofolge. 


Echofolge. Ihre Auswertung und den Vergleich mit 
dem unter den vorliegenden Bedingungen berech- 
neten Verlauf der geometrischen Ausbreitungsdamp- 
fung zeigt Bild 6. 

Durch Extrapolation in Bild 4 kann man den Wert 
der reinen Materialdimpfung in Quarzglas bei 


2; 


= 
0 ha? 50 100 140 


Bild 6. Gemessene (a) und nach Gl. (13) berechnete 
(b) Echofolge im Quarzglas, n: Echo-Nr. f= 
4,8 MHz, a=lcm, Lange der Probe 1=5 cm. 


4,8 MHz auf 5- 1074 dB/cm abschatzen. Ferner mu8 
bei diesen kleinen Extinktionswerten beriicksichtigt 
werden, das bei jeder Reflexion an dem Ende der 
Probe, an das der Quarz angekoppelt ist, ein gerin- 
ger Bruchteil der Schallenergie in elektrische Energie 
verwandelt wird. Aus dem Verhaltnis der Spannung 
des angelegten Sendeimpulses zu der Spannung des 
empfangenen ersten Echos wurde die hierdurch be- 
wirkte Dampfung zu 1,7- 10~? dB pro Reflexion am 
Quarzende ermittelt. Bei der vorliegenden Probe von 
5cm Lange entspricht das einer Daémpfung von 
1,7:107-3 dB/em. Die Schallabstrahlung in die 
Luft betragt etwa 2-10~4dB pro Reflexion, das 
entspricht umgerechnet einer Dampfung von 
4: 107° dB/cm. Zieht man alle diese Werte von dem 
bei 4,8 MHz ftir die Extinktion erhaltenen Wert 
6,0: 10° dB/em ab, so verbleibt fiir die geometri- 
sche Ausbreitungsdémpfung nur ein Betrag, der 
kleiner ist als 4: 10~3 dB/cm, das sind weniger als 
3,5% des nach der Theorie im Nahfeld zu erwarten- 
den Wertes. 


Nach den Untersuchungen von Bass [16] ent- 
spricht die geometrische Ausbreitungsdampfung in 
Wasser dem unter den angegebenen Voraussetzun- 
gen berechneten Verlauf wesentlich besser. Offenbar 
treffen diese Voraussetzungen fiir die Schallausbrei- 
tung im festen Korper nicht zu. Ferner ist denkbar, 
daB bei jeder Reflexion an dem Ende der Probe, an 
das der Quarz angekoppelt ist, eine Phasenkorrektur 
iiber den Querschnitt erfolgt und da8 aus diesem 
Grunde die Divergenz des Schallstrahles verringert 
wird. Die zylindrische Begrenzung der Probe hatte 
keinen EinfluS auf die Schallausbreitung. Dies 
wurde bei einem Versuch deutlich, bei dem die Rand- 
bedingungen durch Vollgiefen der eingeschliffenen 
Rillen mit Paraffinwachs verandert wurden. Dabei 
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konnte keine Veranderung im Verlauf der Echofolge 
festgestellt werden. 

Da nun die geometrische Ausbreitungsdampfung 
im Festkérper sehr viel kleiner ist, als man nach den 
angestellten Berechnungen vermuten sollte, kann 
man sie bei’ Messungen in Festkorpern mit entspre- 
chend héherer Materialdimpfung im vorliegenden 
Frequenzbereich oberhalb 5 MHz meistens vernach- 
lassigen. Die Messungen in Quarzglas dienten zur 
Abschatzung der bei solchen Messungen auftreten- 
den systematischen Fehler einschlieBlich der geome- 
trischen Ausbreitungsdampfung. 


4. Die Verlarmune der Impulse 


Bei den erhaltenen Echofolgen wurde eine mit der 
Anzahl der Reflexionen zunehmende Verformung 
der Impulse beobachtet. Diese Impulsverformung 
geschieht in der Hauptsache bei den Reflexionen an 
dem Ende der Probe, an das der Quarz angekoppelt 
ist. Sie ist um so grofer, je dicker der koppelnde 
Olfilm und je héher die Oberschwingung ist, in 
welcher der Quarz angeregt wird, und je mehr sich 
der Wellenwiderstand der Probe von dem des Quarz- 


kristalls unterscheidet. Aus Bild 7 wird deutlich, daB 


Oltilm 


Probe 


Bild 7. Reflektierende Flachen an dem Ende der Probe, 
an das der Quarz angekoppelt ist. 


Reflexionen an diesem Ende der Probe an drei ver- 
schiedenen Flachen geschehen, an der Riickseite des 
Quarzes I, der Vorderseite des Quarzes II und der 
Wand der Probe III. Zwischen den Flachen II und 
III befindet sich der Olfilm, dessen Dicke zwar klein 
gegentiber der Wellenlange ist, jedoch in jedem 
Falle ausreicht, eine Teilreflexion an der Flache II 
bzw. III zu erwirken. Die Flachen I und II haben, 
wenn der Quarz in der n-ten Oberschwingung an- 
geregt wird, einen Abstand n//2 voneinander. Bei 
der Frequenz, die der Oberschwingung des Quarzes 
genau entspricht, tiberlagern sich die an den Flachen 
I und II reflektierten Wellen mit der gleichen Phase, 
bei einer Verstimmung der Frequenz ist das jedoch 
nicht mehr der Fall. Da Schwingungsimpulse jedoch 
ein Frequenzspektrum endlicher Breite haben, ist 
ihre Verformung unvermeidbar (vergleiche [18]). 
Der Schwingquarz bendtigt als gedampftes mecha- 
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nisches Resonanzsystem eine Einschwingzeit, die um 
so kleiner ist, je besser er an das Medium angepaht 
ist. Bei jeder Reflexion am Quarzende wiederholt 
sich dieser Einschwingvorgang, so da verstandlich 
wird, da die Impulse nach mehreren Reflexionen 
ihre Rechteckform verlieren und eine immer spitzere 
Form annehmen. 

Bei der in Bild 5 dargestellten Echofolge in Quarz- 
glas waren die Bedingungen so gewahlt, daf die 
Impulsverformung sehr gering blieb: der Wellen- 
widerstand des Schwingquarzes ist etwa gleich dem 
des Quarzglases, und der Schwingquarz wurde in 
seiner Grundschwingung (4,8 MHz) angeregt. Zwi- 
schen den planpolierten Flachen des Schwingquarzes 
und der Mefprobe konnte ein homogener Olfilm 
von 0,3 Dicke hergestellt werden. Das Ergebnis 
zeigt Bild 8. Die ersten Echos sind nahezu rechteck- 


Bild 8. Impulsverformung im Quarzglas bei 4,8 MHz. 
a) Direktimpuls und 1. bis 3. Echo, 
b) 20. bis 23. Echo, 
c) §.502 bis ~ 58.; Echo, 
d) 100. bis 103. Echo, 
e) 200. bis 203. Echo, 
f) 300. bis 303. Echo. 


formig, erst vom 50. Echo ab macht sich die Ver- 
formung der Impulse starker bemerkbar. Nach 300 
Reflexionen lassen sich die Echos noch gut von- 
einander trennen. 

Bei den Messungen der Temperaturabhangigkeit 
der Absorption in Metallproben konnte der Olfilm 
nicht so diinn gemacht werden, da er in dieser 
Stérke bei héheren Temperaturen nicht stabil ge- 
wesen wire. Bild 9 zeigt eine fiir die durchgefiihrten 
Messungen charakteristische in einem Aluminium- 


einkristall bei 43,2 MHz erhaltene Echofolge, Bild 10 


ai 


Bild 9. Echofolge im Aluminiumeinkristall (t=25 °C). 
f=43,2 MHz, Lange der Probe 2,8 cm. 
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Bild 10. Auswertung der Echofolge von Bild 9; n:Echo- 
Nr. 


Bild 11. Impulsverformung im Aluminiumeinkristall bei 


43,2 MHz. 

a) 1. bis -3. Echo, 
b) 10. bis 14. Echo, 
c) 20. bis 24. Echo. 
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ihre Auswertung. Die Verformung der Impulse ist 
fiir diesen Fall in Bild 11 gezeigt. 

Um die Verfalschung der Messungen durch Im- 
pulsverformung zu untersuchen, wurde auf die Riick- 
seite einiger Proben ein zweiter Quarz von genau 
der gleichen Dicke angekoppelt. Dadurch verdoppelt 
sich die Zahl der Einschwingvorgange und die Im- 
pulsverformung wird erhdht. Es ergab sich eine 
hohere Mefungenauigkeit, jedoch keine systema- 
tische Abweichung von den Messungen mit gerin- 
gerer Impulsverformung. Messungen mit verschie- 
denen dicken Olfilmen bestatigen dieses Ergebnis. 

Die Impulsverformung ist wie die gesamte Echo- 
charakteristik sehr von der Einstellung der Fre- 
quenz abhangig. Da ein dreifach abstimmbarer 
Sender verwendet wurde, gelang es nur bis zu einem 
gewissen Grade, die optimale Einstellung zu finden, 
womit ein weiterer Teil der MeBungenauigkeit ver- 
ursacht wurde. 


5. MeBergebnisse 


Mit dem beschriebenen Verfahren wurden ver- 
schiedene stranggegossene und stranggepreBte Pro- 
ben aus Aluminium von unterschiedlichem Rein- 
heitsgrad und aus Aluminiumlegierungen untersucht. 
Nur in zwei der betrachteten Proben bildeten sich 
bei Frequenzen unterhalb 20 MHz exponentiell ab- 
klingende Echofolgen aus, wahrend in den ibrigen 
Proben unregelmafige und wesentlich schneller ab- 
klingende Echofolgen auftraten, die auBerdem durch 
starken Storschall zwischen den Echos beeintrachtigt 
wurden. Dieses Verhalten war, wie sich bald heraus- 
stellte, auf die Kornstruktur der Proben zuriick- 
zufthren. 

Die der ausgezeichneten Proben, eine 
AlCuMg-Legierung mit der Zusammensetzung Si 
0,42, Fe 0,33, Cu 3,2, Mn 0,77, Mg 1,28, unter- 
schied sich von den tbrigen durch eine sehr feine 
Kornstruktur. Nach vorangegangener Atzung einer 
Stirnflaiche wurden auf ihr 230 Kristallite pro mm? 
ausgezahlt. Denkt man sich die Korner wiirfelformig 
und alle gleich gro, so wiirden diese Wirfel eine 
Kantenlange von 0,066mm und eine Lange der 
Raumdiagonalen von 0,114 mm haben. Der mitt- 
lere Korndurchmesser liegt zwischen diesen Werten 
und ist hier als das arithmetische Mittel, also 
D=0,09 mm, angenommen. Die durch Streuung an 
den Kristalliten hervorgerufene Extinktion ag, labt 
sich nach einer von Mason [19] angegebenen Theo- 
rie berechnen. Es ist 


eine 


4 J)3 f4 
as, = 8,69 26 OF g SB 
9 4 


v cm 


(15) 


S ist ein fir das verwendete Metall charakteristischer 
Streufaktor, sein Wert ist fir Aluminium S =3-10~4. 
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Fiir f = 20 MHz (D=9-10-* cm, v = 6,26 10° cm/s) 
erhalt man den Wert ag;=0,171 dB/cm, das sind 
25% des gemessenen Wertes (vgl. Bild 12). Fiir 
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Bild12. Absorption in einem Aluminiumeinkristall, 
einer feinkérnigen AlMgCu-Legierung und 
einer Probe aus 99% reinem Aluminium 
(A199) in Abhangigkeit von der Frequenz bei 
Zimmertemperatur. 


tiefere Frequenzen ist der Anteil der Streuung noch 
geringer. Die in dieser Probe gemessene Extinktion 
ist also vorwiegend auf Materialdampfung zuriick- 
zufiihren, wahrend in den nicht naher untersuchten 
grobkérnigen Proben die Streuung die Form der 
Echofolgen maSgeblich bestimmt. 

Die Kornstruktur der anderen ausgezeichneten 
Probe, die aus 99% reinem Aluminium (Al 99) der 


Bild13. Durch Atzen sichtbar gemachte Struktur der 
zwei Stirnflachen einer Aluminiumprobe (Al 99). 
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Zusammensetzung Si 0,29, Fe 0,47, Cu 0,034, Mn 
0,035, Mg 0,015, Zn 0,047 bestand, zeigt Bild 13. 
Auf beiden Stirnflachen ist jeweils eine grofere 
homogene Flache sichtbar. Es sind die Enden eines 
Einkristalls, der sich durch die Probe in Langsrich- 
tung hindurchzieht. Nahere Untersuchungen ergaben, 
dafi die Schalliibertragung fast ausschlieBlich durch 
diesen Einkristall erfolgte. 
SchlieBlich wurden Messungen in einem Einkri- 
stall aus 99,98% reinem Aluminium durchgefihrt, 
der in Bild 3 gezeigt ist. Die Achse der zylindrischen 
Probe ist die 1,0,0-Achse des kubisch flachenzen- 
trierten Kristalls. 
Wie aus Bild 12 hervorgeht, wachst die Material- 
_dampfung in allen Proben quadratisch mit der Fre- 
quenz. Dabei ist die Dampfung in der polykristal- 
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Bild14. Absorption in einer feinkérnigen AlCuMg- 
Legierung in Abhangigkeit von der Tempe- 


ratur. 
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Bild 15. Absorption in 99% reinem Aluminium mit der 
in Bild 13 gezeigten Struktur in Abhangigkeit 
von der Temperatur. 


linen Legierung mehr als 10mal so grof} wie die 
im Einkristall. Da die Schallausbreitung in der 
Probe aus Al 99 undefiniert ist, streuen hier die 
MeBpunkte starker als bei den andéren Proben. 
Dort ist die MeBgenauigkeit auf Grund der oben 
beschriebenen Effekte =10%, mindestens aber 
+0,01 dB/cm. 

Die Temperaturabhangigkeit der drei Proben ist 
in den Bildern 14 bis 16 dargestellt. In der poly- 
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Bild 16. Absorption im Aluminiumeinkristall in Abhan- 
gigkeit von der Temperatur bei verschiedenen 
Frequenzen. 


kristallinen Legierung findet man bei Temperaturen 
oberhalb 300°C einen plotzlichen Anstieg der 
Dampfung, der aber nicht so stark ist, wie der von 
Zucxer [10] in einer ahnlichen Legierung gefun- 
dene. Im Einkristall tritt dieser plotzliche Anstieg 
nicht auf, man findet vielmehr eine kontinuierliche 
Zunahme der Dampfung mit steigender Temperatur. 
In Bild 17 ist die Absorption im Einkristall als 
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Bild 17. Absorption im Aluminiumeinkristall als Funk- 
tion der reziproken absoluten Temperatur. 
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Funktion der reziproken absoluten Temperatur auf- 
getragen. Dieser Abbildung entnimmt man, daB die 
Absorption im betrachteten Frequenzbereich durch 
die Beziehung 

a REE eee | 


Pp = 6)35 104 oes RT 


(16) 
cm 

922 cal/mol 

Tea 


s?cm ! 


= ol: LO ven 


dargestellt werden kann. Fiir Zimmertemperatur 


(T =300 °K) folgt daraus 
a)? = 56 a0 es cmiewe 


Ein Vergleich zwischen diesem Ergebnis und dem 
von Ké [6], [19] mit Torsionsschwingungen bei der 
Frequenz 1 Hz in einem Aluminiumeinkristall gefun- 
denen Zusammenhang zeigt, da der fiir Zimmer- 
temperatur sich ergebende Wert 2=3,2:10 °cm™“1 
weit oberhalb des Wertes liegt, der sich durch Extra- 
polation der gefundenen Frequenzabhangigkeit auf 
1 Hz ergibt. Das bedeutet, da bei tiefen Frequen- 
zen ein zusatzlicher Dampfungsmechanismus mit 

einer tiefliegenden Relaxationsfrequenz auftritt. Wie 

zu erwarten, besteht daher auch keine Ubereinstim- 
mung der in den verschiedenen Frequenzbereichen 
gefundenen Aktivierungsenergien. Bei tiefen Fre- 
quenzen wurde von KE H =5200 cal/mol gefunden, 
hier ergab sich H = 922 cal/mol. 


Mein Dank gilt Herrn Prof. Dr. Dr.-Ing. E.h. 
E. Meyer fur die Anregung zu dieser Arbeit. [hm 
und den Herren Prof. Dr. K. Tamm, Prof. Dr. 
P. Haasen und Dr. M. Evers danke ich fiir wert- 


volle Diskussionen. Die Deutsche Forschungsgemein- 
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schaft unterstiitzte die Arbeit in dankenswerter Weise 
durch eine Sachbeihilfe. — (Eingegangen am 15. Oktober 1960.) 
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UNTERSUCHUNGEN ZUR WECHSELWIRKUNG 
ZWISCHEN SCHALLFELD UND STROMUNG IM WINDKANAL, 
DER MIT POROSEN ABSORBERN UND 
HELMHOLTZRESONATOREN AUSGEKLEIDET IST 


von W. Scuriz 


III. Physikalisches Institut der Universitat Gottingen 


Zusammenfassung 


Die Schalldampfung in absorbierend ausgekleideten Windkanalen zeigt eine starke 
Abhangigkeit von der Richtung und Gré8e der Gleichstrémung. Die Dampfungsanderung 
beruht auf einer Kennwertanderung des absorbierenden Materials und einer Verformung 
der Phasenflachen des Schallfeldes durch die Gleichstromung. Die unter Berticksichtigung 
dieser beiden Effekte berechnete Dampfungsanderung ergibt eine gute Ubereinstimmung 
mit den MeSergebnissen. 

Der Mechanismus der Schallverstarkung in einem Kanal, der mit unbedampften Helm- 

- holtzresonatoren ausgekleidet ist, wird aus der Wechselwirkung zwischen Schallfeld und 
Stromungsturbulenz erklart. Die Umwandlung yon Stromungsenergie in Schallenergie wird 
durch die Synchronisation der Turbulenz bewirkt, die bei geeigneten Werten von Stromungs- 
geschwindigkeit und Schallfrequenz zur Ausbildung einer stabilen Pseudoschallwelle fiihrt. 
Die Kopplung zwischen Pseudoschall und Schall geschieht an den Resonatorhalsen. 


Sommaire 


L’atténuation du son dans des conduits de vent revétus de matériaux absorbants dépend 
trés fortement de la direction et de l’intensité du flux continu du vent. Les variations de 
Vamortissement s’expliquent par les modifications des caractéristiques du matériel absor- 
bant et par la déformation des surfaces de phase du champ sonore produites par le flux 
continu de vent. En tenant compte de ces deux effets, le calcul de la modification d’amor- 
tissement est en bon accord avec les résultats de ’expérience. Woe 

Le mécanisme de l’amplification sonore dans un canal revétu de résonateurs d’Helmholtz : 
non amortis s explique par l’action réciproque du champ sonore et de la turbulence du cou- 
rant (d’air). La transformation de l’énergie du courant en énergie sonore est due a la syn- 
chronisation de la turbulence qui, pour certaines valeurs particuliéres de la vitesse du cou- 
rant et de la fréquence du son, conduit a un état stable d’ondes pseudo-sonores. Le couplage : 
entre le son et le pseudo-son se fait aux cols des résonateurs. t. 


Summary 


Sound attenuation in ducts with absorbent coating of the walls is strongly dependent on 
the direction and velocity of the air flow. The change of attenuation is brought about by 
a change of the characteristic properties of the absorbing material and by the deformation ‘ 
of the phase surface of the sound field caused by the air flow. Calculated values of the 
attenuation taking both effects into account are in good agreement with experimental values. F) 

The sound amplification mechanism in a duct coated with undamped Helmholtz resona- : 
tors is explained in the interaction between the sound field and the turbulence of flow. Con- 
version of flow energy into sound energy is effected by synchronisation of the turbulence. “| 
Under favourable circumstances this synchronisation leads to the formation of a stable iE 
pseudo sound wave. Pseudo sound and sound interact at the resonator necks. 


1. Einleitung 


Untersuchungen iiber die Schalldampfung in ab- 
sorbierend ausgekleideten Windkanalen wurden 
unter anderen von F. Mecuex [1] durchgefiihrt. Sie 
erstreckten sich auf die Abhangigkeit des Frequenz- 
ganges der Dampfung von der Grofe der iiberlager- 
ten Windstromung. Als Dampfungselemente wurden 
im wesentlichen pordse Absorber (Gesteinswolle) 
und Helmholtzresonatoren verwendet. Bei pordsen 
Absorbern zeigte sich eine starke Anderung der 


Dampfung mit wachsender Stromungsgeschwindig- 
keit und eine auffallende Unsymmetrie zwischen 
Schallausbreitung mit und entgegen der Stromungs- 
richtung. Die Dampfungsanderung an pordsen Ab- 
sorbern, die nach Beriicksichtigung der durch Uber- 
lagerung von Schallgeschwindigkeit und Strémungs- 
geschwindigkeit bedingten Wellenlangenanderung 
erhalten bleibt, ist in Bild 1 wiedergegeben. 

Der mit Helmholtzresonatoren abgedeckte Kanal 
weist oberhalb der Resonanzstelle, die sich durch 
sehr hohe Dampfung auszeichnet, einen Frequenz- 
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Bild 1. Anderung der auf die Wellenlinge bezogenen 
Kanaldampfung als Funktion der Windgeschwin- 
digkeit. 


bereich auf, in dem das eingespeiste Schallsignal ver- 
starkt wird. Dabei sind die Frequenz maximaler 
Verstaérkung 7) und die dabei vorliegende Stro- 
mungsgeschwindigkeit durch die Dispersionsformel 


V = (-1+V1+ (4L/con) ») 
miteinander verknupft. 


Dampfungskurven bei verschiedenen Wind- 
geschwindigkeiten sind in Bild 2 wiedergegeben. 
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Bild 2. Beispiele von Frequenzkurven der Dampfung 

bei Helmholtzresonatoren. Die Resonanzfrequenz 
betragt 800 Hz. 


Ziel der vorliegenden Arbeit ist es, durch genaue 
Untersuchung des Schallfeldes und der Turbulenz- 
und Strémungsverteilung im Innern des Kanals die 
Dampfungsainderung bei porésen Absorbern und 
den Mechanismus der Schallverstarkung bei reakti- 
ven Absorbern aufzuklaren. ¥ 


2. Messungen im Kundtschen Rohr 
2.1. Porése Absorber 


Die Anderung der Kanaldimpfung kann auf ver- 
schiedene Weise zustande kommen. Als erster Punkt 
soll der Einflu8 der Strémungsturbulenz auf die 
Impedanz der Steinwolleschicht betrachtet werden. 
Wichtig hierftir ist der hohe Turbulenzgrad der be- 
nutzten Kanalstrémung. Die Mittenturbulenz betragt 
etwa 5%, doch steigen Quer- und Lingskomponente 
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der Turbulenz an den Kanalwanden stark an. Dort 
kann die normal auf der Absorberoberfliche ste- 
hende Turbulenzkomponente in den Absorber ein- 
dringen. Da der Stromungswiderstand poréser Ab- 
sorber, der mitbestimmend fiir die Impedanz ist, nur 
bei kleinen Schallschnellen linear ist, wird durch die 
groBen Wechselgeschwindigkeiten der Turbulenz der 
Stromungswiderstand in den nichtlinearen Bereich 
ausgesteuert. Wegen dieser Nichtlinearitat wirkt 
auch auf eine tiberlagerte kleine Schallschnelle der 
erhohte Strémungswiderstand. Die mittlere Turbu- 
lenzschnelle ist eine Funktion der Windgeschwindig- 
keit; somit la8t sich der Stromungswiderstand r des 
Absorbers als Funktion der Windgeschwindigkeit V 
schreiben: 
(1) 
Mit steigendem Stromungswiderstand der Wand fallt 
die Kanaldampfung, da der Absorber schlechter an 
das freie Medium angepafit ist. Bei einem Kanal, 
dessen Querdimensionen klein zur Wellenlange sind, 
laBt sich aus der Wandimpedanz 

[Z| = eyVo/m Vro tr (V) (2) 
die folgende Naherung fiir die Kanaldampfung an- 
geben: 


r=rot+r(V). 


a~[rotr(V)]~. (3) 


Die Naherungsbedingungen sind bei dem Kanal, der 
zu diesen Untersuchungen benutzt wurde, erfillt. 

Die in Bild 3 dargestellte Anordnung dient zur 
Messung der Impedanz eines pordsen Absorbers 
bei Uberlagerung einer Gleichstr6mung. Das Kundt- 
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Bild 3. MeSkammer am Kundtschen Rohr zur Bestim- 
mung der Impedanz poréser Stoffe bei Wind- 
iiberlagerung. 


sche Rohr wird mit einer Kassette abgeschlossen, die 
mehrere Schichten Gesteinswolle enthalt. Durch 
diinne Bohrungen in zwei gegentiberliegenden Sei- 
tenwanden wird Luft eingepreBt bzw. abgesaugt. Die 
stark bedampften konusformigen Ansatze gewahr- 
leisten eine hinreichende akustische Entkopplung von 
Me8kammer und Luftzufuhr. 

Mit dieser Anordnung wurde die Frequenzkurve 
der Impedanz bei verschiedenen Stromungsgeschwin- 
digkeiten gemessen. Ein Beispiel einer solchen Mes- 
sung ist in Bild 4 in der komplexen Widerstands- 
ebene dargestellt. Die durch die Dimensionen der 
Me8kammer bestimmte Resonanzschleife zieht sich 
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Bild 4. Frequenzgang der Impedanz der Probe bei ver- 
schiedenen Windgeschwindigkeiten. 


mit wachsender  Gleichstromungsgeschwindigkeit 


stark zusammen, das heift der Stromungswiderstand 
und damit die Dampfung werden erhoht. Dabei 
bleibt die Frequenzlage der Resonanz nahezu erhal- 
ten. Die aus diesen MeSkurven entnommene Ande- 
rung des Stromungswiderstandes zeigt Bild 5, 
Kurve A. Bei Verwendung von sehr diinnen Gesteins- 
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Bild5. Gemessene Anderung des Strémungswider- 
standes r. 


wolleschichten (etwa 1mm) ist die Anderung des 
Stromungswiderstandes grofer, da hier eine starkere 
Durchmischung stattfindet (KurveB). Wegen der 
nur sehr geringen Anderung des Imaginarteiles der 
Impedanz (und damit der Resonanzfrequenz der 
Kassette) werden die spater am Kanal ausgefiihrten 
Impedanzmessungen auf den Stromungswiderstand 


beschrankt. 


2.2. Bedimpfte Helmholtzresonatoren 


Zur Erganzung dieser Messungen wurde die Im- 
pedanzanderung an einem bedampften Helmholtz- 
resonator gemessen. Der Hals des Resonators war 
mit einer etwa 0,1 mm starken Gesteinswolleschicht 
abgedeckt, ittber die Luft geblasen wurde. Um die 
Stoérung der Messung durch die Luftzufuhrstutzen 
auszuschalten, wurden diese als Kurzschlufischieber 
ausgebildet, deren harter Abschlufi fiir jede Fre- 
quenz an die Eintrittsoffnung in das Kundtsche Rohr 
transformiert werden konnte. 

Die aus der Halbwertsbreite der Resonanzkurven 
bestimmte Dampfung zeigt ebenfalls einen Anstieg 
mit der Strémungsgeschwindigkeit, der in Verlauf 
und Gréffe den Messungen am pordsen Absorber 
ahnlich ist. Eine merkbare Verschiebung der Reso- 
nanzfrequenz tritt nicht auf. Dies ist auch bei den 
relativ kleinen Strémungsgeschwindigkeiten (etwa 
5 m/s) nicht zu erwarten. 

Bei stark bedampften Helmholtzresonatoren und 
bei pordsen Absorbern ist der Mechanismus der 


Dampfungsanderung sehr ahnlich; demzufolge zei- 
gen auch die von F. Mecuet [1] gemessenen Kurven 
der Dampfungsanderung fiir beide Absorbertypen 
den gleichen Verlauf. 


3. Schalldampfung im Windkanal mit einer 
schallabsorbierenden Wand aus Gesteinswolle 


Der bei den folgenden Messungen benutzte Wind- 
kanal hat eine Lange von 2,50m bei einem freien 
Querschnitt von 35mmx100mm. Die maximale 
Windgeschwindigkeit liegt bei 80 m/s. Die Wind- 
stromung wird von einem |-stufigen Radial-Geblase 
geliefert und in zwei Dampfern von insgesamt 4m 
Lange akustisch beruhigt. Dadurch hat die Stro- 
mung am Kanalanfang einen niedrigen Storpegel. 
Zur Uberlagerung eines Schallsignals sind zwei 
Lautsprecher von je 200 W elektrischer Leistung 
am Kanalanfang angekoppelt, die eine verfiigbare 
Dynamik von etwa 40 dB gewahrleisten. 


3.1. Verkniipfung von Impedanzénderung und 
Kanaldémpfung 


Um die in Abschnitt 2.1 besprochene Impedanz- 
anderung auf die Kanaldampfung zu itbertragen, 
mufte eine genaue Verkniipfung dieser Anderung 
mit der Mittengeschwindigkeit im Kanal gefunden 
werden. Hierzu geniigte, wie schon gezeigt, eine Mes- 
sung des Stromungswiderstandes. Bild 6 zeigt eine 
Skizze der MeSanordnung. Das Messingrohr von 
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Bild 6. MeBanordnung am Windkanal zur Bestimmung 
des _geschwindigkeitsabhangigen Stromungs- 
widerstandes. 


etwa 15 cm Durchmesser dient zur Abschirmung der 
im Schluckstoff verlaufenden Welle. Die Drucksonde 
mit einem inneren Durchmesser von 3 mm mifbt den 
Abfall der normal zur Absorberoberflache verlau- 
fenden, stark gedampften Welle. 

Ebenso wie das Schallsignal konnen auch die mit 
der Turbulenz verkniipften schwachen Druckande- 
rungen aufgenommen werden. Im Absorber verhal- 
ten sich Schall und Turbulenz gleich, das heifit die 
Messungen ergeben fiir beide die gleiche Damp- 
fungskonstante. Bei Uberlagerung einer Windstré- 
mung steigt die Dampfung mit wachsender Ge- 
schwindigkeit stark an, ist allerdings iiber die Tiefe 
des Absorbers nicht konstant. Da die Turbulenz mit 
wachsendem Abstand von der kanalseitigen Absor- 
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beroberflache stark abnimmt, verringert sich auch 
die Aussteuerung des Stromungswiderstandes und 
damit die Schalldampfung. Ein Beispiel fiir diese 
Ortsabhangigkeit der Dampfungskonstante zeigt 
Bild 7. 
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Bild 7. Ortsabhangigkeit der Dampfung im Schluckstoff 
(vergleiche hierzu Bild 6). 


Die weiter oben angegebene Naherungsrechnung 
fiir die Kanaldimpfung (Gl. (2) und (3)) enthalt 
die Impedanz eines homogenen Schluckstoffes. Durch 
die stromungsbedingte Ortsabhangigkeit des Stro- 
mungswiderstandes normal zur Kanaloberflache liegt 
aber bei Windiberlagerung ein inhomogener Stoff 
vor. Die genaherte Berechnung des Reflexionsfaktors 
der Absorberschicht ergab eine Abweichung von 
etwa 4%, wenn die Anderung des Stromungswider- 
standes im Absorber beriicksichtigt wurde. Diese 
Abweichung wird bei der Auswertung durch Ver- 
wendung eines geringfiigig erniedrigten Stromungs- 
widerstandswertes im Rahmen der Mefgenauigkeit 
kompensiert. 

Mit dieser MeBanordnung wurde nun fur die Fre- 
quenzen 0,4, 0,6, 1,0 und 1,4 kHz die relative Ande- 
rung des Stromungswiderstandes r/ry bestimmt. Der 
Verlauf ist fiir alle Frequenzen verschieden, doch 
laBt sich mit Hilfe eines frequenzabhangigen Para- 
meters A, eine empirische Zusammenfassung der 


MeSwerte durch die Gleichung 
rire Ae je (4) 


erreichen. Diese Gleichung und die Mefwerte sind 
in Bild 8 dargestellt. Die Abszisse ist in A,V geteilt, 
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Bild 8. Relative Anderung des Strémungswiderstandes. 


damit die Kurven zusammenfallen. Die empirischen 
Werte fiir A, sind 


Ao 4 = 9: home s/m, 
Aog = 3,5'1073 s/m, 


Ai = 2,5-1078 s/m, 
Ay,4=2,3-1073 s/m. 


Die Werte fiir r/ry wachsen nicht iiber 3,5 an. Hier- 
aus ist die Abweichung der Messungen bei 0,4 kHz, 
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deren A, auch stark von den ibrigen verschieden 
ist, als ein Sattigungseffekt zu deuten. Wegen des 
glatten Verlaufes der Funktion A,={(v) laBt sich 
fiir andere Frequenzen die relative Anderung des 
Strémungswiderstandes leicht interpolieren. 

Berucksichtigt man, da8 mit der Erhohung des 
Stromungswiderstandes im wesentlichen nur der 
Realteil der Impedanz des Absorbers geandert wird, 
so ]aBt sich mit Hilfe von Gl. (3) der Str6mungs- 
widerstand mit der Kanaldémpfung verkniipfen. 
Auch hier ist wieder wichtig, da die Kanalhohe 
klein zur Wellenlange bleibt. Mit Gl. (4) zusammen 
erhalt man fiir die Dampfungskonstante a der 
Kanaldimpfung: 


ane A 4 Vy (5) 


Bei den oben angefiihrten Messungen ist das Rich- 
tungsverhaltnis von Wind zu Schall ohne Einfluf. 
Unterscheidet man nun die Schallausbreitung mit, 
beziehungsweise gegen die Stromungsrichtung durch 
das Vorzeichen von J, so ist in Gl. (5) V durch | V| 
zu ersetzen. Die auf diese Weise berechnete Ande- 
rung der Kanaldampfung, die also durch die Impe- 
danzerhohung der Gesteinswolle hervorgerufen wird, 
ist in Bild 9 als Funktion der Frequenz eingetragen. 


— - - —+— 07 + + + 204 
-80 -60 -40 -20.0 20 40 60 80m/s 


y— 
Bild 9. Relative Anderung der Kanaldampfung, die 
durch die Kennwertanderung bedingt ist. 


Dampfungsanderung —* 


“BOSD, ieee 20 40 60 80mis 
—> 
Bild 10. Reduktion der Dampfungskurven yon Bild 1 
auf konstante Wandimpedanz bei allen Ge- 
schwindigkeiten. 


Mit diesen Werten sind die in Bild 1 gezeigten 
Dampfungskurven korrigiert worden und in Bild 10 
dargestellt. 

Nachdem nun der Einflu8 der Schallwellenlangen- 
anderung und der Impedanzerhéhung beriicksichtigt 
ist, zeigen die Dampfungskurven noch eine Abhan- 


- gigkeit von der Windgeschwindigkeit. Wahrend bei 


Vorwartsausbreitung die so korrigierte Dampfung 
innerhalb der Fehlergrenze als konstant zu betrach- 
ten ist, zeigt sich bei Riickwartsausbreitung ein star- 
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ker Anstieg der Dampfung mit der Windgeschwin- 
digkeit. Dieses Verhalten soll nun im folgenden aus 
- der Wechselwirkung zwischen den Phasenflachen des 
Schallfeldes und dem Profil der Kanalstromung er- 
klart werden. 


3.2. Geschwindigkeitsverteilung und Phasenflachen 
im Kanal 


Zur Erfassung dieser Wechselwirkung ist zunachst 
_ eine genaue Kenntnis der Stromungsverteilung not- 
wendig. Ein Hitzdraht-Anemometer, das wegen sei- 
ner Kleinheit das Stromungsfeld nur wenig stort, 
dient zur Ausmessung der Geschwindigkeitsvertei- 
lung. Gleichzeitig konnen Stromungsfluktuationen 
im Frequenzgebiet von 2 bis 100 000 Hz aufgenom- 
men und phasen- sowie amplitudenkorrigiert wie- 
dergegeben werden. Die Mefanordnung ist absolut 
eichbar. 

Durch die verschiedene Struktur der Kanalwande 
— die Oberseite besteht aus Schluckstoff, der Boden 
ist poliertes Holz — ist das Stromungsprofil unsym- 
metrisch. Der Schwerpunkt der Strémung fallt nicht 
mit der Kanalachse zusammen. 

Eine einfache Beurteilung des Stromungsprofils 
ist mit Hilfe des universellen Geschwindigkeits- 
gesetzes der Rohrstromung [2] moglich: 


(u/U)* = (Y/R). (6) 


Darin ist u die lokale Geschwindigkeit, Y der Wand- 
abstand, U die Mittengeschwindigkeit und 2R die 
Kanalhohe. 

Aus der doppelt-logarithmischen Auftragung des 
Stromungsprofils bestimmt man den Profilexponen- 
ten K. Die so gefundenen K-Werte liegen im Bereich 
der Kanalstrémung bei glatten Wanden. Dies war 
_ zu erwarten, da die Oberflachenrauhigkeit der Ge- 
steinswolle klein zur Grenzschichtdicke ist, die im 
Mittel bei 3 bis 4mm lag. Die gemessenen K-Werte 
sind fiir Absorberseite und Boden des Kanals ver- 
schieden. Wahrend die am Boden gemessenen Werte 
mit den Soll-Werten bei entsprechenden Reynolds- 
zahlen (auf die Kanalhohe bezogen) wbereinstim- 
men, zeigen die absorberseitigen Profile starke Ab- 


weichungen (Tab. I). 


Tabelle I. 
U Peeks. K eae a 
“ee Ga eeabe 1, |2. “Boll. 
11 7,3 - 104 5,26 | 6,8 
45 3,0 - 10° GBT eo 17,7 
55 3,6 - 105 Fad iets 07,8 


Die gemessenen K-Werte liegen stets tiefer als 
die Soll-Werte. Es liegt nahe, als Begriindung hier- 
fiir eine Turbulenzerniedrigung der Kanalstromung 
durch die absorbierende Berandung anzusehen. Da 
der normal zur Absorberoberflache schwingende 


Turbulenzanteil genau wie die Schallschnelle in den 
Absorber eindringt und gedampft wird, ist auf die- 
sem Weg ein Entzug von Turbulenzenergie moglich. 

Verlauft eine Schallwelle in einem bewegten Me- 
dium der Geschwindigkeit V parallel (antiparallel) 
zur Stromungsrichtung, so erhalt man mit der Schall- 
geschwindigkeit cy in ruhender Luft eine resultie- 
rende Phasengeschwindigkeit 


Ca Cg (7) 


Die in Bild 11 gezeichneten Stromungsprofile sind 
sehr vollig, so da® sich als gute Naherung mit der 
Mittengeschwindigkeit V,, eine mittlere Phasen- 
geschwindigkeit entsprechend der Gl. (7) bilden 
laBt. Fur die genaue Untersuchung des Schallfeldes 


‘ Absorbergrenze 
Gy ° 11 m/s 
+45 m/s 
| rv) 955 m/s 
mm 
“ig L 1 4 
06 07 08 09 40 
U/ Vin —— 


Bild 11. Stromungsprofile bei verschiedenen Wind- 
geschwindigkeiten. 


quer zur Ausbreitungsrichtung ist es jedoch not- 
wendig, eine von der Kanalhohe y und der Mitten- 
geschwindigkeit V,, abhangige Phasengeschwindig- 
keit zu definieren: 


c(y, Vm) =cot+V (y, Vm). (8) 


Diese Phasengeschwindigkeitsprofile sind fur ver- 
schiedene Windgeschwindigkeiten in Bild 12 dar- 


(Co + icq 
1404 408 112 61,20 
T T T 


0 aT Ley Ee 
2 Absorbergrenze 
Ne 
y 7 x ° 11 m/s 
\ + 45 m/s 
3 55 m/s 
12 4 
mmLi \ 
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Bild 12. Schallgeschwindigkeitsprofile bei Vorwarts- und 
Riickwartsausbreitung. 


gestellt. Hierbei bezieht sich der obere (untere) 
Abszissenmafistab auf Vorwarts-(Rickwarts-) Aus- 
breitung. 

Diese -Phasengeschwindigkeitsprofile bestimmen 
die Verformung der Phasenflachen des Schallfeldes. 
Wahrend im schallharten Kanal ohne Windiber- 
lagerung unterhalb der Grenzfrequenz eine ebene 
Phasenflache vorliegt, die normal zu den Wanden 
verlauft, tritt bei absorbierender Berandung eine 
Kriimmung der Phasenflache gegen die Wand der- 
art auf, daf} die Phase an der Berandung nacheilt. 


er ee a 
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Diese Verformung der Phasenflache ist das Er- 
gebnis der Wechselwirkung zwischen dem Schallfeld 
im Kanal und den Nahfeldern in der Absorberober- 
flache. Der Aufsetzwinkel der Phasenflache auf die 
Berandung ist mit den Kennwerten des Absorbers 
durch die folgende Gleichung verkniipft [3] : 


tan y= y 2e Wel (9) 
o Hy H+) 


(o Porositat, ¢ =@/0) Dichteanderung im Absorber, 
y =r/o @ Verlustfaktor). 

Mit wachsendem Winkel wy fallt die Kanaldamp- 
fung. 

Durch die Kennwerte des Absorbers ist nun am 
Anfang der Dampfungsstrecke ein Phasenprofil 
Poly) vorgegeben, das mit dem Phasengeschwindig- 
keitsprofil c(y, Vm) in Wechselwirkung tritt. Man 
erhalt daraus ein Phasenprofil, das sich mit wach- 
sendem Laufweg kontinuierlich verandert. Fir Orte 
gleicher Phase gilt dann nach einem Laufweg von 


Ly V 


P=2H die Gleichung 


Cm 


x(p,y) i [yp — Poly) ] [ey + V ty, Vin) | wer, (10) 


In Bild 13 ist diese Gleichung ftir @=22 600 Hz 


Absorber 
Stromung 
—---S 
Schall 

15cm 30cm 60cm 

Absorber 
Stromung 
ene 
Schall 

10cm \ 20cm 42cm 


Bild 13. Berechnete Phasenprofile bei Vorwdarts- und 
Riickwartsausbreitung nach verschiedenen Lauf- 
wegen. 


und V,,=45 m/s bei verschiedenen Laufwegen aus- 
gewertet. Bei Vorwartsausbreitung pragt sich das 
Stromungsprofil dem Phasenprofil immer mehr auf, 
wobei jedoch der Aufsetzwinkel zwischen Phasen- 
flache und Absorberflache nur wenig verandert wird. 

Wesentlich stirker ist der Einflu8{ der Windstro- 
mung bei Schallausbreitung entgegen der Stromungs- 
richtung. Hier ergibt die Rechnung bei kurzem Lauf- 
weg eine Aufhebung der Kriimmung der Phasen- 
flache. Nach einer Strecke von etwa 20cm erhalt 
man jedoch ein Phasenprofil, dessen Kriimmung 
umgekehrt ist. Die wandnahen Anteile des Schallfel- 
des zeigen eine vorlaufende Phase. Diese erhebliche 
Veranderung des Phasenprofils bei Riickwartsaus- 
breitung im Gegensatz zur Vorwartsausbreitung 
steht in Einklang mit den Mefkurven aus Bild 10, 
die auch nur bei Riickwartsausbreitung noch eine 
starke Abhangigkeit der Dampfung von der Wind- 
geschwindigkeit zeigen. 
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3.3. Messung der Phasenfldchen 


Die nun folgende experimentelle Untersuchung 
der Phasenflache wurde im wesentlichen mit einem 
Schallsignal von 600 Hz ausgefihrt. Bei dieser Fre- 
quenz ist die Wellenlange sehr gro8 gegeniiber der 
Kanalhohe, und die verfiigbare Dynamik der gesam- 
ten MefSiapparatur gewdhrleistet einen Stérabstand, 
der fiir das einwandfreie Arbeiten des Phasenmes- 
sers notwendig ist. Die vorwiegende Benutzung einer 
Frequenz ist keine Einschrankung, da der Damp- 
fungsverlauf fiir alle Frequenzen im_betrachteten 
Bereich ahnlich ist. 


Zur Ausmessung der Phasenflachen dient eine 


Schalldrucksonde. Thre Ausfithrung zeigt Bild 14. 


Mikrofon 


Bild 14. T-formige Sonde zur Schalldruckmessung. 


Der horizontale Teil ist stromlinienformig ausgebil- 
det und tragt in der Mitte die mit Gaze abgedeckte 
Schalleintrittséffnung. Diese Offnung liegt etwa 
5 Sondendurchmesser hinter der Sondenspitze. An 
dieser Stelle ist das von der Sonde erzeugte Stro- 
mungsrauschen am kleinsten, und damit die verfiig- 
bare Dynamik grofer. 


Die Verwendung eines Schalldruckempfangers an 
Stelle des Schnelleempfangers (Hitzdraht) ist eben- 
falls aus Griinden der Dynamik geschehen. Das Ver- 
haltnis von Druck zu Schnelle ist bei Turbulenz um 
Grofenordnungen kleiner als bei Schall, so dafi man 
bei gleicher mittlerer Schnelle von Schallfeld und 
Turbulenz mit dem Druckempfanger noch ein relativ 
ungestortes Schallsignal aufnehmen kann. Bei der 
Besprechung des Pseudo-Schalls wird darauf noch 
genauer eingegangen werden (4.2). Diese Druck- 
sonde wird normal zur Absorberoberflache tiber die 
Kanalhéhe bewegt und der Verlauf des Phasenwin- 
kels gemessen. Diese Winkelanderungen sind in 
Ortsanderungen umgerechnet, so da die Kurven 
Punkte gleicher Phase verbinden, wie sie auch in 
Gl. (10) berechnet wurden. 


Bild 15 zeigt bei Vorwartsausbreitung gemessene 
Phasenprofile fiir verschiedene Windgeschwindig- 
keiten. Das Phasenprofil V,, =0 m/s wurde schon in 
Gl. (10) als Anfangsphase 9) (y) benutzt. Den Ver- 
lauf des aus diesen Kurven entnommenen Aufsetz- 
winkels ¢’ zeigt Tabelle II. 


Tabelle IT. 
Vm in m/s PPA ae 25 45 | 55 
ee ee el eta ee 40° | 479 
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Wahrend nach der Rechnung eine Verkleinerung 
des Aufsetzwinkels zu erwarten war, ergibt die Mes- 
sung eindeutig eine Zunahme des Winkels y’ mit 
wachsender. Windgeschwindigkeit, das heift eine 
Abnahme der Kanaldampfung. Bei der Beurteilung 
ist jedoch zu beachten, dafi sich die Einfltisse des 
Strémungsprofils und der Anderung der Absorber- 
Kennwerte bei der Messung nicht trennen lassen. 
Die in 3.1. behandelte Zunahme des Stromungs- 
_ widerstandes im Absorber hat auch eine Aufsteilung 
der Phasenflache gegentiber ruhender Luft zur Folge. 
Die beiden Ursachen fiir die Kriimmung der Pha- 
senflache, Kennwertanderung und Stromungsprofil, 
sind in ihrer Wirkung gegenlaufig, so da} bei dem 
geringen Einflu8 der Stromung auf das Phasenprofil 
bei Vorwartsausbreitung die VergroBerung des Auf- 
setzwinkels durch die Kennwertanderung tiberwiegt. 
Dieses Ergebnis wird durch den Verlauf der Kurven 
in Bild 10 bestatigt, die zeigen, daf} bei Vorwarts- 
ausbreitung schon die Korrektur der Kennwerte zu 
konstanter Damptfung fuhrt. 

Die gemessenen Phasenflachen bei Rtickwartsaus- 
breitung stehen in gutem Einklang mit dem Ergeb- 
nis der Rechnung. Wegen der starken Verformung 
wurde die Anderung der Phasenflachen nach ver- 
schiedenen Laufwegen gemessen. Das Profil nach 
20cm Laufweg ist fur die Windgeschwindigkeiten 
11,25 und 55m/s in Bild 16 dargestellt. Eine ein- 
schneidende Verformung tritt erst bei 55 m/s auf. 
Durch die gegenlaufige Krummung von Phasenflache 
und Stromungsprofil erhalt das Phasenprofil eine 
doppelte Krimmung. Voll ausgebildet sind die Pha- 
senflachen nach einem Laufweg von etwa 40 cm. Wie 
aus Bild 17 ersichtlich ist, unterscheiden sich die 
Profile fir alle Geschwindigkeiten stark und zeigen 
keine Ahnlichkeit mehr mit dem ungestorten Phasen- 


Absorber 
: Stro 
20 seus 
30 
Bite 0 "1 25 45 55 m/s 


Bild 15. Gemessene Phasenprofile bei Vorwartsausbrei- 


tung. 
0 Absorber 
10 Stromung 
20 Schall 
ese Aten 25-42 -55.4m/s 


Bild 16. Gemessene Phasenprofile bei Riickwartsaus- 
breitung nach kurzem Laufweg (etwa 20 cm). 
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Bild 17. Voll ausgebildete Phasenprofile bei Riickwarts- 
ausbreitung (Laufweg > 40cm). 
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profil vor einer absorbierenden Wand. Die in 
Gl. (9) angegebene Beziehung zwischen Aufsetz- 
winkel und Absorber-Kennwert, aus der die Kanal- 
dampfung bestimmt werden kann, ist hier nicht 
mehr giltig. 


3.4. Berechnung der Dimpfungsdnderung 


Eine annahernde Berechnung der Kanaldampfung 
laBt sich aus der Definition der Dampfung « = dN/N 
herleiten. Dabei ist dV die absorbierte Schalleistung 
und WN die Schalleistung in Ausbreitungsrichtung. 
Durch die stromungsbedingte Verformung der Pha- 
senflache wird die stromungsparallele Schalleistung 
N geandert. Als MaB fiir diese Anderung dient der 
wandnormale Energieanteil des Schallfeldes E, , der 
dem Zeitintegral der absorbierten Schalleistung und 
damit der Abnahme der stromungsparallelen Ener- 
gie entspricht und aus der Form der Phasenflache 
entnommen wird. Bei dieser Rechnung wird die von 
Risner [4] angegebene Gleichung benutzt, die den 
stromungsnormalen Anteil eines Schallfeldes be- 
stimmt, dessen Wellenfronten unter dem Winkel a zu 
einer Stromung der Geschwindigkeit V (Machzahl 
M=V'/co) verlaufen. Man erhalt eine von der Mach- 
zahl abhangige Abweichung vom Cosinusgesetz: 

Ee 


(1—M sina) cosa. (11) 


p2 
oc 

Im vorliegenden Fall der gekriimmten Phasen- 
flachen und der inhomogenen Stromung wird tber 
die Phasenflache integriert, wobei M und a Funktio- 
nen der Kanalhéhe y werden. Der Ubergang zu vor- 
eilender Phase in Absorbernahe bei groBen Wind- 
geschwindigkeiten andert das Vorzeichen von cosa 
und damit von £ , doch kann bei der Auswertung 
stets |Z | genommen werden, da in der Herleitung 
von Gl. (11) die Wechselschnelle des Schallfeldes 
mit einem willkurlichen Phasenwinkel in w ¢ enthal- 
ten ist. Das Integral erhalt dann die Form: 


p2 ; | 
fe [1—M(y) sin a(y) ] cos a(y) dy |: 
(12) 
Der Faktor 1/B dient zur Normierung. Bei der nu- 
merischen Auswertung wurde das Integral als 
Summe geschrieben und die Werte fir M(y) und 


a(y) aus den Mefkurven entnommen. Man erhalt 
fur EL, die in Tabelle III angegebenen Relativwerte. 


Tabelle III. 
Vente ae leteon eS ig. ab} 55 oF | n/a 
Bie he 1,005 e270) 0,49. | 2,03 | 


Hieraus laBt sich der nach Bild 10 zu erwartende 
Dampfungsanstieg mit wachsender Windgeschwin- 
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digkeit ablesen. Bei der Berechnung der wirklichen 
Kanaldampfung ist nun noch die Anderung des Ab- 
sorptionsgrades der Wand zu beriicksichtigen. Aus 
den Messungen des Strémungswiderstandes errech- 
net man dafiir die Werte der Tabelle IV. 


Tabelle IV. 
Vat Os | dd ssl? COB arabe a Ise mas 
Avm/4a | 1° | 0,95 | 0,86 | 0,79 | 0,63 


Nach der Definition des Absorptionsgrades erhalt 
man nun fiir die absorbierte Energie Ey; : 


Paes = Ey : Aym . 


Die iiber die absorbierte Energie berechnete Kanal- 
dampfung zeigt bei den verschiedenen Windgeschwin- 
digkeiten gute Ubereinstimmung mit den Direkt- 
Messungen der Ausbreitungsdampfung. Die Abwei- 
chungen liegen unter 5%, sowohl bei Vorwarts- wie 
auch bei Riickwartsausbreitung. 

AbschlieBend Ja8t sich sagen, da bei Vorwarts- 
ausbreitung die Anderung der Kanaldimpfung im 
wesentlichen von der Kennwertanderung bestimmt 
wird. Bei Schallausbreitung entgegen der Stromung 
tritt als gegenlaufiger Effekt gleicher GroSenordnung 
eine Verformung der Phasenflachen auf, die in ihrer 
Wirkung den Abfall der Dampfung tberwiegt, der 
von der Kennwertanderung hervorgerufen wird. 


4. Untersuchung der Schallverstarkung 
bei Helmholtzresonatoren 


4.1. Kurze Zusammenfassung der Schallverstarkung 
Der nun folgende Teil der Arbeit behandelt die 


Schallverstirkung in einem Windkanal, der mit 
Helmholtzresonatoren abgedeckt ist. Wegen der 
hohen Dampfung im Resonanzbereich eignen sich 
Helmholtzresonatoren besonders Dampfung 
ausgepragter Gerduschspitzen. Die Resonatoren sind 
in Dreiergruppen in die Kanaldecke eingelassen und 
auf 800 Hz abgestimmt. Der Abstand zweier Grup- 
pen in Stromungsrichtung betragt etwa 3,5 cm. Die 
Untersuchungen von F. Mecuet [1] haben ergeben, 
da8 ein solcher Windkanal neben den Dampfungs- 
bereichen ein Frequenzgebiet aufweist, in dem das 
eingespeiste Schallsignal verstarkt wird. Die Ver- 
starkung erreicht maximal 10 dB/m. Die Frequenz 17 
maximaler Verstirkung andert sich mit der Wind- 
geschwindigkeit nach folgender Dispersionsformel : 


V= 9 [-1+V(4L/con) v. 


zur 


(13) 


V ist die Stromungsgeschwindigkeit, L der Abstand 
zweier Resonatorgruppen. Die Verstarkungsbereiche 
wiederholen sich bei fester Windgeschwindigkeit mit 
der Laufzahl n. 
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In der oben angegebenen Arbeit ist die Schall- 
verstarkung in Analogie zur Verstirkung elektro- 
magnetischer Wellen behandelt und ein analytischer 
Nachweis der Verstarkungsmoglichkeit gegeben wor- 
den. 

Die vorliegenden Untersuchungen erklaren die 
Schallverstarkung aus der Wechselwirkung des 
Schallfeldes mit den Turbulenz- und Strémungsvor- 
gangen im Kanal. 


4.2. Allgemeine Untersuchungen an Helmholtz- 
resonatoren 


Voraussetzung fiir-die folgenden Untersuchungen 
ist eine genaue Kenntnis der Stromungsvorgange im 
Innern des Kanals und in den Halsen der Resonato- 
ren. Zur Messung wurde das weiter oben beschrie- 
bene Hitzdraht-Anemometer benutzt, mit dem Tur- 
bulenz- und Schallschnelle ohne Storung der Feld- 
verteilung aufgenommen werden kénnen. Im wesent- 
lichen wurden zwei Groen des Strémungsfeldes ge- 
messen: 

i. mittlere Geschwindigkeit lV’, 

ii. mittere Turbulenzschnelle Vv? . 


Die Kernstro6mung im Kanal besitzt bei einer 
Stromungsgeschwindigkeit von 45 m/s einen Tur- 
bulenzgrad von 4% in dem oben _beschriebenen 
Kanal mit Helmholtzresonatoren. 

An der Resonator-Berandung des Kanals steigt 
die Turbulenz stark an und erreicht in unmittelbarer 
Nachbarschaft der Resonatoréffnungen 7 bis 10%. 
An der glatten Wand zwischen zwei Resonatoren 
klingt die Turbulenz etwas ab, so dafi dicht an der 
Resonatorwand eine raumlich periodische Turbulenz- 
struktur vorliegt. Diese Schicht ist im Mittel nur 
7mm tief, wie der genaue Verlauf der Turbulenz- 
verteilung in Bild 18 zeigt. An der scharfen Kante 


der Resonatoréffnung wird starke Turbulenz er- 


Resonatorwand 


2 3 ms 


Vee 


Bild 18. Turbulenzyerteilung unter einem Resonator (a) 
beziehungsweise zwischen zwei Resonatoren (b). 


zeugt. Fiir die folgenden Betrachtungen ist es von 
Bedeutung, da hier Energie der Gleichstromung in 
Turbulenzenergie tiberfiihrt wird. 

Die Frequenzanalyse der Turbulenz zeigt im Zen- 
trum des Kanals keine Bevorzugung der Eigenfre- 


- quenz der Helmholtz-Resonatoren. Mit konstanter 


Bandbreite analysiert findet man einen Abfall der 
Turbulenzschnelle von etwa 5 dB/Oktave (Bild 19). 


———_—_____—— a 
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Dicht unterhalb des Resonatorhalses betragt die 
Uberhéhung der Eigenfrequenz im Mittel 1,6 dB. 
Diese schwache Beeinflussung des Turbulenzspek- 
trums durch die Resonatoren beruht darauf, da die 
Turbulenzschwankungen durch Zahigkeitskrafte und 
nicht durch Druckkrafte hervorgerufen werden. Diese 
fast druckfreien Vorgange kénnen das Federvolu- 
men des Resonators nicht zu Schwingungen anregen. 

Entsprechend der verschiedenen Struktur der 
Resonatorseite und der glatten gegentberliegenden 
Kanalseite ist das Stromungsprofil im Kanal unsym- 
metrisch. AuSerdem zeigt die Kernstromung eine 
schwache raumliche Periodizitat, die im Verstar- 
kungsbereich stark hervortritt (Bild 25). Entgegen 
dem oben Gesagten steht die Beobachtung, da8 der 
Kanal bei bestimmten Windgeschwindigkeiten pfeift. 
In diesem Sonderfall erfolgt eine Anregung der 
Resonatoren durch am Resonatorhals periodisch auf- 


' tretende Wirbel, die ttber Druckschwankungen den 
_ Resonator anregen, der seinerseits die Wirbel stabi- 


lisiert. Bei diesem Vorgang wird Schall erzeugt und 
abgestrahlt. Die Frequenzanalyse der Turbulenz 
zeigt nun eine Uberhéhung der Pfeiffrequenz um 
14.dB. Wird zusatzlich in den Kanal ein Schall- 
signal einer anderen Frequenz eingespeist, so ver- 
schwindet mit wachsender Signalamplitude das Pfei- 
fen vollstandig. Im Turbulenzspektrum tritt jetzt 
die Frequenz des eingespeisten Schallsignals stark 
tberhoht auf. 

Beide Analysen sind in Bild 20 dargestellt. Hier 
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Bild 19. Spektrale Verteilung der mittleren Turbulenz. 
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Bild 20. Turbulenzspektrum im pfeifenden Windkanal 
mit (©) und ohne (+) Schalliiberlagerung. 


wird durch den eingespeisten Schall die Selbsterre- 
gung des Resonators gestort und die Turbulenz- 
erzeugung an der Kante des Resonatorhalses durch 
die Schwingung des Resonators im Takt der Signal- 
frequenz gesteuert. Auf diese Weise laBt sich das 
Frequenzspektrum der Turbulenz am Resonatorhals 
durch ein in den Kanal eingespeistes Schallsignal 
stark beeinflussen. 
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An dieser Stelle ist es notwendig, auf den Unter- 
schied zwischen Schall und der oben gezeigten syn- 
chronisierten Turbulenz genauer einzugehen. Im 
Schallfeld sind Druck p und Schnelle tber die cha- 
rakteristische Impedanz Z miteinander verkniipft: 
p/v =Z. Die Schallwellen breiten sich mit der Schall- 
geschwindigkeit c aus, wobei die Gleichung c? = x RT 
gilt. Darin bedeuten: R allgemeine Gaskonstante, 
T absolute Temperatur, ~ Adiabaten-Exponent. 

Bei der synchronisierten Turbulenz, die ihre Ur- 
sache in Schub- und Zahigkeitskraften hat, ist das 
Verhaltnis p/v um Gro8enordnungen kleiner. Ver- 
glichen mit einem Schallfeld ist der Druckanteil sehr 
gering. So besitzt die Turbulenz auch keine eigene 
Ausbreitungsgeschwindigkeit, sie wandert mit der 
mittleren Stromungsgeschwindigkeit V. Wegen der 
nahezu fehlenden Druckkrafte strahlt die Turbulenz 
auch nicht ab. Monochromatische Turbulenz, die sich 
in der Schnellekomponente nicht von Schall unter- 
scheidet, wird vielfach Pseudoschall genannt. Das 
unterschiedliche Verhalten von Schall und Pseudo- 
schall beziiglich der Ausbreitung und der Druckkom- 
ponente macht es moglich, sie meftechnisch zu unter- 
scheiden. . 

Das Druck-Schnelle-Verhaltnis der Turbulenz in 
Vergleich zu den SchallfeldgroBen wurde genauer 
untersucht. Fiir die Druckmessung wurde die schon 
besprochene Sonde benutzt, deren Formgebung eine 
zuverlassige Messung der Druckschwankungen ge- 
wahrleistet. Die Schnellemessung erfolgte mit dem 
Hitzdraht und anschlieBendem Filter von 10 Hz 
Bandbreite. Das Schallsignal wurde auf die gleiche 
Schnelle wie die gefilterte Turbulenz eingestellt. Die 
Messung erfolgte bei verschiedenen Windgeschwin- 
digkeiten. Unter der Voraussetzung gleicher Schnelle 
wurden fiir das Verhaltnis von Schalldruck zu Tur- 
bulenzdruck P;/Py im schallharten Kanal folgende 


Werte gemessen: 


40 | 60 | m/s 
| 99 | dB 


Strémungsgeschwindigkeit | 20 
20 lg Ps/Pp |. 37 -| 29 

Mit wachsender Stromungsgeschwindigkeit nimmt 
der Druckanteil der Turbulenz zu, da die Stromung 
immer mehr als kompressibel betrachtet werden 
mu. 

Unter diesen Gesichtspunkten ist es auch zu ver- 
stehen, dai nur bei giinstigen Bedingungen eine An- 
regung der Helmholtzresonatoren durch die Turbu- 
lenz stattfindet. Im folgenden wird der Kanal bei 
45 m/s Windgeschwindigkeit betrieben, ein Bereich, 
in dem keine Selbsterregung auftritt. 

Bild 21 zeigt das Ergebnis einer genauen Unter- 
suchung des Synchronisationseffektes. Analysiert 
wurde die Turbulenz mit und ohne Uberlagerung 
des Schallsignals des Verstarkungsmaximums. Die 
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Bild 21. Differenz der Terzbandenergien des Rausch- 
spektrums mit und ohne Schalliiberlagerung 
(Relativwerte). 


Anzeige des Terzfilters ist auf den Energieinhalt der 
einzelnen Terzbereiche umgerechnet und die Diffe- 
renz der Energien mit und ohne Schalliberlagerung 
als Ordinate aufgetragen. Aus der Kurve wird deut- 
lich, daB® bei der Synchronisation ein Teil der Ener- 
gie, besonders der tieferen Frequenzbander, in den 
Bereich der Synchronisationsfrequenz tberfihrt 
wird. Die Schallschnelle tiberlagert sich nicht einfach 
dem weiB®en Turbulenzrauschen, sondern die Stro- 
mungsenergie, die an der Resonatorkante in Turbu- 
lenz umgewandelt wird, wird spektral gefarbt. Durch 
den schwingenden Resonator wird an jeder Resona- 
torkante Pseudoschall erzeugt. 

Zur Messung der Grofe des Synchronisations- 
effektes wurde am ersten Resonator der Verstar- 
kungsleitung die Schallschnelle ohne Windiberlage- 
rung dicht unterhalb des Resonatorhalses bei kon- 
stantem Pegel und verschiedenen Frequenzen ge- 


messen (Kreise) (Bild 22). Durch die 
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Bild 22. Schallschnelle an der Miindung des Resonator- 
halses. Bei 0,8 und 1,3 kHz (x): Schallschnelle 


bei 45 m/s Windiiberlagerung. 
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Dampfung der Resonanzfrequenz ist nur eine 
schwache Uberhéhung der Resonanzschnelle zu fin- 
den. Bei den Frequenzen 800 und 1300 Hz sind die 
im Resonatorhals gemessenen Schnellen bei 45 m/s 
Windiiberlagerung ebenfalls eingetragen (Kreuze). 
Schall und Pseudoschall sind tiberlagert und in der 
Schnellekomponente nicht voneinander zu trennen. 
Die dritte, gestrichelte Kurve in dem Diagramm 
zeigt den Frequenzgang der Turbulenzschnelle ohne 
Schalliiberlagerung. Die Uberhéhung der Schall- 
schnelle durch den Pseudoschall ist wesentlich gr6- 
Ber als die ohne Schalliiberlagerung in diesem Fre- 
quenzbereich vorhandene Turbulenzschnelle. Dies 
la8t sich nur so erklaren, daf} Energie aus den be- 
nachbarten Frequenzbereichen synchronisiert wird. 
Die gemessenen Geschwindigkeitsprofile der Ka- 
nalstro6mung geben einen weiteren Hinweis dafiir, 
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daf} die Energie der synchronisierten Turbulenz der 
Str6mung entnommen wird. Das tiberlagerte Schall- 
feld andert besonders in der Nahe der Resonator- 
offnungen das Stromungsprofil so, daB vor der 
Resonatorwand ein starker Geschwindigkeitsgradient 
und damit groBe Schubspannungen auftreten. In die- 
sem Bereich findet eine erhohte Turbulenzerzeugung 
statt, die fiir den Synchronisationseffekt notwendig 
ist. : 

Da der Pseudoschall sich mit der mittleren Stro- 
mungsgeschwindigkeit bewegt, mufs der synchroni- 
sierte Turbulenzbereich den Isotachen des Kanals 
folgen. Er erfiillt nicht, wie das Schallfeld, den gan- 
zen Kanal. Frequenzanalysen der Turbulenz zeigen, 
daB die Spektrallinie des Pseudoschalls mit wach- 
sendem Abstand von der mit Resonatoren besetzten 
Wand des Kanals sehr rasch abklingt. Die Pseudo- 
schall-Zone ist im Mittel nur 7 bis 10 mm stark. Auf 
die ,,Lebensdauer“ des Pseudoschalls, der mit der 
ungesteuerten Turbulenz in starke Wechselwirkung 
tritt, wird bei der Diskussion der Verstéarkung noch 
genauer eingegangen. 

Die Richtungsverteilung der Turbulenzschnelle 
vor und hinter dem Resonatorhals ist in Bild 23 a 
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Bild 23. Richtungsverteilung der Turbulenz unterhalb 
der Resonatoréfinung (a); Anteil des Pseudo- 
schalles (b). 


dargestellt. Zu dieser Messung wurde die Richtungs- 
selektivitat des Hitzdrahtes ausgenutzt, der die Ge- 
schwindigkeitsanderungen normal zur Drahtachse 
bevorzugt registriert; seine Richtcharakteristik ent- 
spricht der eines Dipols. Der Vergleich der Messun- 
gen mit und ohne Schalliiberlagerung zeigt ein Auf- 
treten von Pseudoschall nur an der Hinterkante des 
Resonatorhalses. Dieser Pseudoschallanteil ist in 
Bild 23 b gesondert aufgetragen. Die Polarisation 
des Pseudoschalles weicht um 40° stromabwarts von 
der Wandnormalen ab. Wahrend bei ruhender Luft 
die Polarisation der Resonatorschnelle am Hals nor- 
mal zur Wand liegt, wird bei Strémungsiiberlage- 
rung die austretende Wechselschnelle in Strémungs- 
richtung ,,verweht*. 

Die oben beschriebene Erzeugung von Pseudo- 
schall ist die Ursache fiir die im Kanal auftretende 
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Verstarkung. Einen Hinweis hierfiir geben schon die 
friiheren Untersuchungen an schwach bedampften 
Helmtholtz-Resonatoren [1]. Es wurde die Ausbrei- 
tungsdampfung im Kanal bei zwei verschiedenen 
Typen untersucht. Dabei war die dampfende Gaze 

i. auf der Kanalseite der Resonatorhalse, 

ii. auf der Federvolumenseite der Resonatorhalse 

angebracht. 
Im ersten Fall erreicht man eine glatte Berandung 


_ des Kanals, die Turbulenzbildung an der Kante der 


Resonatoren wird weitgehend unterdriickt. Dement- 
sprechend zeigt auch die Frequenzkurve der Damp- 
fung keine Besonderheiten. Im zweiten Fall wird 
der Mechanismus der Turbulenzerzeugung nicht ge- 
stort. Die Schallschnelle im Hals wird jedoch ge- 
dampft und damit die Synchronisation der Turbu- 
lenz stark herabgesetzt. Die Dampfungskurve zeigt 
schwach ausgepragte Entdampfungsbereiche, die mit 


_ den Verstaérkungsgebieten der unbedaémpften Reso- 


natoren zusammenfallen. 


4.3. Der Verstarkungsmechanismus 


Mit dem absolut-geeichten Hitzdrahtanemometer 
wurde nun die Turbulenzverteilung normal zur Re- 
sonatorseite des Kanals an verschiedenen Stellen des 
Kanals gemessen, um die Anderung der Turbulenz- 
struktur bei Uberlagerung des Schallsignals genau 
zu untersuchen. Uberlagert wurde die Frequenz des 
Verstarkungsmaximums bei 45 m/s Windgeschwin- 


_digkeit. Aus Griinden der Absoluteichbarkeit sind 


die Messungen breitbandig ausgefiihrt, doch zeigt 

die Frequenzanalyse, daB® die Anderung der Turbu- 

lenzschnelle bei Schalliiberlagerung durch die stark 
iiberhéhte Spektrallinie des Pseudoschalles hervor- 
gerufen wird. 

Diese Turbulenzdifferenz A 2? ist in Bild 24 maB- 
stablich in einen Langsschnitt des Kanals eingetra- 
gen. Daf es sich hierbei um Pseudoschall und nicht 
um Schallschnelle handelt, wird durch folgende 
Punkte belegt: ; 

i. Vor dem ersten Resonator ist A/v?=0; ein 
Schallsignal wird jedoch in den Kanal ein- 
gespeist. 

ii. Der Pseudoschall ist drtlich begrenzt. 

iii. Er wandert mit den Isotachen (Bild 25) der mitt- 
leren Stromungsgeschwindigkeit. 

iv. Infolge der starken Wechselwirkung mit der un- 
gesteuerten Turbulenz verschwindet der Pseudo- 
schall nach etwa 5 cm Laufweg. 

v. Er tritt nur in der Umgebung der Verstarkungs- 
frequenz auf. 

Die Bedeutung von Punkt v. wird bei der Diskus- 
sion der Phasenbeziehungen zwischen Schall und 
Pseudoschall deutlich. 

Die in Bild 24 dargestellten MeBergebnisse sol- 
len nun kurz erlautert werden: Die einzelnen Reso- 
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natoren werden mit k=1,2,... 

Schallausbreitung numeriert. 

k= 1. Vor dem 1. Resonator ist AVo?=0, Pseu- 
doschall ist nicht vorhanden. Im Resonator- 
hals wird Turbulenz synchronisiert, der Be- 
reich AVv?>0 ist nach unten scharf be- 
grenzt. Hinter dem Resonator bewegt sich 
dieser Pseudoschallbereich entlang 
Isotachen mit Windgeschwindigkeit. 

. Die Ausbildung der Pseudoschallbereiche 
hat sich stabilisiert. Vor dem Hals ist der 
Pseudoschall des Resonators k=9 zu fin- 
den, im Resonatorhals wird neuer Pseudo- 
schall erzeugt, hinter dem Resonator sind 
beide Bereiche schwach getrennt zu finden. 
Im weiteren Verlauf lost sich der 4ltere 
Bereich auf, bei 

k=11 wiederholt sich der gleiche Vorgang. 

Hierbei ist folgendes wichtig: Die Pseudoschall- 
bereiche lassen sich lokalisieren und stromabwarts 
verfolgen. Der bei k erzeugte Pseudoschall kann mit 
dem bei k + 1 erzeugten in Wechselwirkung treten. 

Bei der Betrachtung der Wechselwirkung und der 
Phasenbeziehungen wird die unterschiedliche Ge- 
schwindigkeit von Schall und Pseudoschall wesent- 
lich. Es gelten folgende Beziehungen: 


in Richtung der 


einer 


eesti dig- Frequenz | Wellenlange 
Schall c v As = c/v 
Pseudoschall V v Ap = Viv 


Die Phasenlage der Helmholtz-Resonatoren wird 
von der anregenden Schallwelle festgelegt. Die Pha- 
sendifferenz zweier Resonatoren ergibt sich aus 


ihrem Abstand L und der Wellenlange 4 zu: 


miles 
Cc 


Oea2nk (14) 


As 


P+ 


Ahnlich ist die Phasenlage des am k-ten Resonator 
erzeugten Pseudoschalles. Die Phasendrehung der 
Pseudoschallwelle auf dem Weg L ist wegen der 
kleineren Ausbreitungsgeschwindigkeit grofer: 

1g Ly 

SO eg gt 

Wk+1— Wk a, V 
Zur Vereinfachung wird ~, = yx =0 gesetzt. Sollten 
nun der am Ort k und k+1 erzeugte Pseudoschall 


am Ort k+1 die gleiche Phase besitzen, so muh 
gelten: 


(15) 


VYk+1—Pki1=20n Tiel aD ateyotess (16) 
daraus folgt: 
1 ae Een (17) 
c V 


Diese Gleichung zeichnet ein Wertpaar (vy, V) aus 
(periodisch in n), bei dem die Uberlagerung des 
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Bild 25. Isotachen der Stromung im Kanal. 


Pseudoschall 


Bild 26. Phasenbeziehung zwischen Schall- und Pseudo- 
schall (die Schallwellenlange ist verkiirzt) . 
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Bild 27. Anderung des Schalldruckes in verschiedenem 
Abstand von der Resonatorwand (+ — —-+ 


ohne Windiberlagerung). 


Pseudoschalles k und k+1 am Ort k+1 eine Ver- 
groBerung der Amplitude ergibt. Fir andere Fre- 
quenzen (bei konstantem ) ist der Mittelwert der 
Uberlagerung kleiner und ergibt bei einer Phasen- 
differenz x Ausloschung des Pseudoschalles. Aus der 
Gl. (17) erhalt man eine Dispersionskurve fiir die 
ausgezeichneten (vy, ’)-Paare. Die Schallgeschwin- 
digkeit c¢ setzt sich additiv aus der Schallgeschwin- 
digkeit bei ruhender Luft und der Windgeschwindig- 
keit V zusammen: c=cy+V . Damit erhalt man aus 
Gl. (17) die Dispersionsformel 


V=D(-1+V1+(4L/eon)y). (18) 
Sie ist identisch mit der Dispersionsformel, die von 
F. Mecuer auf Grund der Wanderfeldroéhren-Analo- 
gie aufgestellt wurde. Die Herleitung dieser Formel 
aus den Phasenbeziehungen zwischen Schall und 
Pseudoschall hat demgegeniiber den Vorteil groRe- 
rer Anschaulichkeit. Diese Beziehung, die in wenig 
abgewandelter Form auch von F. Mecuex benutzt 
wird, ergibt mit den experimentellen Untersuchun- 
gen des Pseudoschallfeldes eine einfache Erklarung 
der Schallverstarkung. 

Die Wertepaare (V,¥) der Maxima der gemesse- 
nen Schallverstarkung liegen auf dieser Dispersions- 
kurve. 

Fiir die nachfolgende Berechnung der Verstar- 
kungsbereiche ist der Zusammenhang zwischen Fre- 


es 


Umgebung der Frequen- 
zen der Verstarkungs- 
maxima vy, notwendig. 
LaBt man in der Glei- 
chung fiir den Synchro- 
nismus Gl. (17) eine 
Phasenabweichung Ag 
zu, so erhalt man fir 
diesen Zusammenhang 
bei konstanter Geschwin- 
digkeit V 


Ap=2x(2 =) 
UP 
(19) 


Die Phasenlage zwischen Schall und Pseudoschall in 
einem festen Zeitpunkt t=t) ist in Bild 26 dar- 
gestellt. Am 1. Resonator ist y,; =~, =0 gesetzt. Der 
Zusammenhang gilt fiir das 1. Verstarkungsmaxi- 
mum, das bei der benutzten Windgeschwindigkeit 
bei vy; = 1,35 kHz liegt. 

Phasengleichheit von Schall und Pseudoschall 
liegt nur in unmittelbarer Umgebung der Resonator- 
offnungen vor. Hier findet auch die Kopplung zwi- 
schen Schall und Pseudoschall statt. 


Die Starke der Kopplung und die daraus resultie- 
rende Verstarkung ist aus den Schalldruckkurven in 
Bild 27 abzulesen. Die Messungen wurden mit der 
bereits mehrfach erwaéhnten Sonde ausgefiihrt. Zum 
Vergleich ist in Bild 27 der Schalldruckverlauf bei 
der gleichen Frequenz ohne Windiiberlagerung ein- 
gezeichnet. 


Verfolgt man den Verlauf der Druckkurven, so 
findet man in allen Kanalhdhen einen Zuwachs des 
mittleren Schalldruckes um 4% pro Resonator bei 
der Frequenz des Verstérkungsmaximums. Integra- 
tion dieser Zuwachsrate tiber einen Meter Resonator- 
strecke ergibt mit 10 dB eine gute Ubereinstimmung 
mit der bei diesem Wertepaar (vy, ’) gemessenen 
Maximalverstarkung gegentiber ruhender Luft. Die 
Kurven in Bild 27 sind Relativanderungen des 
Schalldruckes. 

Die teilweise bzw. vollstandige Ausléschung des 
Pseudoschalles bei einer Phasendifferenz zwischen 
den Bereichen k und k+1 am Ort k+1 ist aus dem 
Frequenzverlauf der Uberlagerung von Turbulenz 
und Pseudoschall in Bild 28 deutlich zu sehen. Hier 
wurde dicht hinter einem Resonatorhals die Wechsel- 
schnelle mit dem Hitzdraht gemessen und mit 10 Hz 
Bandbreite analysiert, wobei die Frequenz des iiber- 
lagerten Schallsignals synchron mit der Mittenfre- 
quenz des Filters verandert wird. Die Kurve weist 
zwei Maxima auf, die den Verstarkungsbereichen vy, 
und vy, entsprechen. Gleichzeitig ist der Frequenz- 
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Bild 28. Frequenzanalyse der Turbulenz im Pseudo- 
schallbereich. Die Signalfrequenz wird syn- 
chron mit der Filterfrequenz geandert. 


gang der Turbulenz ohne Schalliiberlagerung ein- 
getragen. Die bei vy, gemessene Uberhdhung ist ge- 
ringer, weil hier wegen der doppelten Phasendre- 
hung der Pseudoschallwelle gegeniiber dem ersten 
Verstarkungsbereich die Storung der Pseudoschall- 
welle durch die nichtsynchronisierte Turbulenz gré- 
Ber ist. 

Zusammenfassend Jat sich sagen: Pseudoschall 
wird zwar bei jeder Frequenz an jedem Resonator 
erzeugt, wegen der Phasendifferenz der einzelnen 
Wellenzige ist jedoch eine stabile Pseudoschallwelle 
nur in direkter Nachbarschaft der Frequenzen mit 
Phasensynchronismus zwischen Schall- und Pseudo- 
schallwelle moglich. Damit ist auch nur in diesen 
Bereichen eine Verstarkungsmoglichkeit der Schall- 
welle gegeben. Die Kopplung zwischen der ankom- 
menden Pseudoschallwelle und dem Schallfeld im 
Helmholtzresonator kann tiber die Druckkompo- 
nente geschehen, die der an der Resonatordfinung 
vorbeistromende Pseudoschall hervorruft. Diese im 
Verstarkungsfall in der richtigen Phasenlage auf- 
tretende Druckanderung unterstiitzt die Schwingung 
im Helmholtzresonator. Der Verlauf der Schalldruck- 
kurven, die vor der Resonatoréfinung von einem 
relativen Minimum in ein Maximum ibergehen 
(Bild 27) spricht fiir diese Kopplungsart, die den 
von der Impedanz des Helmholtzresonators vor- 
geschriebenen Druckverlauf vor dem Resonator stort. 
Dieser Mechanismus wird von der raumlichen Perio- 
dizitat der Kanalstromung (Isotachen, Bild 25), die 
ebenfalls im Verstarkungsbereich besonders aus- 
gepragt ist, unterstttzt. 

Mit Hilfe des gemessenen Kopplungsfaktors K 
und dem Zusammenhang zwischen einer Frequenz- 
verschiebung und der daraus folgenden Phasendif- 
ferenz zweier Pseudoschallbereiche 


dp =2n( =n}. K=4% 


VY, 


(19) 


lassen sich die GroSfe und der Frequenzgang der 
Verstarkung im Verstarkungsbereich berechnen. Die 
Uberlagerung der normierten Schallamplitude 


A,=sin@ t 
mit dem Kopplungsanteil der Amplitude des Pseudo- 
schalles, der eine Phasendifferenz Aq@ besitzt, 
A, =K sin(w t+ Aq) 
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ergibt: 
A=sinwt+Ksin(wi+Ag), (20) 


A=K [Fe + cos Ae) sinwt+sinAgcosa@t|. 


Eliminierung des zeitabhangigen Gliedes ergibt fiir 
die Maximalamplitude der Uberlagerung der beiden 
Sinusschwingungen: 


A= (14+ K?4+2K cos Aq)". (21) 


Dieser Amplitudenwert gilt fiir den Ubergang zwi- 
schen zwei Resonatoren. Unter Beriicksichtigung des 
Abstandes L zweier Resonatoren erhalt man bei 
28 Resonatorgruppen auf einem Meter Laufweg der 
Welle fiir die Amplitudenzunahme pro Meter: 


A’ = (14+ K?+2 K cos Aq) 14 (22) 


und nach Einsetzen der Gl. (19) fiir den Frequenz- 
gang der Verstarkung: 


D(v) =280 log |1+ K?+2Kcos2x (.". a) db, 
VY, m 

(23) 

Die Gl. (23) wurde fiir n=1 und vy, = 1,35 kHz bei 

45 m/s Windgeschwindigkeit ausgewertet. Die ge- 

rechnete Kurve ist zusammen mit den gemessenen 

Verstarkungswerten in Bild 29 eingetragen. Die 
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Bild 29. Berechneter und gemessener Verstarkungsver- 
lauf in einem Verstarkungsbereich. 


schlechte Ubereinstimmung bei tiefen Frequenzen 
beruht auf der Annaherung an die Eigenresonanz 
der Helmholtz-Resonatoren und der damit verbun- 
denen Dispersion der Signalwelle. 


4.4. Einige Untersuchungen zur Rickwarts- 
ausbreitung. 


Zum Abschlu8 sollen noch kurz einige Untersu- 
chungen zur Schallausbreitung entgegen der Stro- 
mungsrichtung gebracht werden. Die bei Vorwarts- 
ausbreitung auftretenden, ausgepragten Verstar- 
kungsbereiche sind nicht zu finden. Die Frequenz- 
kurve der Dampfung zeigt schwache Entdampfungs- 
gebiete (1 bis 2 dB), deren Mittenfrequenz auf der 
Dispersionskurve liegt, die nach der Wanderfeld- 
rohren-Analogie fiir Riickwartswellen zu erwarten 
ist [1]. Im Gegensatz zur Vorwartsausbreitung ist 
diese Entdampfung sehr stark amplitudenabhingig. 
Betrachtet man diese Entdimpfung nach den oben 
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diskutierten Gesichtspunkten der Verstarkung, so 
sind verschiedene Unterschiede festzustellen: 

i. Ein synchronisierter Turbulenzbereich lauft ge- 
gen das Schallfeld an. Die Bedingung fiir Pha- 
sensynchronismus fihrt wegen des gréferen 
Laufweges (zu der Phasendrehung der Pseudo- 
schallwelle addiert sich die Phasendrehung der 
Schallwelle zwischen 2 Resonatorgruppen) zu 
einer geanderten Dispersionskurve, die wieder 
mit der aus der Wanderfeldrohrenanalogie her- 
geleiteten identisch ist 


V = 2[1-Vl—(4L/egn) J. (24) 

ii. Wegen der entgegengesetzten Ausbreitung von 
Schall und Pseudoschall ist auch die Phasen- 
drehung beider Wellen entgegengesetzt. Am Ort 
des Resonators, wo die Wechselwirkung zwischen 
Schall und Pseudoschall stattfindet, ist dadurch 
die Verweilzeit bei annahernd gleicher Phase 
gegeniiber der Vorwartsausbreitung verkiirzt. 

iii. Die Schallschnelleverteilung vor dem Resonator- 
hals, die durch die Windstromung beeinflubt 
wird, zeigt eine Unsymmetrie, die bei Riickwarts- 
ausbreitung fiir die Synchronisation der Turbu- 
lenz ungiinstig ist. Eine Messung der Schnelle- 
verteilung bei 30m/s Windgeschwindigkeit ist 


in Bild 30 wiedergegeben. An der — in Wind- 


Schallschne!le—= 
— — Gw — NH 


Resonatoroffnung 


Bild 30. Intensitat der Schallschnelle vor der Resonator- 
offnung. 


richtung gesehen — vorderen Kante des Resona- 
torhalses wird die Turbulenz erzeugt, dort ist 
aber die Schallschnelle wesentlich kleiner als an 
der gegentiberliegenden Kante. 


Nach diesen Uberlegungen la8t sich schon eine 
schwachere Wechselwirkung zwischen Schall und 
Pseudoschall erwarten, die unter Umstanden itber- 
haupt das Entstehen einer Pseudoschallwelle, die 
iiber eine langere Strecke stabil ist, nicht zur Aus- 
bildung kommen lat. Frequenzanalysen der Tur- 
bulenz bei starrer Synchronisation von Signalfre- 
quenz und Filterfrequenz entsprechend Bild 28 zei- 
gen demzufolge bei Rickwartsausbreitung nur eine 
sehr schwache Uberhéhung der Entdaimpfungs- 
bereiche, die nach der Dispersionskurve zu erwarten 
sind. Analysen bei verschiedener Windgeschwindig- 
keit sind in Bild 31 eingetragen. 
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Bild 31. Frequenzanalyse entsprechend Bild 28 bei 
Riickwartsausbreitung und verschiedenen Wind- 
geschwindigkeiten. 


Hieraus wird ersichtlich, daB auch bei Ruickwarts- 
ausbreitung eine Wechselwirkung zwischen Schall- 
feld und Stromung vorliegt, die dem Verstaérkungs- 
mechanismus bei Vorwartsausbreitung entspricht. 
Wegen der wesentlich ungiinstigeren Bedingungen 
kommt jedoch nur eine schwache Entdampfung zu- 
stande. 

Ein Vergleich der vorliegenden Untersuchungen 
zur Schallverstarkung mit den Arbeiten von 
F. Mecuer [1] laBt erkennen, dai beide trotz ver- 
schiedener Ausgangspunkte zu der gleichen Disper- 
sionsformel fiir den Zusammenhang zwischen der 
Frequenz maximaler Verstéarkung und der Wind- 
geschwindigkeit fiihren. Dies beruht auf der Existenz 
der Pseudoschallwelle, die sich mit der mittleren 
Stromungsgeschwindigkeit bewegt. Wahrend jedoch 
F. Mecuet zur Erklarung der Schallverstarkung die 
Analogie zur elektrischen Wanderfeldrohre benutzt, 
wird in der vorliegenden Arbeit ein Verstarkungs- 
mechanismus an Hand der gemessenen Stromungs- 
und Schallvorginge entwickelt. AuBerdem wird aus 
den Phasenbeziehungen zwischen Schallfeld und 
Pseudoschall der Frequenzgang der Verstaérkung in 
einem Verstarkungsbereich berechnet. Im Gegensatz 
zu der Wanderfeldrohren-Analogie haben diese 
Untersuchungen ergeben, dafs die Turbulenz der 
Kanalstr6mung wesentlich fiir das Eintreten der 
Schallverstarkung ist. Die Umwandlung von Stro- 
mungsenergie in Schallenergie geschieht tiber die 
Turbulenz als Zwischenstufe, wobei die von dem 
schwingenden Resonator hervorgerufene Verfor- 
mung des Strémungsprofils die Turbulenzerzeugung 
unterstutzt. 

Bei der Schallausbreitung im Windkanal entlang 
eines pordsen Absorbers werden zwei Mechanismen 
untersucht, tiber die die Stromung die Schalldamp- 
fung beeinfluBt. Die quantitative Behandlung gibt 
bei Vorwarts- und Riickwartsausbreitung im Rah- 
men der MeSgenauigkeit die gleiche Dampfungs- 
anderung, die bei der Direkt-Messung der Kanal- 
dampfung gefunden wurde. 

Herrn Prof. Dr. Dr. E. h. E. Meyer mochte ich an 
dieser Stelle fiir die Stellung der vorliegenden Arbeit 
und fiir die fordernde Anteilnahme bei der Durch- 
fiihrung danken. 
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Sommaire 


Aprés un rappel préliminaire des conditions que doivent satisfaire les méthodes de 
mesures acoustiques, on expose les solutions adoptées pour |’étude des champs sonores de 
trés grande intensité, Les dispositifs expérimentaux comportent des chaines cinématiques de 
déplacement des microphones et une télécommande d’analyse du bruit effectuée par com- 
mutation automatique des bandes de fréquences. 


Zusammenfassung 


Nach einer einleitenden Zusammenstellung der Anforderungen, denen die akustischen 
Mefimethoden geniigen miissen, werden die zur Messung von Schallfeldern sehr hoher Inten- 
sitat angewandten Verfahren dargestellt. Der experimentelle Aufbau enthalt Transport- 
bander zur Bewegung der Mikrophone und eine Fernsteuerung des Gerduschanalysators, 
die durch automatisches Umschalten der Frequenzbereiche erhalten wird. 


Summary 

After a preliminary review of the conditions which must be satisfied by the methods of 
acoustic measurement, the resolutions adopted are applied to the study of sound fields of 
high intensity. The experimental apparatus comprises movable chains of displacement 
microphones and a controlled analysis of the noise is effected by automatic commutation of 


the frequency bands. 


1. Introduction 


Les problémes acoustiques des avions et des 
fusées ont vu leur importance s’accroitre au cours 
des derniéres années. Ceci résulte de la généralisa- 
tion de la propulsion par réaction et de la mise en 
ceuvre de réacteurs de plus en plus puissants. La 
nature et l’intensité des sources de bruit s’en sont 
trouvées profondément modifiées. 

En effet, au spectre acoustique discontinu (spectre 
de raies).émis par les groupes motopropulseurs 
s'est substituée une répartition spectrale continue, 
spécifique de tout bruit d’origine aérodynamique. 
En méme temps, énergie acoustique rayonnée a 
augmenté parfois dans des proportions considéra- 
bles, entrainant dans les structures des résonances 
et des phénoménes de fatigue sonore demandant des 
études particuliéres. 

En ce qui concerne les mesures de bruit, la néces- 
sité est apparue d’asservir le microphone de mesure 
par une chaine cinématique fonctionnant sans inter- 


vention extérieure et conformément a un programme 
imposé par les essais. 

Les procédés basés sur ce principe et les solutions 
adoptées sont décrits dans ce travail aprés un rappel 
préliminaire consacré aux méthodes d’analyse et 
d’enregistrement des pressions sonores. 


2. Principes de la mesure et de I’analyse du bruit 
2.1. Définition du niveau 


Le niveau sonore en décibels, en un point d’obser- 
vation donné, est défini dans l’intervalle de fréquen- 
ces f,, fy a partir de la connaissance de la pression 
sonore efficace p dans cet intervalle: 


N (fas fo) =10 logy P*(fa> fo) dB, 
Po 


Po étant la pression de référence de 20 uN/m? . 


L’intervalle f,. f,, peut contenir en particulier la 
quasi-totalité des fréquences audibles, soit 20 a 
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10000 Hz. La valeur obtenue sera alors le niveau 
global Vg de la pression sonore efficace et est mesu- 
rée soit directement au moyen d’un décibelmétre, 
soit calculée par l’expression 


Ng=10 logio >” 10%/40 dB. 


N, étant le niveau sonore efficace de chaque bande 
de fréquences, si cette information est directement 
accessible. 


2.2. Analyse dun bruit 


L’analyse d’un bruit s’effectue par intervalles 
de fréquences dont la largeur doit étre choisie sui- 
vant la nature du bruit. 

Lorsque le phénoméne a étudier' est constitué 
principalement d’une suite de fréquences discrétes, 
on a recours a l’analyseur du type hétérodyne 
(AF=const.), ou plus généralement a celui a 
contre-réaction sélective (AF/F =const., F étant la 
fréquence centrale de la bande AF analysée) . 

En raison de sa grande sélectivité (de lordre de 
80 dB), ce qui conduit 4 une importante constante 
de temps, l’analyseur hétérodyne ne peut convenir 
que pour |’étude d’un spectre a raies stables, c’est-a- 
dire exempte de modulation en amplitude ou en fré- 
quence en fonction du temps, par exemple: un 
moteur électrique tournant a régime constant. 

Lorsque le spectre acoustique est fluctuant ou in- 
stable, ce qui est le cas général dans nos études, 
Yanalyseur a largeur relative constante est le mieux 
adapteé. 

L’étude d’un bruit a spectre continu s’effectue au 
moyen d’un analyseur a bandes de largeur relative 
constante, et plus ou moins large, par exemple 1/3 
d’octave ou plus généralement une octave. 


2.3. Chaines usuelles de mesures 


La forme spectrale du bruit n’est pas la seule 
caractéristique déterminante du choix des appareils 
de mesure: la méthode est différente suivant que les 
propriétés statistiques du champ sonore restent sta- 
tionnaires ou, au contraire, subissent des variations 
en fonction du temps. 

Dans le premier cas, la technique des mesures est 
particuliérement simple. L’appareillage se compose 
d’un microphone, d’un analyseur fin ou par bandes 
de fréquences suivi d’un mesureur de niveaux 
(Fig. 1). 

Les caractéristiques balistiques de l'appareil de 
lecture sont habituellement ajustées afin d’obtenir 


microphone 
0-2 . 
ot Bes 
adaptateur décibelmetre 
analyseur 


Fig. 1. Analyse d’un bruit stationnaire. 
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une valeur moyenne des fluctuations, ce qui est pra- 
tique et généralement suffisant. 

L’appareillage offre évidemment la possibilité de 
mesurer le niveau global du bruit. Lorsqu’un état 
sonore évolue en fonction du temps, ou des coordon- 
nées des points de mesure (cas des explorations 
continues), la mesure directe devient alors in- 
suffisante pour la description compléte du phéno- 
méne étudieé. 

On a alors recours a |’enregistrement sur ruban 
magnétique suivi d’analyses effectuées en labora- 
toire. Ce procédé abrége la durée des essais sur le 
terrain ou bien, a durée égale, augmente le nombre 
d’informations_recueillies. Le principe des analyses 
différées est séduisant par son apparente simplicite. 
Le recours a l’enregistrement magnétique complique 
cependant les dépouillements en raison des correc- 
tions liées a la courbe de réponse de l’enregistreur ; 
d’autre part, la restitution des niveaux sonores 
absolus est plus difficile. 

Les caractéristiques d’une chaine de mesure acous- 
tique sont définies par la sensibilité absolue du 
microphone, la bande dynamique transmissible, la 
courbe de réponse de la chaine de mesure, y com- 
pris l’enregistreur magnétique, et enfin la loi d’inté- 
gration d’un signal électrique complexe par l’appa- 
reil de lecture. 


2.4. Intensité des sources sonores 


La technique des mesures dépend également de 
Vintensité de la source du bruit. Or, les sources 
sonores d’intensité trés élevée telles que les fusées 
ou les turboréacteurs sont maintenant nombreuses 
en Aéronautique. Lorsque les niveaux sonores dé- 
passent 140dB environ, ce qui est le maximum 
mesurable par un décibelmétre ordinaire, les micro- 
phones usuels ne conviennent plus. 

En effet, on sort alors du domaine des caracté- 
ristiques linéaires des microphones courants et, dans 
les cas extrémes, on risque méme de les endom- 
mager. I] faut alors utiliser des microphones spé- 
ciaux susceptibles de transmettre sans distorsion 
appreciable des signaux de 180 dB et plus. De plus, 
par raison de sécurité, les expérimentateurs ne peu- 
vent se trouver sur les lieux mémes de mesures. 

Les principes de base étant rappelés, nous allons 
décrire les appareillages d’enregistrement et d’ana- 
lyse de niveaux sonores réalisés par le groupe de 


recherches acoustiques de ’O.N.E.R.A. 


3. Appareillages d’enregistrement et d’analyse 
des pressions sonores 


La connaissance du champ de pression sonore 
dans une cabine d’avion ou autour d’un jet exige un 
trés grand nombre de mesures. D’autre part, le 
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temps imparti aux essais est toujours limité. Dans 
ces conditions, en plus de la précision des mesures, 
notre appareillage doit étre 4 méme d’en assurer la 
rapidite et la streteé. 

Ces considérations nous ont conduits a l’enregis- 


trement automatique des niveaux sonores tele- 


commandé soit par lopérateur chargé du micro- 
phone, soit par un dispositif a relais électriques 
asservi au microphone. Le schéma de _ principe 
(Fig. 2). montre la composition de nos chaines de 
mesures. 


mi adil 


Fig.2. Schéma d’analyse et enregistrement a deux 
canaux. 
mi: microphone, 
adl: adapteur de ligne, 
pol: lecture de la tension de polarisation, 
os: tension stabilisée, ; 
am: amplificateur ou affaiblisseur —20 a 
+40 dB, 
adf: adapteur de filtre, 
fo: filtre d’octave, 
Am: amplificateur, 
EN: enregistreur de niveau, 
EM: enregistreur magnétique 2 pistes. 


-Les caractéristiques les plus intéressantes de cha- 
que chaine sont les suivantes: 

i. La courbe de réponse est rectiligne a+ 1 dB 
dans l’intervalle 40 a 10 000 Hz. 

ii. La bande dynamique est comprise entre 40 et 
150 dB (rapportés 4 20 uN/m?). Cette bande sera 
étendue jusqu’a 165, puis 200dB par d’autres 
microphones. 

iii. La polarisation du microphone €électrostati- 
que qui conditionne sa sensibilité est contrdlée en 
cours d’essai. 

iv. L’étalonnage électronique de chaque chaine et 
de l’enregistreur magnétique est assuré par un oscil- 
lateur local incorporé dans le chassis principal. 

v. L’analyse par bandes de fréquences est auto- 
matique, la commutation des octaves étant obtenue 
par un moteur électrique. Le temps d’arrét sur cha- 
cune des huit bandes est préréglé au_moyen d’un 
compte-pose électronique. 

Le caractére aléatoire du bruit conduit d’ailleurs 
a augmenter la durée d’enregistrement des fluctua- 
tions basse fréquence pour faciliter la lecture de la 
valeur moyenne lorsque celle-ci est choisie comme 
valeur significative. 
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vi. Dans le cas ou le microphone est tenu a la 
main, le déclenchement des analyses est télécom- 
mandé par l’opérateur. 

vii. Enfin, la pression globale ayant été simulta- 
nément enregistrée sur bande magnétique, |’analyse 
fine est obtenue par lecture, sur le méme enregis- 
treur, d’une boucle qui porte |’inscription du signal. 
L’analyseur utilisé est a contre-réaction sélective, 
donc a largeur relative constante. Les potentiométres 
d’accord sont entrainés ici par un moteur électrique. 


4. Utilisation des chaines de mesures 


Nous allons décrire a présent les dispositifs par- 
ticuliers associés a l’appareillage électro-acoustique 
en vue de |’étude du champ sonore des jets et des 
avions a réaction. 

Le champ sonore d’un jet est essentiellement ani- 
sotrope. La répartition spatiale de énergie acousti- 
que est, en outre, différente suivant que l’on con- 
sidére le champ sonore trés proche de |’écoulement 
ou, au contraire, le champ lointain. 

Ce phénoméne est di a limportance des sources 
sonores volumétriques distribuées sur le corps du 
jet [1]. Lorsque les distances deviennent grandes 
par rapport a la longueur d’onde des fréquences 
les plus basses du spectre d’émission, |’anisotropie 
du champ cesse d’étre évolutive et obéit, en particu- 
lier, a la loi classique de décroissance de |’énergie 
acoustique en fonction de la distance (sous réserve 
dune atténuation supplémentaire due a l’absorption 
par l’atmosphere, des obstacles, etc.) . 

L’étude du champ sonore proche est nécessaire 
pour la localisation des «anneaux» émissifs associés 
aux bandes de fréquences qui forment le spectre du 
bruit. L’étude du champ lointain rend compte de 
Yorientation de vecteurs de rayonnement acoustique 
maximum liés a ces bandes de fréquences. 

Ces deux études font appel a des techniques 
distinctes. La méthode expérimentale dépend, en 
outre, de la durée des essais qui sont beaucoup plus 
courts pour les fusées que pour les turboréacteurs. 
C’est pourquoi ces deux cas sont traités séparément. 


4.1. Turboréacteurs (poussées statiques) 
4.1.1. Champ sonore proche 


La mesure précise du champ sonore des jets de- 
vient particuliérement délicate lorsque les investiga- 
tions englobent également le champ trés proche. En 
effet, A une faible distance de l|’écoulement gazeux, il 
importe de connaitre avec précision les coordonnées 
des points d’analyse en raison des variations rapides 
dans l’espace de la forme spectrale du bruit. L’inten- 
sité de bruit est d’ailleurs susceptible d’y atteindre 
des valeurs considérables. En plus, la chaleur déga- 


154 


gée impose des précautions pour maintenir la sen- 
sibilité et méme la vie du microphone. On congoit 
dés lors la nécessité de la substitution aux opérations 
manuelles d’un dispositif entiérement automatique. 
En particulier, la translation du microphone est 


obtenue par un chariot automoteur. Celui-ci, guidé 
par un rail, se déplace suivant un rayon vecteur 
ayant pour origine l’extrémité de la tuyére. Aprés 
chaque série d’essais, une orientation nouvelle est 
donnée a ce rayon pour balayer le secteur choisi, 
repéré en coordonnées polaires (de 10° en 10°) 


(Fig. 3). 


Fig.3. Rail et chariot support de microphone (cliché 
SNECMA, Villaroche, Etude ATAR 9). 


Afin d’abréger la durée des essais, la longueur 
totale des rayons qui est de 32m a été partagée en 
deux parties, de longueur 1m a 16m et 16m a 
32 m. Les enregistrements suivant ces deux trongons 
se font simultanément avec nos deux chaines de me- 
sures. L’enregistrement automatique a été adopté 
uniquement pour le parcours 1m a 16m, le 2eme 
microphone étant déplacé manuellement de 16m a 
32m. 


Translation automatique du microphone 


Les détails du chariot sont visibles sur la Fig. 4. 
On y distingue, a gauche, les boitiers des préampli- 
ficateurs cathodynes; a droite, le moteur d’entraine- 
ment du chariot; au dessus, le boitier des relais de 
commande et les logements du microphone et du 
relais thermique de sécurité. 

Ce bilame intervient, en inversant la marche du 
chariot, dans le cas ot la température ambiante 
atteint 50° 1. Cette précaution, qui sauvegarde les 
performances du microphone de mesures, s'est 
avérée précieuse, notamment lors des enregistre- 
ments des pressions sonores pour le régime de plein 
gaz, post-combustion maximum des turboréacteurs. 

En exploitation normale, l’inversion du sens de 
marche du chariot est obtenue au moyen de butées 
fixées aux extrémités du rail-guide. 


' Cette limite de température peut étre relevée par 
l’emploi d’autres microphones en cours d’étude. 
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Fig. 4. Détail du chariot. 


En raison des variations considérables de la pres- 
sion sonore sur les premiers métres comptés depuis 
la tuyére, on a prévu deux vitesses de translation du 
chariot: une vitesse lente depuis 1m (début des 
mesures) jusqu’a 4m, et une vitesse sensiblement 
double pour le reste du parcours. Le changement de 
vitesse s’effectue automatiquement lors du passage 
du chariot par le point a distance R=4m de lori- 
gine. 

Le repérage de la position du microphone sur le 
rail est indiqué par le chariot lui-méme, qui envoi 
des «tops» dans l’appareillage d’enregistrement, au 
moment de passage par les points choisis de rayons 
R=1m, 2m, 4m, 12m, 16m comptés depuis 
Pextrémité de la tuyére. 


Déplacement manuel du microphone 


Ce mode opératoire a été adopté pour l’explora- 
tion du champ sonore sur les rayons de 16m a 
32 m. 

Le microphone, associé a son préamplificateur 
cathodyne, est fixé sur un support déplacé par un 
opérateur muni de protége-oreilles. 

Lors des mesures, le microphone se trouve dans 
le plan horizontal passant par l’axe du jet, soit a la 
méme hauteur que le second microphone faisant 
partie du chariot. 


Technique des mesures 


La conduite des mesures, semblable pour chaque 
rayon, nécessite comme nous l’avons vu deux chaines 
d’analyse et Lutilisation simultanée d’un enregi- 
streur magnétique a deux pistes. 

Au début de l’essai, le microphone monté sur le 
chariot se trouve au point 16m. Au méme point, le 
second microphone est tenu par l’opérateur. Les 
deux analyses des niveaux sonores par octave de 


he 
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fréquences se font en méme temps, tandis que 
lenregistreur magnétique est mis en marche et fonc- 
tionnera sans interruption jusqu’a la fin de l’essai. 

Ce procédé a pour avantage de comparer la sen- 
sibilité des deux microphones dans les conditions 
réelles de leur utilisation. 

A la fin de chaque analyse, durant laquelle la 
rotation de commutateur de bandes des analyseurs 
par octave se fait d’une maniére automatique, un 
dispositif a relais déclenche la mise en marche du 
moteur de chariot. 

L’opérateur se déplace vers le point R = 24m. 

Parvenu au point R=12m, le chariot s’immo- 
bilise et, par un jeu de relais, déclenche la rotation 
du commutateur de bandes des analyseurs, |’opéra- 
teur étant alors au point R= 24 m. 

Le méme processus se répéte aux points R=8m 
et 32 m. 

Aprés les enregistrements, le chariot se déplace 
vers le point suivant, et l’opérateur revient lente- 
ment vers le point R=16m, en envoyant un top 
dans l’enregistreur magnétique, au passage par le 
point R = 24 m, puis sort du champ. 

A la fin de Panalyse a R= 1 m, si le bilame «limi- 
therme» a permis de venir si prés, le chariot revient 
automatiquement en arriére grace a la butée d’in- 
version du sens de marche, sinon cette inversion est 
déclenchée par le limitherme lui-méme. I] poursuit 
sa course d’une maniére continue et envoie dans les 
enregistreurs des «tops» au moment de passage par 
les points précédemment analysés. La Fig. 5 montre, 
a titre d’exemple, les enregistrements obtenus par ce 
procédé. 

‘Au retour, au point R=16m, le chariot reste 
immobile. 


Fig. 5.Avion RA 29 
Variation des niveaux so- 
nores de 1 4 16m. Azumit 
160°, rapporté a l’axe de 
l’écoulement. 


Enregistrement continu de 
niveau sonore de 1 a 16m 
de la bouche d’un RA 29 
(angle 160° de V’axe du 
jet), analyse suivant 8 oc- 
taves. 
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L’opérateur replace alors son microphone au 
méme point et l’on répéte l’analyse des niveaux sono- 
res. Cette précaution a pour objet de vérifier la sen- 
sibilité du microphone et d’éviter ainsi tout risque 
dincertitude. 

Cette derniére opération clot |’exploration des 
niveaux sonores du rayon étudié. 


Précision des mesures 


Les coordonnées des points d’analyses sur l’aire 
d’expériences, de méme que leurs azimuts, sont 
déterminées avec précision suivant le programme 
dessais. Les erreurs dues uniquement aux positions 
du microphone sont négligeables dans notre procédé 
expérimental. 

Aprés correction des lectures brutes pour tenir 
compte de la sensibilité des microphones, de la 
courbe de réponse des chaines de mesure et du mode 
dintégration du signal électrique par |’enregistreur 
de niveau, l’erreur maximum obtenue est de 


= 155-dB. 


4.1.2. Champ sonore lointain 


Le principe des enregistrements continus est main- 
tenu. Par contre, les lignes des mesures sont situées 
a présent sur des demi-cercles concentriques en rai- 
son de la symétrie axiale des jets. 

L’opérateur, muni du microphone, suit le tracé 
jalonné de 10° en 10° et est accompagné dans sa 
progression par le camion-laboratoire contenant les 
appareils d’enregistrement. Le microphone est pré- 
senté pour commencer a |’azimut coincidant appro- 
ximativement avec la direction de rayonnement 
acoustique maximum. A ce point, on _ effectue 
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Fig.6. Avion de transport a réaction. 
Répartition des niveaux sonores 
suivant un demi-cercle de R=60 m. 
Régime de décollage. 


Enregistrement continue de niveau 
sonore sur un demi-cercle de rayon 
60 m, analyse suivant 2 octaves. 
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Panalyse complete des niveaux sonores par octave 
de fréquences et l’on situe au mieux le signal dans la 
bande dynamique de l’enregistreur magnétique. La 
premiére opération a principalement pour objet de 
verifier, lors des restitutions, la courbe de réponse 
de l’enregistreur magnétique. Les mesures débutent 
a l’azimut de 180° (axe de l’écoulement) lorsque la 
turbulence du jet le permet. Aux points de repéres 
matérialisés par des jalons, lopérateur envoie un 
«top» dans lenregistreur magnétique. Un exemple 
d’enregistrements obtenus est donné par la Fig. 6. 


4.2. Fusées 


L’étude acoustique des fusées exige des techniques 
particuliéres. La durée de fonctionnement est tres 
différente: dans le cas des fusées, elle est comprise 
entre quelques secondes et plusieurs 
D’autre part, contrairement au cas des turboréac- 
teurs, le nombre de mises a feu des fusées est égale- 
ment limité. Ces deux conditions restrictives impo- 
sent des méthodes de mesure distinctes de celles que 
Yon vient d’exposer a propos des turboréacteurs. 


minutes. 


L’abandon des explorations continues, du moins 
pour des poussées statiques, conduit a recourir a un 
grand nombre de chaines de mesures. La technique 
usuelle consiste a utiliser un enregistreur magnéti- 
que a 14 pistes ou davantage et un nombre égal de 
microphones de mesures. Le nombre de micropho- 
nes peut d’ailleurs étre supérieur a celui des voies 
d’enregistrement si la fusée est susceptible de fonc- 
tionner pendant plus d’une dizaine des secondes. I] 
est alors possible d’associer 4 chacune des voies de 
lenregistreur magnétique deux microphones (ou 
plus) que l’on insére successivement par commuta- 
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tion. Le cadrage du signal microphonique, a |’inté- 
rieur de la bande dynamique de lenregistreur 
magnétique, nécessite un systéme spécial d’asservis- 
sement a réponse rapide. 


5. Etude du bruit des avions a réaction 


Les problémes acoustiques posés par les avions 
de transport a réaction peuvent étre groupés en trois 
catégories: 

i. Etude du traitement acoustique de la cabine des 
passagers, 

ii. Etude de la répartition de la pression sonore 
sur la structure des avions, 

iii. Etude des niveaux sonores produits au sol 
lors de décollages, passages ou atterrissages. 


5.1. Etude de la pression sonore dans la cabine 


La théorie de base et les méthodes expérimentales 
mises au point en vue de sélection des matériaux 
acoustiques pour les cabines d’avions ont déja eté 
publiées par l’un de nous [2], [3]. On va décrire 
ici les techniques de mesures que nous avons adop- 
tées pour |’étude des niveaux sonores dans les avions 
en vol, par exemple dans le SE 210 «Caravelle». 

La disposition des appareils dans la cabine de cet 
avion est visible sur la Fig. 7. On y distingue a 
droite les deux racks d’analyse des niveaux sonores 
et a gauche l’enregistreur magnétique. Plus au fond 
on apercoit deux opérateurs munis de microphones. 
Les deux mesures que l’on effectue simultanément ré- 
pondent a des buts distincts. L’une a pour objet de 
déterminer |’état acoustique dans la cabine mesuré 
a la hauteur de la téte d’un passager assis, confor- 
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Fig. 7. Mesure de bruit dans la cabine «Caravelle». 


mément aux spécifications des clauses techniques. La 
seconde opération consiste a déplacer le microphone 
d'une manieére continue prés de l’une des surfaces de 


la cabine. Cette exploration a pour but de déceler 


tout défaut de continuité du traitement acoustique 


-" 


de la surface étudieée. 

On remarquera un dispositif a balais fixé sur 
Yembase du microphone, destiné 4 maintenir cons- 
tant lintervalle entre le microphone et la structure. 

La technique des mesures est différente dans ces 
deux cas. Le premier opérateur commence par en- 


Fig. 8. a) Effets de la discontinuité du traitement acoustique. 
b) Variation des niveaux sonores le long du plafond. 
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registrer en phonie la position du point de mesure, 
puis, agissant sur un contacteur incorporé dans le 
préamplificateur cathodyne du microphone, met en 
marche l’appareillage d’analyse des niveaux sonores 
par octave de fréquences qui s’effectue ainsi en 
direct. Le phénoméne est enregistré en outre sur 
Pune des pistes de l’enregistreur magnétique, en vue 
éventuellement d’analyses fines du bruit. 


Pour la seconde opération, l'une des octaves de 
lanalyseur est insérée dans le circuit de la deuxiéme 
chaine de mesures. D’autre part le courant micro- 
phonique est envoyé sur la seconde piste de |’en- 
registreur magnétique. Le manipulateur part de 
Yarriére de la cabine, puis au cours de sa marche 
envoie périodiquement des «tops» en agissant sur 
Pinterrupteur du courant microphonique incorporé 
dans le boitier du préamplificateur. Ces tops sont 
enregistrés simultanément sur le bathygraphe et 
Yenregistreur magnétique en vue de localiser les 
différents endroits de la surface explorée. Dans ce 
cas, l’analyse complete s’obtient a partir de la resti- 
tution du signal enregistré sur le ruban magnétique. 


Cette méthode d’exploration conduit a des résul- 
tats trés intéressants, comme on peut s’en rendre 
compte a l’examen de la Fig. 8a. Dans cette expé- 
rience, le microphone a été déplacé horizontalement 
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le long de la cloison de la cabine et des parois du 
fuselage. Le «top» indique l’endroit ot la cloison 
rejoint le fuselage. L’important accroissement des 
niveaux sonores observés a ce point est da a l’inter- 
ruption du traitement acoustique, défaut auquel il 
faut évidemment remédier. La Fig. 8b indique la 
répartition des niveaux sonores prés du plafond de 
la cabine. 


5.2. Relevés de la pression sonore sur la structure 
d'un avion a réaction 


La structure d’un avion a réaction est soumise, en 
aval des jets, a un champ sonore dont l’intensité peut 
atteindre des valeurs trés élevées. Les contraintes 
dynamiques qui s’y développent sont alors suscep- 
tibles d’engendrer des phénoménes de fatigue sonore 
de la structure. L’étude systématique des effets de la 
fatigue sonore a conduit SUD-AVIATION a réaliser 
une installation comprenant toute la partie arriére 
de la «Caravelle» (Fig. 9). Cette recherche a pour 
but d’établir la relation entre les valeurs des con- 
traintes, mesurées par des jauges extensométriques 
et la pression sonore excitatrice. 


La répartition de la pression acoustique sur la 
structure est étudiée soit par la méthode dite «point 
par point», soit en faisant appel a l’exploration con- 
tinue. 


5.2.1. Méthode «point par point» 


Cette méthode convient dans le cas d’une surface 
relativement petite telle que la «chapelle» de l’avion. 
Le dispositif expérimental comprend un certain 
nombre de microphones. Chaque microphone asso- 
cié a son préamplificateur cathodyne est fixé sur la 
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Fig. 9. Structure arriére 
du fuselage «Caravelle». 


Fig. 10. Mesure de pression sonore sur la structure 
arriére «Caravelle» — Chapelle. 


structure au moyen de ventouses de miroitier 
(Fig. 10). L’adhérence ainsi réalisée et la souplesse 
d’emploi se sont révélées a |’usage trés satisfaisantes. 
La liaison entre le boitier porte-microphone et les 
ventouses est assurée au moyen d’un pivot réglable. 
Les cables microphoniques sont reliés 4 une boite 
de commutation ayant deux sorties, une par voie 
d’analyse. Les deux chaines de mesures sont utilisées 
simultanément. Les analyses par octave de fréquen- 
ces sont effectuées directement. L’enregistrement 
magnétique est utilisé en laboratoire pour des analy- 
ses fines du bruit. 


5.2.2. Relevés continus de la pression 
sonore 


L’étude précise de la répartition de la pression 
sonore sur une grande surface courbe telle que le 
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fuselage d’un avion, conduit a adopter le déplace- 
ment continu du microphone de mesures a distance 
constante de la structure: distance trés faible, de 
Yordre de quelques millimetres. C’est pourquoi, 
lorsque cet intervalle varie, on observe des varia- 
tions de pressions sonores qui s’expliquent par le 
phénoméne de doublage de pression lié a toute sur- 
face réfléchissante. 


Le dispositif expérimental réalisé rappelle par 
~son principe celui déja décrit a propos de |’étude du 
bruit des jets. Les microphones, cette fois-ci au 
nombre de deux, sont également portés par un cha- 
riot entrainé sur les éléments du rail du méme type. 
Cependant, la longueur des éléments du rail, ses 
supports, l’entrainement et la constitution du chariot 
sont différents. En effet, le rail se compose d’éléments 
de 0,5 m a 1 m de long afin de s’adapter a une sur- 
face courbe donnée. Les fourchettes réunissant les 
' éléments-du rail sont fixées sur la poignée des ven- 
_touses de miroitier, de sorte que le tout est directe- 
ment collé sur la structure (suivant une ligne de 
moindre courbure) grace a la succion exercée par 
les ventouses. 


Le chariot est entrainé par la rotation de la vis 
couplée a un moteur électrique. Le chassis du moteur 
se fixe sur la structure également au moyen de ven- 
touses. Le chariot est équipé de deux bras portant 
des tricycles. Au centre de chaque tricycle est fixé 
le microphone. Des rotules a deux degrés de liberté 

_ sont intercalées entre les bras et les tricycles afin que 
ces derniers puissent rouler sur une surface de cour- 
bure quelconque. Des ressorts de rappel appuient 
les tricycles sur la structure. Le chariot porte, en 
outre, les préamplificateurs cathodynes des micro- 
phones, l’inverseur de fin de course et un dispositif 
a microcontacteur qui commande la mise en marche 
automatique et l’arrét du chariot. 


Lors des essais, le chariot est équipé successive- 
ment de bras courts (Fig. 11), puis de bras longs. 
Pour une position donnée du rail, il est ainsi pos- 
sible de déterminer la variation des pressions sono- 
res suivant quatre lignes. Seule la mise en marche du 
chariot est télecommandée depuis le camion-labora- 
toire. Toutes les autres opérations telles que les 
arréts successifs, les déclenchements des analyses de 
pression sonore sont effectuées par le mécanisme 
automatique lui-méme, conformément au_ pro- 
gramme qui a été fixé avant les essais. 

La technique des enregistrements est identique a 
celle décrite a propos de l'étude du bruit des jets. 


5.3. Bruit causé par le survol des avions 


Le bruit causé par les avions a réaction survolant 
les zones résidentielles, riveraines des aérodromes, 
pose un probléme aigu en raison de la géne causée 
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Fig. 11. Chariot porte-microphones, bras courts. 


aux populations. Cela conduit a rechercher des pro- 
cédures de décollage donnant au sol lémission 
sonore minimum. Cette étude repose essentiellement 
sur les analyses de pression sonore au sol en fonc- 
tion du régime des turboréacteurs et de l’altitude 
des avions. On en déduit les courbes de variations 
qui permettent de définir la trajectoire de décollage 
appropriée a chaque type d’avions. 


La méthode de relevés d’évolution des niveaux 
sonores est particuliérement simple (Fig. 12). Au 


Fig. 12. Mesure de bruit lors d’un décollage «Cara- 
velle». 


moment du passage de l’avion a la verticale du point 
de mesure, les enregistreurs recoivent un «top» en- 
voyé dans les chaines de mesures par l’opérateur qui 
se tient prés du microphone. Un exemple d’analyse 
des nivéaux sonores par octave de fréquences est 
donné par la Fig. 13. 


Lorsqu’on dispose de deux chaines de mesures, il 
est bon de les utiliser simultanément avec des co- 
efficients d’amplification différents. Le cadrage du 
signal microphonique dans la bande dynamique de 
Penregistreur graphique est alors plus sir. 
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Fig. 13. Avionde transport 
a réaction. 
Trajectoire de décollage. 
Variations des niveaux so- 

nores au sol. 

Analyse des niveaux sono- 
res recus au sol lors d’un 
décollage. 
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6. Conclusions 


Les techniques de mesures acoustiques, décrites 
dans ce travail, ont été employées avec d’excellents 
résultats dans les conditions expérimentales ren- 
contrées dans le domaine de |’ Aéronautique. 

La streté des enregistrements des niveaux sono- 
res, la rapidité des opérations obtenue par |’écono- 
mie des moyens manuels, la simplicité des condi- 
tions d’emploi ont permis de réduire la durée des 
essais et le temps d’immobilisation des avions. 
L’abondance d’informations ainsi recueillies s’est 
avérée particuliérement précieuse pour la solution 
rapide des problémes acoustiques des avions SE 210 


«Caravelle». (Regu le 6 Décembre 1960.) 
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PoNemer me inhioCHiER KETTENLEITER ZUR UNTERSUCHUNG 
DER MECHANISCHEN SCHWINGUNGSVORGANGE 
IM INNENOHR 


von R. Orrincer* und H. Hauser ** 
Mitteilung aus dem Institut fiir Nachrichtentechnik der Technischen Hochschule Stuttgart 


Zusammenfassung 


Die Schwingungsvorgange im Innenohr, und hierbei besonders die der Basilarmembran, 
werden an Hand eines vereinfachten Modells beschrieben. Die dreidimensionalen Differen- 
tialgleichungen von Druck und Schnelle der Fliissigkeit in den Schneckenkandlen samt ihren 
Randbedingungen werden aufgestellt. Da letztlich nur die Auslenkung der Basilarmembran 
sowie der auf sie wirkende Druck von Interesse ist, kénnen die dreidimensionalen Diffe- 
rentialgleichungen zu eindimensionalen vereinfacht werden. Diese enthalten als einzige 
Variable nur noch die Ortskoordinate der Langsachse der Schneckenkanile. Die Krafte, die 
infolge der senkrecht zur Langsachse schwingenden Fliissigkeitsmasse auf die Basilar- 
membran wirken, konnen auch in den eindimensionalen Differentialgleichungen beriicksich- 
tigt werden, da diese Schwingung in erster Naherung gleichphasig mit der Basilarmembran 
erfolgt. Die Auflosung der eindimensionalen Differentialgleichungen wird mit einem elek- 
trischen Kettenleiter, der aus 65 Gliedern besteht, durchgeftihrt. Die MeBergebnisse am 
ausgefiihrten Kettenleiter werden mit bekannten Messungen, die direkt am menschlichen 
Gehor durchgefiihrt wurden, verglichen. Die Auslenkung der Basilarmembran bei Erregung 
durch Druckimpulse wird gemessen, und daraus die Geschwindigkeit der Wanderwelle auf 
der Basilarmembran ermittelt. Zuletzt wird noch die Auslenkung der Basilarmembran fiir 
den Fall untersucht, da sie ahnlich wie bei einem in der modernen Chirurgie gebrauch- 
lichen Verfahren, der sogenannten Fensterung, erregt wird. 


Summary 


Vibrations of the inner ear and especially of the basilar membrane are described with 
the help of a simplified model. Three-dimensional differential equations for pressure and 
particle velocity of the liquid in the spiral ducts are formulated together with suitable 
boundary conditions. Since only the pressure exerted on the basilar membrane and the 
distortions of the latter are of interest, the three-dimensional differential equations may be 
reduced to one-dimensional differential equations. The only variable is now the position on 
the axis of the spiral ducts. The vibration of the liquid normal to the axis of the spiral duct 
exerts forces on the basilar membrane, which may be included in the one-dimensional 
differential equations, since these oscillations are approximately in phase with the vibrations 
of the basilar membrane. The one-dimensional equations are solved with the help of an 
electric chain network consisting of 65 elements. Measuring results obtained with such 
networks are compared with results obtained from measurements in the human ear. The 
distortions of the basilar membrane are measured by exciting with a pressure impulse and 
the velocity of the travelling wave over the basilar membrane is found. Lastly the distortions 
of the basilar membrane are also studied for excitation not from the oval window but from 
another artificial opening, as used in modern surgery. 


Sommaire 


En utilisant un modéle simplifié on décrit les phénoménes oscillatoires dans Voreille 
interne et plus particuliérement ceux de la membrane basilaire. Les équations différentielles 
a trois dimensions relatives 4 la pression et a la vitesse du liquide dans les canaux coch- 
léaires sont établies en tenant compte de leurs conditions aux limites. Comme en fin de 
compte ce qui nous intéresse ce sont les élongations de la membrane basilaire et la pression 
exercée sur elle, nous pouvons en simplifiant ramener les trois équations a trois dimensions 
a une seule équation différentielle 4 une seule dimension. Cette derniére contient comme 
unique coordonnée l’abscisse curviligne comptée sur l’axe des canaux cochléaires. Cette 
équation différentielle 4 une seule dimension permet de tenir compte des forces exercées 
sur la membrane basilaire par les masses liquides qui oscillent perpendiculairement a l’axe 
longitudinal des canaux; ces oscillations sont en effet en phase, en premiére approximation, 
avec la membrane basilaire. On résout cette équation 4 une dimension a l’aide d’un réseau 
en chaine composée de 65 articles. Ce réseau une fois construit fournit des mesures que l’on 
compare a celles que donne directement |’étude de l’audition humaine. On mesure l’élonga- 
tion de la membrane basilaire excitée par des impulsions de pression et on en déduit la 
vitesse de l’onde transitoire sur la membrane basilaire. Enfin, on étudie encore les élonga- 
tions de Ja membrane basilaire dans un cas analogues a celui qui se présente dans une opé- 
ration de chirurgie actuellement pratiquée et qu’on appelle «Fensterung» (ouverture d’une 
fenétre) . 


* Jetzt Mitarbeiter der Telefunken G.m.b.H., Back- ** Jetzt Assistent am Institut fiir Theorie der Elek- 
nang. trotechnik an der Technischen Hochschule Stuttgart. 


162 


1. Einleitung 


Fir das Verstandnis des Hoérvorgangs ist die 
Kenntnis der Erregung der Sinneszellen von grofer 
Bedeutung. Der Sitz der Sinneszellen im Innenohr 
ist das Cortische Organ, das fest mit der Basilar- 
membran verbunden ist. Das Schwingungsverhalten 
der Basilarmembran bei Erregung durch ein Schall- 
ereignis ist somit bestimmend fiir die Erregung der 
Sinneszellen. 

Diese Arbeit beschrankt sich auf die Untersuchung 
des Schwingungsverhaltens der Basilarmembran, be- 
zieht aber die mechanischen Schwingungsvorgange 
des gesamten Innenohres in die Rechnung ein. 

Naheliegend sind direkte Untersuchungen im 
Innenohr, die aber auerst schwer durchfihrbar 
sind. Uns sind nur Untersuchungen von BEKEsy 
[1], [2] bekannt. Sie wurden mit Funkenknallen 
und mit Einzeltonen im Frequenzbereich bis etwa 
3kHz durchgefiihrt. Es ist jedoch wiinschenswert, 
einerseits das Schwingungsverhalten des Innenohrs 
fiir Einzeltone im gesamten Tonfrequenzbereich zu 
kennen, andererseits ist das Verhalten fiir beliebige 
Schallereignisse von Interesse. So hat zum Beispiel 
FLanaGan in einer kurzlich erschienenen Arbeit [11] 
aus den Messungen fiir Sinustone von BEKEsy mit 
Hilfe der Fouriertransformation das Impulsverhal- 
ten des Innenohrs berechnet. 

Uber die geometrischen Abmessungen des Innen- 
ohrs, uber das Massenverhaltnis der schwingenden 
Teile sowie tiber die Steifigkeit der Basilarmembran 
bei statischer Belastung [3] liegen umfangreiche 
und sichere MeBergebnisse vor. Deshalb ist es nahe- 
liegend, aus diesen Daten entweder mit Hilfe eines 
Modells oder mit mathematischen Hilfsmitteln auf 
das dynamische Verhalten des Innenohrs zu schlie- 
Ben. 

Die Herstellung eines mechanischen Modells, wel- 
ches das Schwingungsverhalten fiir den gesamten 
Tonfrequenzbereich richtig wiedergibt, ist nur be- 
dingt moglich [12], [13]. Der Grund hierfiir ist in 
erster Linie die grofe Steifigkeitsanderung der Basi- 
larmembran im Verlauf ihrer Lange. Ihre Nachbil- 
dung mit kiinstlichen Materialien sto8t auf erhebliche 
Schwierigkeiten. Bei der rein mathematischen Be- 
handlung ist es einfach, die Differentialgleichungen 
der Bewegungsvorgiinge aufzustellen [4], [5]. Ihre 
exakte Auflésung ist jedoch, auferst umfangreich. 
Zwistocki [4] hat vereinfachte Differentialgleichun- 
gen aufgestellt und ihre Auflosung fur rein sinus- 
formige Erregung durchgefitihrt. Die Auflosung die- 
ser von ZwisLock! aufgestellten Differentialgleichun- 
gen fiir nicht sinusformige Vorgange, etwa fiir Im- 
pulse, ist wiederum nur unter grof{em Rechenauf- 
_ wand mdglich [11]. Baucn [6] hat deswegen diese 
Differentialgleichungen mit Hilfe eines elektrischen 
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Netzwerkes dargestellt, und somit Losungen fiir be- 
liebige Schallvorginge gewonnen. Dabei zeigte es 
sich, da} die Ergebnisse die man auf Grund der von 
Zwistocki getroffenen Annahme erhielt, mit einigen 
Messungen von BExEsy [2] nicht ibereinstimmten. 
In dieser Arbeit werden deshalb Differentialglei- 
chungen der Bewegungsvorginge aufgestellt, denen 
weniger grobe Naherungen zugrunde liegen. Ihre 
Auflosung wird hier mit Hilfe eines elektrischen 
Netzwerkes durchgefihrt, da eine analytische Be- 
handlung auf zu groSe Schwierigkeiten stoft. 


2. Die Mechanik des Schneckenkanals und ihre 
Anniaherung durch ein vereinfachtes 
zweidimensionales Modell 


Wird das Trommelfell durch ein Schallereignis 
erregt, so wird dies tiber die Gehorknoéchelchenkette 
und die Steigbiigelplatte auf das ovale Fenster tiber- 
tragen. Das ovale Fenster ist der Abschlu8 des 
Innenohrs, das als schneckenformiger Kanal im Kno- 
chen eingebettet liegt. Bild 1 zeigt den Schnecken- 
kanal in schematischer Form als geradlinigen, recht- 


Basilar- 
membran 


Kanal] 
Kanal 2 


0) l 


Bild 1. Vereinfachtes -zweidimensionales Modell des 


Innenohres und seine Abmessungen. 


eckigen Kanal. Samtliche AuSenwande, abgesehen 
von der Vorderfront, stellen eine starre Knochen- 
wand dar. Der ganze Kanal ist mit der Perilymphe, 
einer Flissigkeit geringer Viskositat ahnlich der des 
Wassers, angefiillt. Der Schneckenkanal ist seiner 
Lange nach durch eine elastische Zwischenwand, die 
Basilarmembran, in einen oberen und unteren Kanal 
geteilt. In Wirklichkeit spannt sich die Basilarmem- 
bran nicht iiber die ganze Kanalbreite. Ein Teil der 
Zwischenwand wird durch eine feste Knochen- 
lamelle, die Schneckentrennwand, gebildet, an welche 
die Basilarmembran fest angewachsen ist. Die Zwi- 
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schenwand hort kurz vor dem Schneckenende auf. 
Kanal 1 und 2 sind durch dieses Loch, das Heliko- 
trema, miteinander verbunden. Am Schneckenein- 
gang sind die Offnungen beider Kanale, das ovale 
und das runde Fenster, durch Membrane abgeschlos- 
sen. 

Wird die Membran des ovalen Fensters in Schwin- 
gungen versetzt, so werden diese tiber die Lymph- 
flissigkeit auf die Basilarmembran tibertragen. Er- 
regt man zum Beispiel die Membran des ovalen 
Fensters mit einem kurzen StoB, so wird die Basilar- 
membran zunachst in der unmittelbaren Umgebung 
des ovalen Fensters ausgelenkt. Die Auslenkung 
pflanzt sich auf der Basilarmembran in Richtung 
der 2-Achse, das heiBt zum Helikotrema hin, wie 
eine Wanderwelle fort. Wie in dieser Arbeit gezeigt 
wird, nimmt die Geschwindigkeit dieser Wander- 
welle immer mehr ab, je mehr sie sich dem Heliko- 
trema nahert. Auf dem Weg entlang der Basilar- 
-membran andert sich die Form der Wanderwelle in 
der Weise, da die zeitliche Auslenkung der Basilar- 
membran um so langsamer verlauft, je mehr sich 
_ die Wanderwelle dem Helikotrema nahert. 

Die Nachgiebigkeit der Knochenwand und die 
Kompressibilitat der Lymphfliissigkeit ist gegentiber 
der Nachgiebigkeit der beiden Fenstermembrane 
vernachlassigbar klein. Eine Auslenkung der Mem- 
bran des ovalen Fensters mu sich deshalb voll iiber 
die Membran des runden Fensters ausgleichen. 

Die Darstellung in Bild 1 gibt aus Griinden der 
mathematischen Behandlung die tatsachlichen Ver- 
haltnisse sehr vereinfacht wieder. Die Querschnitts- 
form wurde zu einem Rechteck vereinfacht, der 
Schneckenkanal seiner Lange nach zu einem gerad- 
linigen Kanal aufgebogen. In Bild 1 zeigen die bei- 
den Skizzen links unten den Verlauf der Hohe der 
Kanale nach [4]. Die in den Skizzen eingezeichnete 
Gerade, stellt fiir beide Kanile eine gute Naherung 
dar. Der daraus folgende symmetrische Aufbau ver- 
einfacht die mathematische Behandlung. Eine wesent- 
lich einschneidendere Vereinfachung ist die bereits 
vorher erwahnte Annahme, daf sich die Basilarmem- 
bran tiber die ganze Breite der Kandle spannt. Fer- 
ner ist angenommen, dafi die Kanalbreite 6 konstant 
ist. Nimmt man weiterhin an, dafi die Auslenkung 
» der Basilarmembran fiir konstantes 2 von z unab- 
hangig ist, das heif®t, da die Basilarmembran in 
ihrer gesamten Breite gleichmafig schwingt, so ist 
das an sich dreidimensionale Schwingungsproblem 
auf ein zweidimensionales zuritickgefiihrt. Samtliche 
auftretenden Groen, der Flissigkeitsdruck, die 
Geschwindigkeit der Fliissigkeitsteilchen sowie der 
Basilarmembran, sind nur Funktionen von x und y , 
nicht aber von z. 

Die Vereinfachung der Abmessungen und der 
Mechanik des Innenohres durch das zweidimensio- 
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nale Modell in Bild 1 stellt eine recht grob erschei- 
nende Naherung dar. Bei den bisherigen Lésungs- 
versuchen [4] ging die Vereinfachung jedoch noch 
einen Schritt weiter, indem das Schwingungsproblem 
schon im Ansatz auf ein eindimensionales reduziert 
wurde. 

Dem endgiiltig ausgeftihrten elektrischen Ketten- 
leiter wird nicht das zweidimensionale Modell von 
Bild 1 sondern das vereinfachte dreidimensionale 
Modell nach Bild 4 zugrunde gelegt. Da die mathe- 
matische Behandlung des dreidimensionalen Modells 
wesentlich umfangreicher und untbersichtlicher ist, 
wird die Behandlung des zweidimensionalen Modells 
vorangestellt. Das dabei gezeigte Losungsverfahren 
laBt sich ohne weiteres auch auf das vereinfachte 
dreidimensionale Modell anwenden. 


3. Die Differentialgleichungen der Bewegungs- 
vorgange im zweidimensionalen Modell und ihre 
angepaBbte Form zur Auflosung durch einen 
elektrischen Kettenleiter 


Die Flissigkeitsbewegungen in den beiden Schnek- 
kenkanalen sind fiir die Auslenkung der Basilar- 
membran mitbestimmend. Bei der mathematischen 
Behandlung miissen daher die Differentialgleichun- 
gen fur die Flussigkeitsbewegungen aufgestellt wer- 
den. Bei ihrer Auflosung treten die Impedanzen der 
Basilarmembran und der beiden Fenstermembrane 
sowie der Knochenwand als Randbedingungen in 
Erscheinung. Wie vorne erwahnt, wird keine voll- 
standig analytische Aufldsung dieser Differential- 
gleichungen angestrebt. Sie werden vielmehr unter 
Beriicksichtigung der Randbedingungen und erlaub- 
ter Naherungen auf ein System von linearen Diffe- 
rentialgleichungen zuriickgefuhrt, die durch einen 
verhaltnismaBig einfachen Kettenleiter dargestellt 
werden konnen. Dies bedingt zweckmafigerweise 
eine Form der Darstellung, die von dem sonst in 
der Theorie der Differentialgleichungen tblichen 
Losungsgang abweicht. Die erhaltenen Differential- 
gleichungen beschreiben nicht mehr die Vorgange im 
gesamten Modell, sondern verkntipfen im wesent- 
lichen nur noch die wirklich interessierenden Gro- 
fen, némlich die Auslenkung der Basilarmembran 
und den Druck der Flissigkeit unmittelbar auf der 
Oberflache der Basilarmembran. 


3.1. Die zweidimensionalen Differentialgleichungen 
der Flissigkeitsbewegung und thre Riickfiih- 


rung auf eine eindimensionale Form 


Das Verhalten von Druck und Schnelle in einer 
Flissigkeit wird durch die Kontinuitatsgleichung 
und die Euterschen Bewegungsgleichungen vollstan- 
dig beschrieben. Die Kontinuitatsgleichung fiir in- 
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kompressible Flissigkeiten hat bei zweidimensiona- 
len Vorgangen die Form [7] 


—=+354=0. (1) 


Dabei sind v, und v, die Komponenten der Ge- 
schwindigkeit der Fliissigkeitsteilchen in z- bezie- 
hungsweise y-Richtung. Physikalisch besagt die Kon- 
tinuitatsgleichung, da} die in der Zeit dé in ein 
Raumelement einstromende Flissigkeitsmenge gleich 
der in derselben Zeit dt ausstromenden Flissigkeits- 
menge ist. 

Die Evterschen Bewegungsgleichungen haben fiir 
kleine Schwingungsamplituden (bei zweidimensio- 
nalen Vorgangen) die Form [7] 


Op Ov, on 
Be UME ae On ence 
p) ) 


Dabei sind p der Flissigkeitsdruck, @ die Dichte der 
‘Flissigkeit und ry der mechanische Reibungswider- 
stand der Flissigkeit je Volumeneinheit. Physika- 
lisch besagen die Euterschen Bewegungsgleichungen, 
da das Produkt aus Masse und Beschleunigung 
eines Fliissigkeitsteilchens gleich der Resultierenden 
der angreifenden Krafte ist. 

Fir die Auflosung durch einen Kettenleiter ist es 
zweckmaBig, sich auf sinusformige Vorgange zu be- 
schranken. Da die Losung im Kettenleiter fiir alle 
Frequenzen giiltig ist, bedeutet dies keine Einschran- 
kung der Allgemeinheit. Die Gln. (2a) und (2b) 


vereinfachen sich dadurch zu 


Op Zor =O. (ora) op +z9vy=0. (3b) 
Ox Oy 
Dabei ist 
Z=M+joe (4) 


der komplexe Reibungswiderstand je Volumenein- 
heit. Die Groen p, v; und vy sind jetzt keine Augen- 
blickswerte mehr, sondern komplexe Groen. Wegen 
einer einfacheren Darstellung werden dafiir keine 
neuen Symbole eingefthrt. Die Kontinuitatsglei- 
chung (1) kann daher unverandert fir die sinus- 
formigen Vorginge tibernommen werden. 

Die Gln. (1) und (3) beschreiben das Verhalten 
von Druck und Geschwindigkeit an jeder Stelle x, y 
des Kanals. Interessant ist in erster Linie die Aus- 
lenkung 7 der Basilarmembran beziehungsweise ihre 
Auslenkungsgeschwindigkeit 7 . In den Gln. (1) und 
(3) bedeutet dies, dafS bei der Geschwindigkeits- 
komponente v,(z,y) nur der Wert an der Stelle 
y =0, namlich v, (x, 0) von Interesse ist. 

Mit der in Bild 1 gewahlten Komponentenrich- 
tung ergibt sich 

vy (2,0) = (2). (5) 
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Ebenso wird fiir den Druck nur die Loésung fir 
p(x,0) vorgenommen. Die folgenden Uberlegungen 
werden zunachst nur fir Kanal 1 durchgefiihrt. Am 
Schlu8 des Abschnittes werden dann die entspre- 
chenden Gleichungen fiir Kanal 2 angeschrieben. 

Integriert man Gl. (1) tuber y von y=0 bis 
y= +h(x), das hei®t von der Basilarmembran bis 
zur gegentiberliegenden Knochenwand von Kanal 1, 
so ergibt sich 

h (@) 


[PEP ey soy h) —v,(x,0) =0. 
Ox 
rt 


Da_starre Knochenwand vorausgesetzt wird, ist 
v,(x,h) =0. Durch Vertauschen der Reihenfolge 
von Integration und Differentiation erhalt man aus 
der vorigen Gleichung 

h (a) 


a [ec(ew) dy — vz (x, h) dh —v,(x,0) =0. 
Ox, da 


0 
ZweckmaBigerweise fihrt man hier die tber den 
ganzen Kanalquerschnitt gemittelte Geschwindig- 
keitskomponente in x-Richtung 

h (x) 


% 1 
Bea) = 5. [ relay) dy (6) 
0 
ein. 
Nach einigen Umformungen erhalt man damit 


dv, (x) 


h(x) { + © [Be(2) — ve (a, h) vy (2, 0) =0. 


(7) 
Fir den in Bild 1 unten eingezeichneten Verlauf 
h(a) gilt dh/de=0,02 <1. Da die Viskositat der 
Lymphflissigkeit verhaltnismaBig gering ist, wird 
v, an der Knochenwand nicht wesentlich verringert. 
Damit ist 0,(2) und v;(2, h) von derselben Grofen- 
ordnung, so da das zweite Glied in Gl. (7) ver- 
nachlassigt werden kann. Mit Gl. (5) erhalt man 
damit die gesuchte Form fiir die Kontinuitatsglei- 
chung im Kanal 1: 


waza 


h(a) P22) 4 (2) <0. (8) 
x 


Um die Kontinuitatsgleichung in dieser Form an- 
schaulich interpretieren zu konnen, betrachtet man 
einen Scheibe, die durch zwei zur x-Achse senkrecht 
stehende Schnitte an der Stelle x und x + dx aus dem 
oberen Kanal ausgeschnitten wird. Die Gl. (8) be- 
sagt nun, daB die an der Stelle x in die Scheibe ein- 
tretende Fliissigkeitsmenge gleich der Fliissigkeits- 
menge ist, die einmal an der Stelle x+dx aus der 
Scheibe austritt, und die zum andern, durch die 
Bewegung der Basilarmembran bedingt, nach unten 
aus der Scheibe austreten kann. 

Die Bewegungsgleichungen (3) werden in 4hnli- 
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cher Weise zunachst fiir Kanal 1 so umgeformt, dah 
nur noch die Groen 7 (x), 0,(x) und p(x, 0) auf- 
treten. Es ware naheliegend, analog der mittleren 
Geschwindigkeit v,(x) einen mittleren Druck p(x) 
einzuftihren. Dieses Verfahren scheidet jedoch aus, 
da wir die auf die Basilarmembran wirkende Druck- 
differenz brauchen. 

Integriert man Gl. (3b) tiber y von 0 bis y und 
differenziert anschlieBend partiell nach x, so erhalt 


“man 


Op (x,y) 
ar 


Sp(x,0) 


y 
vy fee fon dy=0. 
0 


(9) 
Eliminiert man aus den Gln. (3 a) und (9) das erste 
Glied, so erhalt man eine einzige Bewegungsglei- 
chung: 

-Op(x, 0) 


Ox 


y 
- + 29 Vz (2, Yy) <a fetey) dy. 
6 (10) 
Durch Mittelung von Gl. (10) tiber y kann v,(z, y) 
gemas Gl. (6) durch die mittlere Geschwindigkeit 


_ 0, (az) ersetzt werden: 


(11) 
h (2) y 


| vy (054) dy,| dy =0. 


fe. 


Die ersten beiden Glieder der Gl. (11) enthalten 


nur noch die gewunschten GrofSen. Das dritte Glied 


enthalt die GroBe v, (x, y). Der weitere Losungsgang 
wird wesentlich dadurch bestimmt, daf fiir v, (x, y) 
folgende Annahme getroffen wird: 


Vy (2, Y) = vy (x, 0) Fly). (12) 


| Dabei soll F(y) eine von x unabhingige Funktion 


ICO CCTCDCEL Ltt Et TE EE 


sein. Bei sinusformigen Vorgangen soll F(y) zudem 
eine reelle GroBe sein, die von der Frequenz unab- 
hangig ist. Anschaulich ausgedrickt bedeutet diese 
Annahme, dafs die Flissigkeit senkrecht tiber der 


“¥=0 


Bild 2. Schematische Darstellung des zeitlichen Ge- 
schwindigkeitsverlaufes der Kanalflissigkeit 
senkrecht zur Langsachse bei Erregung durch 
einen DruckstoB. Die ausgezogene Kurve stellt 
den Geschwindigkeitsverlauf der Basilarmem- 
bran an einer beliebigen Stelle x dar. Die ge- 
strichelten Kurven zeigen die Geschwindigkeit 
der Fliissigkeitsteilchen an derselben Stelle x 
iiber der Basilarmembran. 
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Basilarmembran gleichphasig mit dieser in y-Rich- 
tung schwingt. Hat etwa bei einer impulsformi- 
gen Erregung die Auslenkungsgeschwindigkeit 
= —v,(t, y=0) der Basilarmembran den in Bild 2 
schematisch aufgezeichneten Verlauf (ausgezogene 
Kurve), so soll vy(t,y) an einer Stelle y den in 
Bild 2 gestrichelt aufgezeichneten Verlauf haben. 
Der zeitliche Verlauf der Geschwindigkeit v, an einer 
beliebigen Stelle y unterscheidet sich also von 7 nur 
um den Faktor F(y), der von t und x unabhangig 
ist. Nach Gl. (12) ist F(0) =1 und, wegen der star- 
ren Knochenwand, F(h) =0. 

Die Berechtigung der oben getroffenen Annahme, 
die durch Gl. (12) ausgedriickt wird, laBt sich durch 
einen Vergleich mit den Untersuchungen iiber Kanal- 
wellen von Lams [8] einsehen. Da sich die dort ge- 
gebenen Verhaltnisse nicht ohne weiteres auf den 
Schneckenkanal iibertragen lassen, wurde von Hav- 
SER in einer unverOffentlichten Arbeit eine Abschat- 
zung der infolge dieser Annahme moglichen Fehler 
vorgenommen, Dabei zeigte es sich, da die auftre- 
tenden Fehler zu yernachlassigen sind. 

Uber den Verlauf der Funktion F(y) wurden bis 
jetzt keine Voraussetzungen getroffen. Vertauscht 
man im dritten Glied von Gl..(11) die Differentia- 
tion mit der Integration, so erhalt man wieder unter 


Beriicksichtigung, dap <1 ist, fiir diesen Aus- 
iS 


druck: 


h(w) y ; 


“02 lays.) f [Pod andy]. 03) 
Ow 20) 


Fir die weitere Rechnung geht demnach F(y) nicht 
als eine Funktion von y ein, sondern lediglich der 
Zahlenwert, der sich durch zweimalige Integration 
gemaf} GI. (13) ergibt. Nach den Untersuchungen 
von Lamp [8] bietet sich hier ftir F(y) als Nahe- 
rung der in Bild 3 dargestellte Verlauf an. Damit 


1 


js 
0 


h 


0 y 
Bild 3. Das Verhaltnis F der Geschwindigkeit v,(y) der 
Fliissigkeitsteilchen an einer Stelle x, bezogen 
auf die Geschwindigkeit v,(0) der Basilarmem- 
bran, in Abhangigkeit von y . 


lat sich der Integralausdruck auswerten, und man 
erhalt fir Gl. (11) mit Gl. (5) als neue Bewegungs- 
gleichung fiir Kanal 1: 


Pee lee Oye aie) \i.0 oy (14) 
da 3 


Um die Bewegungsgleichung in dieser Form an- 
schaulich interpretieren zu konnen, betrachtet man 
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wiederum eine Scheibe im Kanal 1] an der Stelle x 
mit der Dicke dx. Das erste Glied in Gl. (14) stellt 
die durch Reibung und Massebeschleunigung an der 
Scheibe in 2-Richtung angreifende Scheinkraft dar. 
Das zweite Glied stellt die z-Komponente der Kraft 
dar, die durch Druckanderung entlang der x-Achse 
auftritt. Hatte man analog der mittleren Geschwin- 
digkeit 0,(x) den mittleren Druck 
h(a) 


pie) = 5 [ply ay 
0 


eingefilhrt, so hatte man statt Gl. (14) direkt aus 
der Bewegungsgleichung in 2-Richtung (3a) die 


Beziehung 


zy0_(z) + PW). 9 
dx 
erhalten. 
Der Sinn dieser umfangreichen Rechnung aft sich 


also darin zusammenfassen, da der mittlere Druck 
(p(x) durch den Ausdruck p(x, 0) + “2 n(a) ersetzt 


wurde, das heiBt durch den Druck p(z,0) an der 
Basilarmembran und deren Geschwindigkeit 7 (2). 
Es sei an dieser Stelle erwahnt, daB in den friheren 
Arbeiten [4], [5] die Abhangigkeit von y véllig 
vernachlassigt wurde, so daf} sich die hier durch- 
gefuhrte Rechnung eribrigte. In unseren Grofen 
ausgedriickt bedeutet dies, dai p(v) =p(x,0) ge- 
setzt wurde. 

Denkt man sich einen Schnitt senkrecht zur 
x-Achse durch die beiden Kanale in Bild 1 gelegt, 
so ist wegen der starren Knochenwand die Summe 
der Flissigkeitsmengen, die durch diese Schnittflache 
in die beiden Kanale einstromen, gleich Null. 

Dies bedeutet 


+h 0 h 
[r(x y) dy = fre(x,y) dy + [ vz(x, y) dy=0. 
—h —h 0 (15) 


Da im folgenden die Vorgange im Kanal 1 und 2 
gleichzeitig betrachtet werden, setzen wir zur Unter- 
scheidung die Indizes 1 und 2. Setzt man analog zu 


Gl. (6) 


0 
1 iy 
j | oy = 8222). (16) 
—h 
ferner zur Vereinfachung 
ile) ey, 
Vz,9(2) =v2, (17) 


Py (x, 0) = Pils 

P(x, 0) = po; 
so erhalt man aus-Gl. (15) als zusatzliche Konti- 
nuitatsgleichung 


Vy +o=0. (18) 
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Die zur Kontinuitatsgleichung (8) und zur Bewe- 
gungsgleichung (14) entsprechenden Gleichungen 
fiir Kanal 2 erhalt man, indem man die Integratio- 
nen analog von —A(x) bis 0 durchfiihrt. Mit den 
Bezeichnungen (17) ergeben sich damit 

pool women. (19a) 
fur Kanal 1 me 


d Ry lee ; 
z = Dj |. = Oe 
OES a [P+ 3 i| ( } 
ne = 0: (20 a) 
fiir Kanal 2 ns 4 
Zq 0p b= (Ps ee i) =0. (20b) 
dx 3 


Die Gln. (19) und (20) beschreiben die Vorgange 
fiir beide Schneckenkanile getrennt. Diese Vorgange 
bestehen jedoch nicht unabhangig voneinander, son- 
dern werden vielmehr durch die zusatzliche Konti- 
nuitatsgleichung (18), die Impedanz der Basilar- 
membran sowie durch weitere Randbedingungen 
miteinander verkniipft. 


3.2. Die Impedanz der Basilarmembran und das 
Verhalten am Helikotrema als Randbedingung 


Bei dem im Abschnitt 2 eingefiihrten Modell nach 
Bild 1 gilt fiir die Bewegung der Basilarmembran 
an der Stelle x die Beziehung [4] 
dy dy Sew 
Stig ee b et Leg 
B diz (Py — P2) 6g — TB dia 
Dabei ist 


my die Masse der Basilarmembran pro Langen- 
einheit, 

rg der Reibungswiderstand der Basilarmembran 
pro Lingeneinheit, 

cp die Nachgiebigkeit der Basilarmembran pro 
Langeneinheit und 

bp die Breite der Basilarmembran. 


Nach Bild 1 ist fiir das hier behandelte verein- 
fachte zweidimensionale Modell by, = 6, wahrend in 
Wirklichkeit die Breite der Basilarmembran kleiner 
als die Kanalbreite ist. 

Um Gl. (21) anschaulich interpretieren zu kén- 
nen, betrachtet man in Bild 1 einen Streifen der 
Breite dx der Basilarmembran an der Stelle 2. Das 
Produkt aus der Masse dieses Streifens mg dx und 
der Beschleunigung 4) ist gleich der Summe der auf 
ihn wirkenden Krafte. GemaB der rechten Seite von 
Gl. (21) setzt sich diese Summe aus folgenden Teil- 
kraften zusammen: 

1. Aus der Kraft (p,—p,:) bg dx, die von der 

Druckdifferenz der Fliissigkeit zwischen obe- 
rem und unterem Kanal herriihrt. 


2. aus der Reibungskraft rg dx, die durch die 


m 


(21) 
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innere Reibung der Basilarmembran verursacht 
wird, und 


3. aus der elastischen Riickwirkkraft dz, die 
cB 
durch die Auslenkung der Basilarmembran ent- 
steht. 
Beschrankt man sich wiederum auf sinusférmige 
Vorgange, so erhalt man aus Gl. (21) die Beziehung 


: Ea 
(Pi — Pe) ba = [iomtret 2). (22) 
Joc, 


Definiert man die Impedanz der Basilarmembran je 
Langeneinheit durch 


=e oles oe (23) 


2B 


das heift durch den Quotienten aus Druckdifferenz 
und Schnelle, so erhalt man damit aus Gl. (22) 


(24) 


2B= (rai mn + : : ) 
bz joc 
Die GI. (23) stellt zu den Gln. (19) und (20) eine 
weitere Beziehung dar. 

Wie sich spater zeigen wird, treten in den Diffe- 
rentialgleichungen fur das elektrische Netzwerk Gro- 
Ben auf, die zu den Geschwindigkeiten v, und vy und 
den Fliissigkeitsdriicken p, und p, in einer analo- 
gen Beziehung stehen. Man muf daher aus den vor- 
handenen Gleichungen die Schnelle der Basilarmem- 
bran 7 eliminieren. Mit Gl. (23) ergibt sich aus den 
Gln. (19) und (20) 


i dv, 


“he ee ae (25 a) 

fur Kanal 1 . ‘ te ; : 
2g Vyr d bi aid (Pr P| 0, 

oe 3 

(25 b) 
pave eure lL 20: (26 a) 

fir Kanal 2 4 ; Le : 
Zo Va [p. 26 (Pe ps =0. 

daz 2B 
(26 b) 


Die Gln. (25) und (26) enthalten nur noch die Gro- 
fien, die fur den Anschlu8 an das elektrische Netz- 
werk notig sind. Allerdings enthalten die Gleichun- 
gen fiir Kanal 1 noch die GroBe py und umgekehrt. 
Die weiteren Umformungen bestehen nur noch aus 
der vollstandigen Trennung der Variablen von Ka- 
nal | und 2. 

Addiert man Gl. (25b) und (26b), so wird 
unter Berticksichtigung von Gl. (18) 


dp; , dps 
—++—2 =0. 
dx dx 
Integriert man diese Gleichung tiber x, so erhalt man 
Py +P2=2puH- (27) 
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Da am Helikotrema p,=py, ist, stellt die Integra- 
tionskonstante py den Flissigkeitsdruck am Heliko- 
trema dar. Fiihrt man noch die Abkiirzung 


O,, pee (28) 


ein, so ergibt sich mit Gl. (27) die endgultige Form 
der Kontinuitatsgleichungen und Bewegungsglei- 


chungen 
d ” 
ee hs be) 9s. <9) 
fur Kanal 1 a 
2903 +——[(1 +4) (p:—pu)] =9, 
dx 
(29 b) 
dvs 2 
at pay Pe PM =0, (80a) 
fiir Kanal 2 se 


fo 2 +4 [(1 +4) (P2~ Pu) ]=9. 
(30 b) 


Ferner gilt noch die bereits erwahnte Beziehung 


Vy +Vg=0. (18) 
Die Gln. (29) und (30) sind zur Darstelluag durch 
einen elektrischen Kettenleiter besonders geeignet. 
Sie sind vollig gleichartig aufgebaut, was wegen der 
Symmetrie des Modells auch zu erwarten war. 


3.3. Verbesserung der gefundenen Losung mit einem 
vereinfachten dreidimensionalen Modell 


Bei der bisherigen Rechnung wurde aus Griinden 
der einfacheren mathematischen Behandlung an- 
genommen, da} sich gema dem Modell von Bild 1 
die Basilarmembran tiber die ganze Breite des 
Schneckenkanals spannt. Das an sich dreidimensio- 
nale Problem konnte dadurch auf ein zweidimensio- 
nales zurickgefiihrt werden. Im Verlauf der Rech- 
nung wurde weiterhin angenommen, dai die senk- 
recht zur Kanalachse schwingende F'liissigkeit, tiber 
den ganzen Kanalquerschnitt gesehen, gleichphasig 
mit der Basilarmembran schwingt. Mit dieser An- 
nahme wurden die zweidimensionalen Differential- 
gleichungen auf eindimensionale zuriickgefiihrt. Man 
kann der Rechnung auch ein vereinfachtes drei- 
dimensionales Modell, wie in Bild 4 aufgezeichnet, 
zugrunde legen, bei dem sich die Basilarmembran 
nur tber einen Teil der Kanalbreite spannt. Behalt 
man die Annahme bei, da die Flissigkeitsteilchen 
an einer Stelle x im ganzen Kanalquerschnitt in y- 
und z-Richtung gleichphasig schwingen, so laBt sich 
dieses dreidimensionale Problem unmittelbar auf ein 
eindimensionales zuriickfiihren. Mit einem ahnlichen 
Rechnungsgang lassen sich auch hier fiir die inter- 
essierenden Groen eindimensionale Differential- 
gleichungen aufstellen, die gleich aufgebaut sind wie 
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Helikotrema A 


Basilarmembran 


-h 
l= 325mm 
125mm jesmm  somm 
| VA 0,6mm [oY 
Sa eb . 
h 0 H ar : ! 
0 l 0 i 
t AS5M™M 10mm 
{ |06mm wgee = STE 
is vA by y 
Bg 2 ! 
0 l 0 l 
x 
Bild 4. Vereinfachtes dreidimensionales Modell des 


Innenohres und seine Abmessungen. 


die im Abschnitt 3.2. abgeleiteten Differentialglei- 
chungen (29) und (30). Ihre recht umfangreiche 
Herleitung wurde in einer unveroffentlichten Arbeit 
von Hauser durchgefthrt. Eine kurze Zusammen- 
fassung dieser Arbeit wird im Anhang gegeben. [hr 
Ergebnis lautet: 


dv; 2 bp Zz ( 

Sr (Pa PH) 0s 
Fir Kanal 1 dx by hz 
bzw. 2 d 

Zo viet a [(1+ 8) (pie—pu)]=0. 

i (A 9) 
Ferner gilt auch hier die Beziehung (18) : 
V4 +2.=0. 


2 ' 
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Dabei wurden die in den Abschnitten 3.1. und 3.2. 
eingefihrten Groen in derselben Bedeutung iiber- 
nommen. 

Nach Bild 4 ist ferner 

bp die Breite der Basilarmembran, 

bx die Breite des gesamten Schneckenkanals. 


Weiterhin ist 
f= (+13 32)2 ! (s+13 58) 22. 
bx ZB 
10) 


16 bx 24 
Fur den zweidimensionalen Fall ist bg/bk=1. Es 
ist dann =a, und damit gehen die Gln. (A9) in 
die Gln. (29) und (30) iiber. Wie aus Bild 4 hervor- 
geht, ist by/bx eine Funktion von z, die stets kleiner 
als 1] ist. Fir die Auflosung durch einen elektrischen 
Kettenleiter sind die Gleichungssysteme (29), (30) 
und (AQ) vollig gleichwertig. Man erhalt in beiden 
Fallen dieselbe Netzwerksanordnung, lediglich die 
einzelnen Schaltelemente haben teilweise verschieden 
groke Werte. Da die Gln. (A9) den tatsachlichen 
Verhaltnissen naher kommen, wird die Auflésung 
fiir sie durchgefiihrt. Lediglich aus Griinden der bes- 
seren Ubersichtlichkeit, wurde die einfachere Ablei- 


tung ftir das zweidimensionale Modell vorangestellt. 


4. Darstellung der Bewegungsvorgange durch 
einen elektrischen Kettenleiter 


4.1. Die Netzwerksanordnung mit zundchst noch 
unbestimmten Zweipolen 


Die Differentialgleichungen (AQ9) beziehungs- 
weise (29) und (30) lassen sich durch einen elektri- 
schen Kettenleiter mit konzentrierten Schaltelemen- 


ten nicht exakt beschreiben. Auf Grund der end- 


———<$<$<—<_— 4 ie 


Xy-1 


xy Xe Xn=l Xx 


Bild 5. Der elektrische Kettenleiter zur Darstellung der Bewegungsvorgange im Innen- 
ohr, aufgebaut aus allgemeinen Zweipolen. 


i 
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lichen Anzahl von Schaltelementen konnen mit 
einem solchen Kettenleiter naturgema8 nur Diffe- 
renzengleichungen dargestellt werden. Wahlt man 
die Differenzen hinreichend klein, so konnen damit 
die .Differenzengleichungen in Differentialgleichun- 
gen tbergefthrt werden. 

Stellt man von dem in Bild 5 aufgezeichneten Ket- 
tenleiter die Knoten- und die Maschenregel an der 
Stelle x, auf, so erhalt man fiir den Index 1, bezie- 
hungsweise 2 die Differenzengleichungen 
1. Knotenregel : 


Aly 2(%) = ye (ay+1) —I12 (a) 


= — (U12(%) —Ux) Yq(m), (31 a) 
2. Maschenregel: 
A[ (1+ 6 (%)) Ui2(%) — Ux) ] 
= (14 B(a+1)) (U12(a@+1) — Ug) 
= (1+8(%)) (U12(%) — Ux) 
= —Siio (aya) Zy. (31 b) 


Dabei sind U,(x,) und U,(x,) die Spannungen vom 
Punkt At, beziehungsweise Az,, zum Bezugspunkt 
M, J,(%) und /,(x,) die Strome im oberen bezie- 
hungsweise im unteren Langszweig. 

Aus Symmetriegriinden ist 


I, (a) +Ip(a) =0. (32) 


Der Strom im Mittelleiter ist also Null. In der schal- 
tungsmaBigen Ausfiihrung kann daher der Mittel- 
leiter wegfallen. 

Laf®t man die Abstande Ax=2,+1 +2, sowie die 
Groen Z; und Y,(x,) in der Weise immer kleiner 
werden, daf} die Quotienten 


4 


Gena Yq (a) 
Ax 


und’ Y,(z) G29 

n al ) Re 

konstant bleiben, so gehen unter Vernachlassigung 
der GroRen 2. und hoherer Ordnung die Differenzen- 
gleichungen in Differentialgleichungen tber: 


_ Knotenregel: 


Siete) + Yq (2) (Uy2(e) —Un)=0, (33a) 


Maschenregel: 
Zilia(e) +“ [1 +A(2)) (Us2(@) —Un)] =9. 
(33 b) 
Vergleicht man die Gln. (A 9) mit den Gln. (33), so 
sieht man, da} die mechanischen und die elektrischen 
Differentialgleichungen dieselbe Form haben. Dabei 
steht an Stelle der mittleren Schnelle v;2 der Strom 


J,2 und an Stelle des Fliissigkeitsdruckes p; 9 die 
Spannung U;». Fiir den Druck py am Helikotrema 


_ steht die Spannung Uy. Verwendet man zwischen 


den mechanischen und elektrischen Groen die Ana- 
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logiebeziehungen 


Pi2g=Ny Ux2, 
34 
De dts | oe 
Ny 
so ergeben sich durch Koeffizientenvergleich der 
Gln. (A9) und (33) zwischen den mechanischen 
und elektrischen Grofien die Beziehungen 


=, (35) 
Np Ny 
= 2b 
a= 5 5 36 
barred bx h zp ( ) 


Gibt man der Gro8e f(x) im Kettenleiter den Ver- 
lauf gema Gl. (A 10), so erhalt man als weitere 
Beziehung 


(37) 


S| 5S! 


Z,, Yq und Y,/f sind demnach komplexe Zweipole. 
Ihre Frequenzabhangigkeit und ihre Darstellung 
durch Widerstande, Spulen und Kondensatoren, 
wird im nachsten Abschnitt behandelt. 

Beim schaltungsmafigen Aufbau des Kettenleiters 
konnen mit Riicksicht auf die Anzahl der Schalt- 
elemente die Abstaénde Az nicht beliebig klein ge- 
wahlt werden. Bei einer Unterteilung in endlich viele 
Abschnitte, haben dann die Schaltelemente die end- 
lichen GroBen 


und (38) 


Z, = Ue Az 
Y,=Y, Az 
Nach den Untersuchungen von Baucu [6] geniigt es, 
den Abstand 4x=0,05 cm zu wiahlen. Die Lange / 
der Basilarmembran wird demnach in 64 Abschnitte 
unterteilt. Bei dieser Unterteilung ist die GroBe Ax 
in der GroS®enordnung der Querschnittsabmessungen 
der Kanadle. Wie aus den Untersuchungen [6] be- 
kannt ist, erfahrt ein Schallereignis bei seiner Fort- 
pflanzung langs der Basilarmembran auf einer sol- 
chen Wegstrecke Ax keine nennenswerte Formande- 
rung. Der so aus konzentrierten Schaltelementen 
aufgebaute Kettenleiter gibt also die kontinuierli- 
chen mechanischen Vorgange mit gentigender Ge- 
nauigkeit wieder. 

Der Aufbau des Netzwerkes in Bild 5 wurde so 
vorgenommen, daf} ein anschaulicher Vergleich mit 
dem Modell leicht méglich ist. Die obere Halfte des 
Netzwerkes entspricht dem Kanall, die untere 
Halfte dem Kanal 2. Die Stelle x=0 entspricht dem 
Ort des ovalen und runden Fensters. Um den Hor- 
vorgang richtig nachzubilden, mu, wie in Bild 5 
eingezeichnet ist, zwischen P, und M, das heifit am 
ovalen Fenster eingespeist werden. Der Innenwider- 
stand Z, entspricht der mechanischen Impedanz der 
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Steigbiigelplatte und der Membran des ovalen Fen- 
sters. Die Spannung Uy entspricht der erregenden 
Kraft auf die Steigbiigelplatte. Der AbschluBwider- 
stand Z, entspricht der Impedanz der Membran des 
runden Fensters. 

GemaB8 der Analogiebeziehung (34) kann durch 
Spannungsmessung zwischen den Punkten A1,, und 
M beziehungsweise Az, und M der Fliissigkeitsdruck 
an der oberen und unteren Seite der Basilarmem- 
bran abgelesen werden. Die Spannung zwischen den 
Punkten Ai, und Ag, entspricht der Druckdifferenz 
Pi — P2=Np(U, — Us). : 

Setzt man nach Gl. (23) p,—po=%p3% und nach 
Bild5 U,—U,= —2AI/Y,, so erhalt man mit Gl. 
(36) und (38) die Beziehung 


(39) 


Der Strom durch einen Querzweig an der Stelle x, 
ist also ein direktes Maf fur die Schnelle der Basilar- 
membran an derselben Stelle. Legt man, wie in 
Bild 6 eingezeichnet, in den Querzweig einen Mef- 
kondensator Cy, so lat sich an diesem eine Span- 
nung Uy = —1/Cy { AJ, dt abgreifen. Aus ihr labt 
sich gemaf der Beziehung 


(40) 


64 65 


Ry 


———o——— > 


Bild 6. Der ausgefiihrte elektrische Kettenleiter zur 
Darstellung der Bewegungsvorgange. 


4.2. Der schaltungsmaBige Aufbau der eingefiihrten 
Zweipole 


Der Kettenleiter in Bild 5 wurde zunachst aus 
allgemeinen Zweipolen aufgebaut. Soll der Ketten- 
leiter die Vorgange im Innenohr richtig wieder- 
geben, so muB8 die Frequenzabhangigkeit seiner 
Zweipole den Beziehungen (35) bis (37) geniigen. 
Setzt man in Gl. (35) bis (37) ftir zy) und zg die 
Beziehungen (4) und (24) ein, so erhalt man fir 
die Schaltelemente die gesuchten Frequenzfunktio- 
nen 


Zi=aiphe= aie Atrg+j@ Beh Azo, (41) 
Np Ny Np Ny 
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1 1 1 h bx - 1 
panes. = eee =Tp+j@——x (42 
Vane y gare nga 2 Ar bp?” Le I 
hb 1 
2 Ax bp? Lae 2 Ax bp? 
j @ ny ny -—— =. 
h bx 
| 
133 
| 
Poe ee eee h bp Tro +j @X 
Var Vege ie tieh se 48 
(43) 
1343-8 
 ( Ay 
Npny At 48 ‘ 


Der komplexe Widerstand Z, besteht demnach aus 
der Reihenschaltung eines Ohmwiderstandes 


Ry = Ax To 
Tip Thy 
und einer Induktivitat 
I= wae Axo. (44) 
Np Ny 


Da beim schaltungsmaBigen Aufbau der Mittelleiter 
weggelassen wird, kénnen, wie in Bild 6 dargestellt, 
die Leitwerte Y, von oberem und unterem Kanal zu 
einem Leitwert 3 Y, zusammengefakt werden. Die- 
ser Leitwert 4 Y, besteht nach Gl. (42) aus der Rei- 
henschaltung eines Ohmwiderstandes 


hb 
Ree 45 
= i ite At bp? > oe 
einer Induktivitat 
liebe 
Less epee 46 
qB ea he be mp (46) 
und einer Kapazitat 
CypH=Npm x by” (47) 
bx 


Der komplexe Leitwert Y,/ besteht aus der Reihen- 
schaltung eines Ohmwiderstandes 


bx 
1343 — 
pre che (13 al 
& Npn, Az 48 i 
und einer Induktivitat 
(13+35) 
Le = a en 
Hp TghNe 48 


Nach Untersuchungen von Zwistocki [4] und Baucn 
[6] ergeben sich infolge der geringen Viskositat der 
Lymphfliissigkeit fiir die Widerstinde R, undR,x so 
geringe Werte, daB sie bereits durch die Wicklungs- 
widerstande der Induktivitaten LZ; und Lx realisiert 
sind. Bild 6 zeigt den Kettenleiter nach Bild 5, wobei 


| 


1 
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hier die Zweipole nach den Beziehungen (44) bis 
(48) durch die einzelnen Schaltelemente ersetzt sind. 
Um die Auslenkung der Basilarmembran direkt mes- 
sen zu konnen, wurde die Kapazitaét C,, aufgespal- 
tet in eine MeBkapazitat Cy von geeigneter Grove 
und eine Restkapazitat. 

Der komplexe Widerstand Z, kann auf Grund der 
Untersuchungen von BExKEsy [9] durch die Reihen- 
schaltung aus dem Ohmwiderstand R, und der Kapa- 


| zitat C,, der komplexe Widerstand Z, durch die 


Kapazitat C, dargestellt werden. Dabei entsprechen 
die Kapazitaten C, und C, der Nachgiebigkeit des 
ovalen und runden Fensters. Der Widerstand R, ent- 
spricht dem Reibungswiderstand bei der Bewegung 
der Steigbiigelplatte. 

Mit Gl. (34) werden Druck und Spannung bezie- 


_hungsweise Schnelle und Strom in eine Analogie- 


beziehung gesetzt. Auf Grund dieser Wahl entspricht 
eine Induktivitét einer schwingenden Masse, eine 
Kapazitat einer elastischen Feder und ein ohmscher 
Widerstand einem Reibungswiderstand. Das bedeu- 
tet fiir die Schaltelemente in Bild 6, da8B zum Bei- 
spiel die Induktivitat L,x der senkrecht schwingen- 
den Kanalfliissigkeit entspricht. Diese Analogie- 
betrachtung darf jedoch keinesfalls dazu verfiihren, 
eben diese Induktivitat L,x proportional der in der 


_ Abschnittsbreite Ax in y-Richtung schwingenden 


Flissigkeitsmasse zu setzen, denn dann mite 


_ Lax ~ Avg sein, wahrend nach der Dimensionie- 


rungsvorschrift in Gl. (48) Lyx ~@/Az ist. Dieser 


' scheinbare Widerspruch hat seine Ursache darin, dah 


_ das Netzwerk in Bild 6 nicht die allgemeinen Diffe- 
_ rentialgleichungen fiir das dreidimensionale Modell 


auflost, sondern nur die vereinfachten Gln. (A 9). 


5. Bemessung der Schaltelemente auf Grund der 


geometrischen Abmessungen sowie aus 
physiologischen und psychologischen Messungen 


In den Analogiebeziehungen (34) zwischen den 


 mechanischen und elektrischen GroSen wurden die 


Transformationskonstanten n, und n, eingefiihrt. 


| Will man aus dem elektrischen Kettenleiter die me- 


chanischen Gréfien, zum Beispiel die Auslenkung der 
Basilarmembran, in ihrem absolut richtigen Wert 
messen, so miissen n, und n, zahlenmafig festgelegt 
werden. Ist etwa der Druck auf die Steigbiigelplatte 
des ovalen Fensters bekannt, so 1aBt sich bei bekann- 


| tem n, die Eingangsspannung Uy des Kettenleiters 


bestimmen. Diese Spannung ruft an den MeBkonden- 
satoren Cy eine Spannung Uy hervor, aus der sich 
bei bekanntem n, nach GI. (40) die Auslenkung 7 


der Basilarmembran in ihrer wirklichen Grofe er- 


_ mitteln lat. Betrachtet man die Bemessungsvor- 


schriften in Gl. (44) bis (48) fiir die elektrischen 


Schaltelemente, so sieht man, da hier nur das Pro- 
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dukt n,n, auftritt. Fiir die Bemessung des Netz- 
werkes braucht daher lediglich die Gré8e n,n, zah- 
lenmaSig festgelegt zu werden. Bei absoluten Mes- 
sungen kann dann noch eine der Transformations- 
groRen willkirlich gewahlt werden. 

Im ausgefiihrten Netzwerk wurde das Produkt 
Ny Ny so festgelegt, daB sich fir die Langsinduktivitat 
der Wert 

1G = 15 mH 


ergab. Nach Gl. (44) ergibt sich dann mit 
Ax = 0,05 cm und der Flissigkeitsdichte @ = 1 g/cm? 
fiir die Transformationskonstante 


(49) 


dyns 
Oo) 


Np Ny = 3,33 (50) 


cm 
Bild 4 zeigt den angenaherten Verlauf von h, bx 
und bg in Abhangigkeit von x. Danach ergibt sich 
gemaB Gl. (48) fiir die Querinduktivitat L,x der 
in Bild 7 aufgezeichnete Verlauf. 

Nach [4] ergibt sich fur die aquivalent schwin- 
gende Masse mp je Langeneinheit der Basilarmem- 
bran die in Bild 8 dargestellte Abhangigkeit von z. 
Mit Gl. (46) erhalt man daraus die ebenfalls in 
Bild 7 dargestellten Induktivitaten Lyp . 

Die Nachgiebigkeit der Basilarmembran konnte 


mH 
50 


40 


‘oB 


0 10 20 mm 30 


x — 


Bild 7. Die Induktivitaten im Querzweig in Abhangig- 
keit von der Ortskoordinate x. 


0 10 
x 


20 mm 30 


Bild 8. Die aquivalent mitschwingende Masse mp der 
Basilarmembran nach [4]. 


Wie 


von BkxkEsy [3] am menschlichen Gehor gemessen 
werden. Die Messung erfolgte so, daf} bei einer kon- 
stanten Druckdifferenz zwischen oberem und unte- 
rem Schneckenkanal die dabei auftretende Volumen- 
verschiebung je Langeneinheit der Basilarmembran 
ermittelt wurde. Die daraus resultierende Nachgiebig- 
keit je Langeneinheit cg ist in Bild 9 aufgetragen. 
Nach Gl]. (47) la8t sich daraus die in Bild 10 auf- 
gezeichnete GroBe der Querkapazitaten C,p berech- 
nen. 


cm2 

dyn 

10° 
5 


2 


a) a ea 7 


2-4 
10717 
5 t 
0 10 20 mm 30 
x—— 
Bild 9. Die Nachgiebigkeit cp der Basilarmembran nach 


[3]. 


20 mm 30 


0 10 
xX 


Bild 10. Die Kapazitaten Cgp im Querzweig. 


Der Ohmwiderstand R, erhielt nach [9] einen 
Wert von 9,3k®. Nach den Untersuchungen von 
BExEsy [1] schwankt die Nachgiebigkeit der Fen- 
stermembran bei verschiedenen gemessenen Gehor- 
organen etwa um den Faktor 10. Fir die beiden 
Kapazitéten C, und C, wurde ein mittlerer Wert 
von 100 nF eingebaut. 

Die bisherigen Groen des Kettenleiters, das heibt 
die Langs- und Querinduktivitaéten sowie die Quer- 
kapazitaten konnten verhaltnismafSig leicht ermittelt 
werden, da von den hierzu notwendigen geometri- 
schen Abmessungen, von der Masse mg und der 
Nachgiebigkeit cg direkte Messungen am mensch- 
lichen Geh6r vorlagen. Wesentlich schwieriger ist die 
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Bestimmung der Ohmwiderstaénde Ry im Quer- 
zweig, da von dem gema Gl. (45) dazu notwendi- 
gen Reibungswiderstand rp bis jetzt keine direkten 
Messungen bekannt sind. Andererseits liegen iiber 
das Schwingungsverhalten des gesamten Innenohres 
bei Erregung durch Sinusténe verschiedene Messun- 
gen von Bexksy [1], [2] vor. Fuhrt man dieselben 
Messungen am elektrischen Kettenleiter durch, so 
sollen die hierbei erhaltenen Mefergebnisse mit den 
direkt ermittelten moglichst gut tbereinstimmen. 
Damit. ergibt sich aus den Messungen am fertigen 
Kettenleiter riickwirkend eine Vorschrift zur Bemes- 
sung der Ohmwiderstainde. Zum Verstandnis dieses 
Vorganges sollen daher zunachst diese oben erwahn- 
ten Messungen erlautert werden. 

Wird die Membran des ovalen Fensters sinus- 
formig erregt, so erfahrt die Basilarmembran eben- 
falls eine sinusformige Auslenkung. In Bild 17 zeigt 
Kurve 1 die Amplitude dieser Schwingung in Ab- 
hangigkeit von x bei einer erregenden Tonfrequenz 
von 2 kHz. 

Demnach schwingt die Basilarmembran nicht in 
ihrer ganzen Lange mit konstanter Amplitude, viel- 
mehr hat die Schwingungsamplitude an einer be- 
stimmten Stelle ein Maximum. Dieser Ort der maxi- 
malen Erregung verschiebt sich zum Helikotrema 
hin, wenn die Frequenz des erregenden Sinustones 
erniedrigt wird. Dabei bleibt die Form der Ampli- 
tudenkurve im wesentlichen erhalten. BEKEsy [1] 
hat am menschlichen Gehor an verschiedenen Stellen 
der Basilarmembran die Frequenz ermittelt, die an 
der betreffenden Stelle die maximale Erregung her- 
vorruft. Das Ergebnis dieser Messung ist in Bild 11, 
Kurve 3 aufgezeichnet. In Bild 11 zeigt Kurve 2 die 
Ergebnisse, die aus psychologischen Messungen er- 
mittelt wurden [10]. Die Ubereinstimmung der Er- 


Bild 11. Die Frequenz des erregenden Sinustones in 
Abhangigkeit vom Ort x auf der Basilarmem- 
bran, an dem die betreffende Frequenz die 
maximale Auslenkung hervorruft. Kurye | 
zeigt den am elektrischen Kettenleiter gemesse- 
nen Verlauf, Kurve 2 den Verlauf nach psycho- 
logischen Messungen [10], Kurve 3 den Ver- 
lauf nach direkten Messungen am Gehérorgan 


[1]. 
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gebnisse derselben Messung am elektrischen Ketten- 
leiter mit den beiden oben genannten Mefergebnis- 
sen dient uns als eines von drei Kriterien zur Be- 
stimmung der Ohmwiderstande im Querzweig. 

Als zweites Kriterium dient folgender Versuch, 
den BEKEsy am menschlichen Gehor ausgefthrt hat 
[2]. Hierbei wurde die Basilarmembran wiederum 
sinusformig erregt. Nach beendetem Einschwingvor- 


- gang wurde die Erregung am ovalen Fenster sprung- 


wna 


-haft abgeschaltet. Der entstehende Ausschwingvor- 


gang wurde fir verschiedene Tonfrequenzen jeweils 
an der Stelle maximaler Auslenkung der Basilar- 
membran beobachtet. BEKEsy fand, daf die Basilar- 
membran dort etwa mit der Frequenz ausschwang, 
mit der sie vorher erregt worden war. Bild 12 zeigt 
als Beispiel einen solchen Ausschwingvorgang. Fir 
das Verhaltnis A)/A, zweier aufeinander folgender 
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Bild 12. Das Ausschwingverhalten der Basilarmembran 
nach dem Abschalten einer sinusformigen Er- 
regung. Die Messung wurde fiir verschiedene 
Tonfrequenzen durchgefthrt, wobei jeweils an 
der Stelle der maximalen Erregung gemessen 
wurde. Im unteren Diagramm zeigen die ein- 
gezeichneten Punkte Messungen am elektri- 
schen Kettenleiter, das schraffierte Rechteck 
den yon G. vy. BExEsy [2] ermittelten Bereich. 


Amplituden erhielt BEKEsy an allen von ihm unter- 
suchten Stellen als Streubereich aus verschiedenen 
Messungen den Wert 4 bis 6. In Bild 12 zeigt das 
schraffierte Rechteck diesen Bereich. Es ist offen- 
sichtlich, da8 der Reibungswiderstand rpg dieses Am- 
plitudenverhaltnis direkt beeinflu8t. Dementspre- 
chend wird dieses Kriterium bei der Bemessung des 
Ohmwiderstandes Ry, besondere Bedeutung haben. 

Als drittes Kriterium dient ein weiterer Versuch 


4 von BExKEsy [2], bei dem das Verhdltnis der Steig- 


biigelschwingungsamplitude 7,4 zur Schwingungs- 
amplitude 73 der Basilarmembran an der Stelle der 
maximalen Erregung gemessen wurde. Das Ergeb- 
nis ist in Bild 13, Kurve 2, aufgezeichnet. 

Zur Bestimmung der Ohmwiderstinde Ray wurde 
in einem langwierigen Verfahren der Verlauf der 
Widerstandswerte in Abhangigkeit von x systema- 
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Bild 13. Das Verhaltnis der Steigbiigelschwingungs- 
amplitude zur Schwingungsamplitude der 


Basilarmembran bei sinusférmiger Erregung, 
gemessen am Ort der maximalen Auslenkung. 
Kurve 1 zeigt die MeBergebnisse am elektri- 
schen Kettenleiter, Kurve 2 die direkten Mes- 
sungen nach [2]. 


tisch so lange variiert, bis die oben angeftihrten Kri- 
terien moglichst gut befriedigt waren. Bei diesem 
Verfahren zeigte es sich, da bei einer bestimmten 
Anderung von R,, der Einflu8 auf die MeSergeb- 
nisse bei Kriterium 2 viel groBer war, als bei den 
beiden tibrigen Kriterien. Es ist leicht einzusehen, 
da8 zur Erfillung von Kriterium 1 in erster Linie 
die Induktivitaten und die Kapazitéten mafgebend 
sind. Weiterhin zeigte es sich, da} alle drei Krite- 
rien nicht gleichzeitig erftllt werden konnen. Bei der 
endgiltigen Dimensionierung der Ohmwiderstande 
gema}} Bild 14 wurde in erster Linie auf die Erfil- 
lung von Kriterium 2 geachtet. In den Bildern 11 
und 13 zeigt jeweils Kurve 1 den am endgiiltig aus- 
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Bild 14. Die Widerstande Rgp im Querzweig in Abhan- 
gigkeit von der Ortskoordinate x. 


gefiihrten elektrischen Kettenleiter gemessenen Ver- 
lauf. In Bild 12 ist der ebenfalls am elektrischen 
Kettenleiter gemessene Verlauf des Amplitudenver- 
haltnisses zweier aufeinander folgender Amplituden 
beim Ausschwingen aufgezeichnet. Es sei noch er- 
wahnt, dai die Messungen am elektrischen Ketten- 
leiter genau so wie die Messungen von BEkEsy er- 
gaben, daf} die Basilarmembran an der Stelle des 
Maximums etwa mit der Frequenz ausschwang, mit 
der sie erregt worden war. Es ist einzusehen, dai 
dieses Verhalten nicht von den Widerstanden Rg, 
sondern in erster Linie von den Spulen und Kapa- 
zitaten im Querzweig beeinfluft wird. Die Tatsache, 
dai der elektrische Kettenleiter dieses spezielle Ver- 
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Bild 15 a,b,c. Die zeitliche Auslenkung der Basilar- 
membran bei Erregung am ovalen Fenster 
durch Druckimpulse gleicher Amplitude und 
verschiedener Dauer t;. Der Verlauf des er- 
regenden Druckimpulses ist jeweils in der 
rechten oberen Ecke aufgezeichnet. Parameter 
der Kurven ist der Ort x, an dem der Bewe- 
gungsablauf gemessen wurde. 
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halten der Basilarmembran richtig wiedergibt, hat 
seine Ursache in der Beriicksichtigung der in y-Rich- 
tung schwingenden Kanalfliissigkeit. Diese Beriick- 
sichtigung fiihrte zu den Induktivitaten Lax, die 
ihrer GroBe nach die Gesamtinduktivitat im Quer- 
zweig wesentlich bestimmen. 

Nach denselben Gesichtspunkten wie bei der Be- 
stimmung von R,p wurde auch der Widerstand 
Ry =40 Q bestimmt. Er entspricht dem Reibungs- 
widerstand der Kanalfliissigkeit beim Durchfluf 
durch das Helikotrema. 

Es ist an dieser Stelle angebracht, ein Wort uber 
die erzielte Ubereinstimmung zwischen den direkten 
Messungen und denen am elektrischen Kettenleiter 
zu sagen. Auf der einen Seite wurden bei der Dimen- 
sionierung des elektrischen Kettenleiters sehr ver- 
einfachte Modellvorstellungen zugrunde gelegt, und 
im Verlauf der Rechnung weitere einschneidende 
Vereinfachungen eingefiihrt. Andererseits ist wegen 
der Streuung zwischen den einzelnen Gehororganen 
die Genauigkeit der direkten Messungen begrenzt. 
Unter Beriicksichtigung dieser Umstiande ist die er- 
zielte Ubereinstimmung gut. 


6. Einige MeBergebnisse tiber das Schwingungs- 
verhalten der Basilarmembran bei Erregung mit 
Druckimpulsen und Sinusténen 


Das Verhalten des Gehors bei impulsformiger 
Erregung ist sowohl bei psychologischen als auch bei 
physiologischen Messungen von besonderem Inter- 
esse. In diesem Abschnitt wird die Auslenkung der 
Basilarmembran bei Erregung durch sogenannte 
Druckimpulse untersucht. Dabei wurde am Eingang 
des Kettenleiters ein Spannungsimpuls mit dem zeit- 
lichen Verlauf nach einer Gaufschen Fehlerkurve 


—/2 (t/t,)® (51) 


Up =tige 


eingepragt. Die Bilder 15 a, b, c zeigen jeweils in 
der rechten oberen Ecke den Verlauf eines solchen 
Impulses. Die in Gl. (51) eingefiihrte Impulsdauer 
t; wurde so definiert, da ein Rechteckimpuls mit der 
Amplitude i) und der Breite ¢, dieselbe Energie wie 
der zugehorige Druckimpuls nach Gl. (51) beinhal- 
tet. Dieser Rechteckimpuls gleicher Energie ist in 
Bild 15 © gestrichelt eingezeichnet. Die Bilder 15 a, 
b, c zeigen den zeitlichen Verlauf der Auslenkung 
der Basilarmembran an verschiedenen Stellen bei 
Erregung durch solche Druckimpulse gleicher Am- 
plitude und verschiedener Impulsdauer 4;=55 us, 
140 us, 570 us. Dabei ist die Auslenkung der Basi- 
larmembran auf einen willkiirlichen Wert bezogen, 
der jedoch fiir simtliche drei Abbildungen derselbe 
ist. Parameter der Kurven ist die Stelle 2 auf der 
Basilarmembran, an der die Auslenkung gemessen 
wurde. 
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Betrachtet man den zeitlichen Verlauf der Aus- 
lenkung an den verschiedenen Stellen, so ist zunachst 
zu bemerken, dafi dieser schon in der Nahe des ova- 
len Fensters nicht mehr dem erregenden Druckimpuls 
gleicht. Die Ursache dafiir ist in erster Linie die Stei- 
figkeit der beiden Fenstermembranen. Im elektri- 
schen Kettenleiter sind sie durch die im Langszweig 
liegenden Kapazitaten C, und C, dargestellt, die fiir 
tiefe Frequenzanteile eine Differentiation bewirken. 
Die Folge davon ist das Uberschwingen der Basilar- 
membran. 

Die Bilder 15 zeigen, da die Auslenkung deutlich 
am ovalen Fenster beginnt und von da aus in Form 
einer Wanderwelle tber die Basilarmembran hinweg- 
lauft. Ein markanter Punkt fir das Vorriicken der 
Wellenfront ist der Augenblick, in dem die Basilar- 
membran zum ersten Mal von oben nach unten durch 
ihre Nullage schwingt. Ein wichtiges Ergebnis dieser 
Untersuchungen ist, daB die Geschwindigkeit, mit 
dem sich dieser Punkt der Wellenfront zum Heliko- 
trema hin bewegt, von der Impulsdauer ¢; unabhan- 
gig ist. 

In Bild 16 zeigen die eingezeichneten Punkte den 
Mittelwert der so definierten Wellengeschwindigkeit 
vw aus den Messungen fur die drei untersuchten 
Impulsbreiten. Der Verlauf dieser Punkte ]a8t sich 
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Bild 16. Wellengeschwindigkeit vy der Wanderwellen 
auf der Basilarmembran in Abhangigkeit von 
der Ortskoordinate x. Die eingezeichneten 
Punkte zeigen die aus Messungen am elektri- 
schen Kettenleiter ermittelten Werte. Die ein- 
gezeichneten Kreuze wurden aus direkten Mes- 
sungen nach [2] ermittelt. 


gut durch die in Bild 16 eingezeichnete Gerade an- 
nahern. Da der Ordinatenmafstab logarithmisch ge- 
wahlt wurde, bedeutet dies, da} die Wellengeschwin- 
digkeit vw zum Helikotrema hin exponentiell nach 
folgender Beziehung abnimmt: 


vw (x) =vw(0) e 7%, (52) 
Dabei ist 
vy (0) =55 m/s, (53) 
%y=10mm. 


Selbst die Anfangsgeschwindigkeit von 55 m/s ist 
noch viel kleiner als die Schallgeschwindigkeit im 
Wasser von 1450 m/s. 
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Einen ahnlichen Versuch hat BEKEsy am mensch- 
lichen Gehor direkt vorgenommen [2], indem er das 
ovale Fenster durch einen Funkenknall erregte. In 
Bild 16 zeigen die eingezeichneten Kreuze die aus 
dieser Messung ermittelte Wellengeschwindigkeit fiir 
den Funkenknall. In dem kleinen Bereich, fiir den 
die direkten Messungen vorliegen, ist die Uberein- 
stimmung mit den Messungen am elektrischen Ket- 
lenleiter gut. 

Im Verlauf ihrer Ausbreitung auf der Basilar- 
membran verandert die Wanderwelle ihre Form in 
zweifacher Hinsicht. Einmal nimmt die Dauer der 
Auslenkung um so mehr zu, je weiter die Wander- 
welle zum Helikotrema hin vorriickt. Bei einer etwas 
groben Betrachtungsweise laBt sich dieser Vorgang 
so deuten, da der Impuls eine Tiefpabkette gemaB 
Bild 6 durchlauft, deren Grenzfrequenz von links 
nach rechts schrittweise abnimmt, so daB die hoheren 
Frequenzanteile des Impulsspektrums immer mehr 
unterdruckt werden. Bemerkenswert ist noch, da in 
der Nahe des Helikotremas die Dauer der Auslen- 
kung von der Dauer ¢; der erregenden kurzen. Druck- 
impulse nahezu unabhangig ist. Weiterhin verandert 
sich die Amplitude der Wanderwelle im Verlauf 
ihrer Ausbreitung. Sie hat an einer bestimmten 
Stelle ein Maximum, das um so naher beim Heliko- 
trema liegt, je groBer die Dauer des erregenden 
Druckimpulses ist. 

Da das Maximum in Abhangigkeit von x sehr flach 
verlauft, kann der Ort des Maximums nur sehr un- 
genau bestimmt werden. Zu erwahnen ist noch, dah 
die Basilarmembran an den Stellen, die rechts von 
dem eben definierten Amplitudenmaximum liegen, 
periodisch ausschwingt, wahrend sie an Stellen links 
davon nur einmal tberschwingt. Die Frequenz des 
periodischen Ausschwingens entspricht etwa der Fre- 
quenz des Sinustones, der gemaB Bild 11 an dieser 
Stelle seine maximale Erregung hat. 

Bild 17 zeigt die Amplitude der Auslenkung der 
Basilarmembran bei Erregung durch einen Sinuston 
der Frequenz 2 kHz in Abhangigkeit von der Orts- 
koordinate x. Die Kurve 1 zeigt den Verlauf der 
Amplitude bei normaler Einspeisung am ovalen Fen- 
ster gemas Bild 6. Die Kurven 2 bis 5 zeigen den 
Verlauf der Amplitude, wenn die Einspeisung zwi- 
schen verschiedenen Punkten x, im Langszweig von 
Kanal 1 und dem Bezugspunkt M erfolgt. Der Punkt 
der Einspeisung war bei Kurve 2: 2,=12,5 mm; 
bei Kurve 3: 2,=17,5mm; bei Kurve4: 2,= 
22.5mm; bei Kurve 5: x, =28,5mm. Der Anlaf 
zu solchen Messungen war ein bestimmter operativer 
Eingriff, die sogenannte ,,Fensterung“, der in der 
modernen Chirurgie an erkrankten menschlichen 
Gehdrorganen schon Ofters vorgenommen wurde. 
Diese Operation ist dann notwendig, wenn infolge 
einer Versteifung der Membran des ovalen Fensters 
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Bild 17. Amplitude der Auslenkung der Basilarmem- 
bran in Abhangigkeit von x bei Erregung durch 
einen Sinuston der Frequenz f=2 kHz. 
Kurve 1 zeigt den Verlauf der Amplitude bei 
normaler Einspeisung am ovalen Fenster. Die 
Kurven 2 bis 5 zeigen den Verlauf der Ampli- 
tude bei Einspeisung zwischen verschiedenen 
Punkten x, im Lingszweig von Kanall und 
dem Bezugspunkt M. 

Kurve 2: x,=12,5 mm, 
Kurve 3: x,=17,5 mm, 
Kurve 4: x,=22,5 mm, 
Kurve 5: x,=28,5 mm. 


das Innenohr durch den ankommenden Schall nicht 
mehr in Schwingungen versetzt werden kann. Das 
Ziel der Operation ist, fiir den erregenden Schall 
einen anderen Zugang zum Innenohr zu schaffen. 
Hierbei wird der Schneckenkanal an irgend einer 
Stelle durch das Felsenbein hindurch angebohrt. Die 
diinne Haut, die das Innenohr zur Knochenmasse hin 
abschlieBt, darf bei dieser Bohrung nicht verletzt 
werden. Sie verhindert das Auslaufen der Peri- 
lymphe und dient als neue Fenstermembran, tber 
die der erregende Schall die Basilarmembran in 
Schwingungen versetzen kann. Erstaunlicherweise 
ergab sich hierbei nicht, wie vielleicht zu erwarten 
gewesen ware, eine Verschiebung der Tonhdhen- 
empfindung des Patienten. Lediglich die empfundene 
Lautstarke gegeniiber dem gesunden Ohr war gerin- 
ger. Diesen operativen Eingriff naherungsweise 
nachzubilden, war das Ziel der Untersuchungen von 
Bild 17. Dabei ergab sich, daB der Ort der maxima- 
len Erregung der Basilarmembran vom Ort der Ein- 
speisung vollig unabhangig ist. Lediglich die GroBe 
des Maximums nahm um so mehr ab, je naher die 
Einspeisung am Helikotrema erfolgte. Da der Ort der 
maximalen Erregung ein Mafi fiir die empfundene 
Tonhohe ist, wird verstandlich, daB die Tonhohen- 
empfindung durch den operativen Eingriff nicht ver- 
andert wird. 

Wir danken Herrn Prof. Dr. rer. nat. Dr.-Ing. 
E.h. R. Fecprkevter und Herrn Dr.-Ing. E. Zwicker 
fiir viele wertvolle Ratschlage. Unser Dank gilt der 
DFG fiir die Unterstiitzung der Arbeit, ferner der 
Firma Siemens & Halske, die uns freundlicherweise 
die ae fiir den Aufbau des Kettenleiters 


zur Verfiigu®® stellte. 
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Anhang 


Die Differentialgleichungén der Bewegungsvorgange 
im dreidimensionalen Modell. 

Die Kontinuitatsgleichung und die Bewegungsglei- 
chungen fiir inkompressible Fliissigkeiten haben im 
dreidimensionalen Raum fiir sinusformige Vorgange die 
Form 


divw=0, (A 1) 
Sp (9,2) + 2% Vx (x, y, z) =0, 
Ox 
) 
SPU) + sy oy(c,y,2)=0,b (A) 
gz. 
Op eS) Zp Vz (2, y, Z) =0. 
Duar 


Mittelt man die Kontinuitaétsgleichung (A 1) iiber den 
Querschnitt Q,=h by des Kanals 1 im Modell von 
Bild 4, so erhalt man unter Beriicksichtigung der Rand- 
bedingungen 


i ¢ a Ovz (a) Raley 
kig_| J avvat= a 5 be eee (A 3) 
a.) 
Dabei ist 
hae | Hi vs (ay. 3) df Mental 
h bx f 
(Q1) 


die tiber den Kanalquerschnitt Q, gemittelte Geschwin- 
digkeitskomponente in z-Richtung. Gl. (A 3) entspricht 
der Kontinuitatsgleichung (8), in die sie fiir by =bx 
tibergeht. Fiihrt man dieselbe Mittelung fiir die erste 
der Gleichungen (A 2) durch, so erhalt man unter 
Beriicksichtigung von Gl. (A 4) 


c [ fee me) af +29 %2(2)=0. (AS) 
Qi 


3 
A 
Ox 


Die Reihenfolge von Integration und Differentiation 
konnte wie in Abschnitt3. vertauscht werden, da 


w< 1 und we <1 ist. Das erste Glied enthalt den 
x x 


mittleren Druck im Kanalquerschnitt Q,. Wie bei der 
Rechnung fiir das zweidimensionale Modell wird er 
durch den Druck unmittelbar auf der Basilarmembran 
ersetzt. Aus der zweiten der Gleichungen (A 2) erhalt 
man durch Integration iiber y 

y 

p(x, y,z) =p(x, 0,2) —29 [ vy(x, ys 2) dy. 
ty 


Durch Mittelung iiber Q, erhalt man daraus 


bx 
1 1 1 
i [ [rena a= [eee 02) ds t95 5X 
(Q,) 0 


«ff [foverns aya. 6) 


(Gy 


Wie in der bereits erwahnten Arbeit von H. Hauser 
gezeigt wird, gilt folgende Naherung: 
bx bz 


£3 foe. 0, z) dewet | p(z, 0, z) dz = p(z); 
bx bz 
6 0 (A 7) 


‘ 
eo 
: 
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p(x) ist der gesuchte Druck auf der Basilarmembran, 
gemittelt iiber ihre Breite bp. 

Fiir die Umformung des zweiten Gliedes von Gl. 
(A 6) wird wie im Abschnitt 3.1. die Annahme getrof- 
fen, daB die Fliissigkeit in der Flache Q, gleichphasig 
mit der Basilarmembran schwingt.’ Entsprechend Gl. 
(12) wird der Ansatz eingefihrt: 


vy (x, ¥y, z) = —7 (x) Fly, z) . 


Damit erhalt man fiir das zweite Glied auf der rechten 


Seite von Gl. (A 6) 


2 og [ [|] ro ay] af. 


Nimmt man fiir F(y,z) die stiickweise definierte Funk- 
tion 


O25 b,, 
y bxe-z 
F(y, z)-= > = ; 
(y, 2) eae 
h 
fiir 0S7S 5: bp<zS bx, 
UN (Diez 
F if 
(yz) ( i) Pome. 
fir Sy Sh, bg <zS br, 


so erhalt man fiir den obigen Ausdruck 

; ch bp 

——{13—— +3}. 

Zo (x) 48 ( oe 

Damit ist der tiber Q, gemittelte Druck in Gl. (A 6) 

durch die gewiinschten Groen p(x) und 7 (x) ersetzt 

und aus Gl. (A 5) erhalt man die Bewegungsgleichung 
in der gewiinschten Form 

so ta(2) + © |p(e) +20 (13 72 +3) 1G) 

dx 48 bx 


0% 

(A 8) 
Gl. (A 8) entspricht Gl, (14), in die sie fiir bg =bx 
ubergeht. Der weitere Rechnungsgang ftir die Bewe- 
gungsgleichung und die Kontinuitatsgleichung verlauft 


vollkommen gleich wie in Abschnitt 3. In den Bewegungs- 
gleichungen steht an Stelle von h/3 der Ausdruck h/48 


(13 bp/bx +3), in den Kontinuitatsgleichungen ist je- 


weils das zweite Glied mit dem Faktor bp/bx multipli- 
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ziert. Das Ergebnis ]a$t sich also direkt hinschreiben: 


d b 2 
Pg . he (p1,2—pu) =0, 
fiir Kanal 1 : e 
beziehungs- as ee ee 
NES: Zo viet [(1+ 6) (p12—pH)]=0. 
(A 9) 
Dabei ist 
1 bp Zyh 1 ( bp 
== 13 +3 13 L3)a.(A10 
B 24 ( bx ZB 16 bx ae ) 
(Eingegangen am 5. Januar 1961.) 
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Summary 


The Fr6uticu-Sack method is applied to study the compressibility of two phase mixtures 
and an unified explanation of the behaviour of various types of mixtures (sintered materials, 
foams, emulsions, dispersions etc.) is given. The expressions used by Wooo, MacKenziz, 
Kepeue and others are deduced under suitable approximations. 


Sommaire 

La méthode Fréuurtcu-Sack est appliquée a l’étude de la compressibilité des mélanges 
sous deux phases et donne une explication du comportement de différentes sortes de mélan- 
ges (matériaux granulés ou agglomérés, mousses, émulsions, milieux dispersifs, etc.). On en 
déduit par des approximations convenables les formules utilisées par Woop, MacKenzin, 
Kerner et d’autres auteurs. 


Zusammenfassung 


Mit Hilfe der Froénuicu-Sacx-Methode wird die Kompressibilitat von Zwei-Phasen- 
Mischungen untersucht und das Verhalten verschiedener Mischungstypen (gesinterte Mate- 
rialien, Schaume, Emulsionen, Dispersionen etc.) erklart. Die von Woop, MacKenzie, 


Kerner und anderen benutzten Formeln werden mit geeigneten Naherungen abgeleitet. 


1. Introduction 


The propagation of sound in two phase mixtures 
like emulsions, dispersions, polycrystalline aggre- 
gates, sintered materials, foam etc., are of impor- 
tance for practical applications. Recently there have 
been several discussions of the velocity of sound in 
various types of two phase mixtures (Uricx [1], 
MacKenzie [2], Caampre [3], Kerner [4]). Woop 
[5] introduced a very simple assumption of the 
additivity of the volume compression, that is, if 1; 
and vp are the volumes of the inner and the outer 
phases, /; and fy are their compressibilities, the 
compressibility 6, of the total equivalent medium 
is given by 

Be=0 Bi+ (1-9) Bos O=vj/ (vj + v9) (1) 


and this expression agrees well with several obser- 
vations (McKinney, Eprtman and Marvin [6], 
Auuinson [7]). The “elementary“ derivation of this 
result (for example, CuamBre [3]) on the basis, 
that the total volume change under the applied pres- 
sure is the sum of the volume changes of the indivi- 
dual phases under the same pressure, is somewhat 
fallacious since the pressure is an extensive property 
having the same value at all points only in a fluid. 
For example, the “elementary” derivation of eq. 1 
leads to the result that a solid with a vacuous cavity 
has infinite compressibility (see however eq. 5) 
just like a liquid with a cavity. 


In this note, the Fr6nticu-Sacx self-consistent 


approximation (FRéuLicHh and Sack [8]), which 
has been successful in explaining the viscoelastic 
properties of dispersions (OLtpRoyp [9], RazacopaL 
[10]), is applied to calculate the compressibility of 
the general two phase system. The various expres- 
sions that have been given by other workers can 
then be discussed as corollaries. 


2. Frohlich-Sack model for the two phase system 


In the FréuticH-Sack method, one confines the at- 
tention to a region close to the particle of the dipersed 
phase. Though the immediate environment of the par- 
ticle of the inner phase is the material of the outer 
phase, the medium at some distance from the particle 
is essentially the smoothed-out equivalent medium. 
This is true to a high degree if the concentration 0 
of the inner phase is small. So one considers a 
sphere of radius r=a of the inner phase (bulk mo- 
dulus k; = 4; + 2 4;/3 , shear modulus 4;) surrounded 
by the outer phase (ky, 49) to a radius r=rg such 
that a/ro>=0. This is surrounded to a large 
distance R, (R>ry), by the equivalent medium 
(hess): 

This composite sphere differs from a fully con- 
tinuous sphere of the equivalent medium of the same 
radius R only in the portions at r < ry and the lead- 
ing term in the contribution of this part to the bulk 
behaviour under a hydrostatic pressure at r= R may 
be expected to be of the order of the ratio of the 
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two volumes 7 9?/R*. So, to get the equivalent me- 
dium, the Frouuicu-Sack postulate is that the bulk 
displacements of the two spheres under a hydrosta- 
tic pressure at r=R must be identical to the terms 
including (r)/R)*. ; 

Using Lamr’s well known results (Love [11]), 
the problem is easily solved. With spherical sym- 
metry, the general equation of equilibrium with 
zero body force is 


0 /Ou 2u\ — 
Oe ee “)=0 


with the arbitrary solution u=Ar+Br~*. One also 
has the radial stress 


rr=(A+2q) Se 
T 


In the composite sphere, one assumes the general 
solutions u;=A;r+B,r-2 in the inner phase, 

‘Ug=Agr+Byr * in the outer phase and 
Ue=A.r+B,r ? in the equivalent medium. Within 
the homogeneous sphere of the equivalent medium 
the solution will be wu,’ = A,’ r’ and B,” will be zero, 
as otherwise the displacement will be infinite at the 
origin r=0. 

The A, B constants are to be evaluated from the 
conditions of the finiteness of the displacement at 
r=0, the continuity of the displacements and pres- 
sures at the interfaces r =a and r =rg and the FR6u- 

uicn-Sack postulate at r=R. The first condition 
gives B;=0 while the last yields A, =A, and 
B.=0. The continuity of the displacements and 

pressures at r =a leads to 

ba 


A,a=Aja+Bya’; 
(3A; +2 Mj) Ai= (3 4p +2 Mo) Ap—4 mo a? By 


; (2 a) 
Ay = Aj(3 ki +4 fg) /(3 ky +4 Mo) 5 
“or By = A; (3 ko — 3 ke) | (3 ko + 4 Mo) - 
The similar conditions at r = ry results in 
Ay= Ae (3 lke +4 ftp) | (3 ky +4 Ho) § 
By = Agta (oe ky — 3h) (3 ko t+ 4 bo) - (2b) 


On eliminating A)/By from eqs. (2a) and (2b). 
pas ky ki +3 ky Mo (1 —9) ote $ kj tg 


: (3) 
ki(1—9) +k A+ % to 
The equivalent density is given by 
0. =9.91 + (1—9) a (4) 


and the sound velocity is obtained from eqs. (3) 


and (4). 
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3. Discussion 


The general expression eq. (3) is identical with 
the result obtained for two phase systems by Kerner 
[4]. Kerner claims his method to be an improve- 
ment over the Frouticu-Sack technique, but it ap- 
pears from our result that this is not so. 

When the internal phase is gaseous, k; may be 
taken as zero compared to ky and eq. (3) then re- 
duces to 


ke *= hg -*(h— 0) 8+ 3 ug 1 (1 —0) (3) 


which is the expression of Mackenzir [2] for sinter- 
ed materials, foams etc. 

When the outer phase is a liquid, “) tends to zero 
and one gets eq. (3) 


hy! Sh 6+ 118) (6) 


which agrees with the assumptions of Woop [5] and 
is valid for emulsions and solid-in-liquid dispersions. 
As mentioned earlier, the “elementary” deviations of 
Woop’s assumption are not free from objections 
regarding the range of their validity. From our 
results, it is not surprising that, when the outer 
phase is a fluid, Woop’s equation represents the 
experimental observations quite well, as noticed by 
Urick [1], Aturnson [7] and others. Other effects 
come into play only when the sound wavelength 
approaches the particle size (Urick and Ament [12]), 
or when the sound frequency interacts with the re- 
sonant frequency of the particle (for example S1z- 
BERMAN [13]), or when the visoelastic behaviour of 
the dispersion medium becomes prominent, none of 
which is considered in the theory outlined here. 

The rigidity of two phase mixtures can be analys- 
ed similarly using the results obtained in another 
connection by Gooprrr [14] and the extension of 
the present discussion to multiphase systems is also 
straightforward. 

The author thanks Professor R. S. Krisunan for 
his guidance and encouragement and the Council 
of Scientific and Industrial Research for the award 
of a Senior Research Fellowship. 

(Received 21th January 1961.) 
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Ein momentan-logarithmischer Verstarker fiir das Tonfrequenzgebiet 


Viele Halbleiterdioden weisen tiber mehrere Dekaden 
des Stroms eine exponentielle Kennlinie auf und schei- 
nen deshalb zum Bau eines logarithmischen Abschwa- 
chers gut geeignet [1]: gibt man den Diodenstrom 
durch eine Stromquelle vor, so hangt die Spannung an 
der Diode logarithmisch vom angelegten Strom ab. Bei 
der technischen Realisierung einer solchen Anordnung 
treten jedoch mehrere Probleme auf. Selbst Dioden, 
deren Kennlinie weitgehend der theoretisch zu erwar- 
tenden entspricht, weisen keine rein exponentielle Ab- 
hangigkeit des Stroms von der angelegten Spannung 
auf. Hinzu kommt im allgemeinen die Stérung durch 
einen mehr oder weniger gro8en Bahnwiderstand, der 
mit der Diode in Reihe geschaltet erscheint. SchlieBlich 
bereitet die erforderliche Stromsteuerung der Dioden 
einige Schwierigkeiten. Im folgenden wird gezeigt, wie 
man mit verhaltnismaig einfachen Mitteln diese Pro- 
bleme lésen kann. 


1. Das Prinzip 


Der logarithmische Verstarker soll beide Stromrich- 
tungen verarbeiten (speziell bei Sinusstrémen) ; es sind 
daher zwei gleiche Dioden antiparallel zusammen- 
zuschalten !. Die Kennlinien der Dioden ohne Beriick- 
sichtigung des Bahnwiderstandes haben die Form [1] 
(Bild 1 a): 


Diode 1: er 


L —eulur_] . 
‘ I, a 
Diode 2: — 7h =e “Ur_], 
T 
Summe: —= T+, 


Is Is 


I, Sattigungsstrom, 


= etl/ur _ e Ulur : 


uy Temperaturspannung, 
up=k T/e=26 mV bei Zimmertemperatur. 


fe 


Bild 1 a. Diodenpaar zur Logarithmierung. 


Der Stérterm —1 hat sich bei der Parallelschaltung 
kompensiert; die gewiinschte exponentielle Charakteri- 
stik fiir die eine Stromrichtung wird jetzt jedoch durch 
1 Nahere Ausfiihrungen dariiber, was man hier unter 
dem Logarithmus einer Sinusfunktion verstehen muf, 
bei [2]. 


den zweiten Term fiir die entgegengesetzte Stromrich- 
tung gestért. Die Storung nach Auflésung der obigen 
Gleichung ist 


otter 5 [1+ VI ta 
det enti toe ts ee 
is I 


Fiir einen Fehler von weniger als 1 dB & 12% in J er- 
gibt sich die Bedingung 
(PORN irs. 


Bei den in dem beschriebenen Gerat verwendeten Gold- 
drahtdioden OA180 der Firma Telefunken ist 
I;~1uA, der maximale Strom betragt etwa 1 mA; 
daraus errechnet sich dann eine Dynamik von 51 dB, die 
wohl fiir die meisten Zwecke zu gering ist. Eine wesent- 
siche Verbesserung wird dadurch erreicht, da} die bei- 
den Dioden jeweils in Sperrichtung vorgespannt werden 
(Bild 1 b). Nun wird der Gesamtstrom 


wis = e(U—%) Uy _ @— (UU) /Uey 

s 

beziehungsweise 
I 


= ellur fia ee U/ur 

Is : 
Der wirksame Sittigungsstrom /;’ erscheint jetzt um 
den Faktor e~“/“t yerringert und die Dynamik entspre- 
chend erhoht; fiir 30 dB Dynamikzunahme zum Beispiel 
wird ein uw) =90 mV bendtigt. 


mit J,’ =I, e~ “Mr. 


Bild 1 b. Dioden mit Vorspannung. 


Die so vorgespannten Dioden sind bei den genann- 
ten Werten fiir Anfangs- und Sattigungsstrom in der 
Lage, etwa 80 dB der Eingangsspannung zu verarbei- 
ten, wenn es gelingt, in diesem ganzen Bereich den 
Eingangsstrom yorzugeben. Im allgemeinen wird man 
einen Spannungsverstarker mit hohem Innenwiderstand 
als Generator verwenden, dessen Urspannung immer 
groB ist gegeniiber der an den Dioden liegenden Span- 
nung. Eine Abschatzung zeigt, daB8 dann ein Verstar- 
ker mit der maximalen Leerlaufspannung 3kV und 
dem Innenwiderstand 3 MOhm nétig ist, der natiirlich 
nicht gut zu realisieren ist. 

In dem ausgefiihrten Geraét wird diese Bedingung 
daher auf eine andere Weise erfiillt. Die Dioden wer- 


im 
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den zwischen Eingang und Ausgang eines einfachen 
Rechenverstarkers geschaltet; die Eingangsspannung 
am Verstarker ist dann um einen Faktor 1/u (—u 
Leerlaufverstarkung) kleiner als die an den Dioden 
liegende Spannung, die praktisch gleich der Ausgangs- 
spannung ist. Der steuernde Verstarker braucht dann 
nur eine Leerlaufspannung und einen Innenwiderstand 
aufzuweisen, die um den gleichen Faktor kleiner sind. 
Die zur Bemessung wichtigen Gleichungen werden an 
Hand der Skizze (Bild1c) hergeleitet. Der Ausgang 
laufe leer, dann ist 


Ue +up = — MM Ue— Rj ig =Uag 
1] F : 
vem l+u [up (ie) + Rite] 
, IR , 
Ug = he | (ie) ey = up (ie). 


Additiy zu up ist ein Term Rg ie hinzuzufiigen, wenn 
wir den Bahnwiderstand Rg der Dioden beriicksichtigen 
wollen. Durch geeignete Wahl der Verstéarkung uu und 
des Innenwiderstandes Rj ist es jedoch méglich, diesen 
storenden Einflu8 zu kompensieren, wie man leicht an 
der letzten Gleichung abliest. 


_ Bild 1c. Prinzipschaltbild des Rechenverstarkers 


(—w Leerlaufverstarkung). 


2. Ausfiihrung und Eigenschaften des Gerates 


Fiir den zu verwendenden Verstarker mu man ein 
geniigend hohes ~ (10?) und einen nicht zu groBen 
Innenwiderstand (etwa 1kOhm) fordern. Auferdem 
mu er absolut stabil bleiben, selbst wenn Eingang und 
Ausgang miteinander verbunden werden. Diese Bedin- 
gungen erfillt in einfachster Weise ein Pentodenver- 
starker mit galvanisch angekoppeltem Katodenfolger. 

Als Dioden wahlt man Typen mit kleinem Bahn- 


3 0A180 
Bild 2. Schaltbild des logarithmischen Verstarkers. 
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widerstand und Sattigungsstrom. Giinstig sind an und 
fiir sich Siliziumkleinflachendioden, weil sie Sattigungs- 
stroéme in der Grofenordnung 10~° bis 1072 uA auf- 
weisen, so daf} sich im allgemeinen eine zusatzliche Vor- 
spannung eribrigt. Leider sind sie fiir unsere Zwecke 
nur bis zu Frequenzen von einigen kHz zu gebrauchen; 
dartiber hinaus nimmt die erzielbare Dynamik sehr 
schnell ab. Wesentlich besseres Verhalten zeigen Gold- 
drahtdioden, die ebenfalls kleine Bahnwiderstande 
haben, jedoch einen héheren Sattigungsstrom als Sili- 
ziumdioden. Mit den schon erwahnten Dioden OA 180 
wird bei einer Dynamik von 60dB eine obere Grenz- 
frequenz von 50 kHz erreicht, die sich bei einem breit- 
bandigeren Verstarker noch etwas hinausschieben lieBe. 
Bei 20 kHz erhalt man sicher 75 dB Dynamik, bei tie- 
feren Frequenzen noch etwas mehr. Die untere Grenz- 
frequenz von 20 Hz wird durch den zur Ankopplung 
des Steuerverstarkers erforderlichen Kondensator be- 
dingt (Schaltbild Bild 2). 

Soll der Verstarker dazu benutzt werden, exponen- 
tiell abklingende Schwingungen zu logarithmieren, so 
darf man bei periodisch wiederholtem Ausklingvorgang 
die Wiederholungsfrequenz nicht zu hoch wahlen. Durch 
die ersten Schwingungen des Nachhallvorgangs wird 
ndmlich unter Umstaénden der Koppelkondensator um- 
geladen; diese Ladung sollte bis zum Eintreffen des 
nachsten Schwinguhgszuges abgeflossen sein. 

Einige Details der Schaltung werden an Hand yon 
Bild 2 erlautert. Die beiden Rohren des Rechenverstar- 
kers sind durch eine kleine Glimmlampe gekoppelt. 
Durch eine Zenerdiode im Katodenkreis der Pentode 
wird erreicht, daB Eingang und Ausgang des Rechen- 
yerstarkers gleiches Potential aufweisen. Mit dem Reg- 
ler in diesem Katodenkreis wird der korrekte Arbeits- 
punkt eingestellt. Durch einen Widerstand im Katoden- 
kreis der letzten Rohre kann der Innenwiderstand des 
Verstaérkers auf den Wert gebracht werden, der zur 
Kompensation des Bahnwiderstandes der Dioden er- 
forderlich ist. Um die Vorspannung der Dioden zu ge- 
winnen, bedienen wir uns eines einfachen Schaltungs- 
kniffes, so daB wir ohne zusatzliche (und erdfreie!) 
Spannungsquelle auskommen. Der Katodenstrom der 
letzten Stufe (3 mA) flieSt durch eine Ge-Flachendiode, 
die in Durchlafrichtung gepolt ist, und erzeugt an ihr 
die gewiinschte Sperrspannung ftir beide Dioden (hier 
300 mV). Der differentielle Widerstand der Diode bei 
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Bild 3. Die Kennlinie des logarithmischen Ver- 
stirkers. Bei 0 dB ist Up=2 V effektiv. 
20kHz, ——-— 5kHz. 
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diesem Strom (etwa 15 Ohm) ist wesentlich geringer 
als der Quotient aus Spannung und Strom; daher 
andert sich die Vorspannung nur wenig mit dem Kato- 
denstrom. 

Die Temperaturabhangigkeit der Anordnung ist 
durch die Anderung von wy und J, der Dioden mit der 
Temperatur gegeben, wobei sich die Empfindlichkeit 
des Verstarkers, das heifft die Ausgangsspannungs- 
anderung fiir ein gegebenes Spannungsverhdltnis am 
Eingang (zum Beispiel in mV/dB), nur mit up=k T/e 
andert, wahrend die Variation yon Js die erzielbare 
Dynamik beeinfluBt. Zwar steigt der Sattigungsstrom 
I;, der durch thermische Paarbildung von Elektronen 
und Defektelektronen bedingt ist, stark mit der Tem- 
peratur an, doch ist die gewahlte Vorspannung der 
Dioden so hoch, daB die Dynamik bis zu Temperaturen 
yon 35 °C noch nicht merklich beeintrachtigt wird. Die 


Bild 4. Exponentiell abklingende Schwingung der Fre- 
quenz 15 kHz und der Abklingzeit 0,56 ms; im 
unteren Bild die logarithmierte Schwingung; 
Zeitmastab 0,5 ms/Skalenteil. Die Gesamtbreite 
entspricht einem Abfall um 77 dB. 
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relative Anderung der Empfindlichkeit ist im iibrigen 
gleich der relativen Anderung der absoluten Tempera: 
tur, das heift fiir 10 Grad Temperaturanderung be- 
kommen wir 3,3% Fehler in der Empfindlichkeit, das 
ist ein Wert, der im allgemeinen zu tolerieren ist. Eine 
Thermostatierung der Dioden eriibrigt sich. 

Fir stationdéren Betrieb zeigt Bild 3 die Ausgangs- 
spannung, als Spitzenwert gemessen, in Abhangigkeit 
yon der logarithmisch aufgetragenen Eingangsspan- 
nung. Erst fiir eine Dynamik > 75 dB_unterschei- 
den sich die Kurven fiir 5 und 20 kHz. Die Photogra- 
phien (Bild 4) schlieBlich zeigen Nachhallvorginge; 
im oberen Bild sieht man die ausklingende Schwingung 


Bild 5. Logarithmierte Nachhall-Schalldruckkurve eines 
mit Terzrauschen angestofenen Modellhallrau. 
mes; ZeitmaRstab 0,1 s/Sktl., Nachhallzeit 0,6 s. 
Terzband mit 10 kHz Schwerpunkt. 


eines dimpfungsbehafteten elektrischen Schwingkreises. 
im unteren die logarithmierte Spannung. Bild 5 zeig! 
die logarithmierte Nachhall-Schalldruckkurve eines mit 
Terzrauschen angestoBenen Modellhallraums. 


Der Deutschen Forschungsgemeinschaft haben wit 
fiir eine Sachbeihilfe zu danken. 
(Eingegangen am 13. Januar 1961.) 
K. H. Zimmermann und 


H. D. Gruscaxa 
Ill. Physikalisches Institut 
der Universitat Gottingen 
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A note on phase distortion and hearing 


The subject of phase distortion in sound reproduction 
is relatively little understood. Many times it has been 
claimed that pure phase distortions are detectable. 
Other authors, however, contend that the audibility of 
phase distortions can always be explained by the associ- 
ated amplitude response. The discussion has concerned 
not only the reproduction by a loudspeaker but the re- 
sponse of amplifiers as well: objectionable phase distor- 
tion in the audible range has been thought to originate 
in the amplifier response at ultra-high frequencies. 


Only a few experiments have been dedicated to the 
problem of detecting phase changes that are not accom 
panied by changes in the amplitude response. In suck 
cases the usual result was that pure phase distortior 
had only a slight or no effect at all upon the reproduc 
tion. The following remarks may be of value in this 
connection; the views to be developed have arisen from 
a study of the so-called “residue” effect [1]. 

In the category of continuous sounds tonal signal: 
are important, because they have a fundamental pitch 
that can clearly be distinguished from the timbre. (This 
distinction is not so clear in noisy signals.) For tonal. 


: 
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although not necessarily harmonic, signals the so-called 
linear phase rule has been formulated [2]. This rule 


_ states that the timbre of a sound does not change, when 


the phases of the components are shifted by amounts 
that are linearly dependent on frequency. Thus when 
the phase p(f) of the component with frequency f is 
changed by an amount 


Ag(f) =a+6f 


no change of the timbre will occur. The evidence for 
this rule comes from experiments on signals of which 


(a and 6 being constants) 


‘the components have equal frequency spacings. Such 


signals, be they harmonic or not, have a steady sound 
quality. 

The term 6 / in the formula is easily understood. In 
fact it represents a simple shift in time equal for all 
components. Fo a steady sound this is, of course, not 
noticeable. The important term is the constant a, ex- 
pressing a constant phase shift of all components. The 
phases of the components can all be changed by the 
same amount without bringing about a change in 


timbre. The experiments, on which this statement is 


based, are fundamentally carried out as follows. A har- 
monic sound complex is set up by a modulation pro- 
cess. A carrier of, say, 1600c/s is modulated by a 
signal that contains components of, e.g., 200, 400, and 


600 c/s. Note that the components of the modulating 


signal are to be harmonically related. The modulating 
signal thus has a fundamental frequency of 200 c/s. If 
now the carrier is also a harmonic of 200 c/s, the result- 
ing signal will be purely periodic. It then consists ex- 
clusively of harmonics of 200 c/s, namely 1000, 1200, 
1400, 1600, 1800, 2000, and 2200c/s. This signal 
sounds as a sharp tone of a pitch of 200 c/s. 

If the carrier frequency (prior to modulation) is now 
shifted a very small amount in frequency, all compo- 
nents of the resulting signal incur the same frequency 
shift. We can describe this by saying that, with respect 
to their former positions, the components have acquired 
a common phase shift that increases proportionally 
with time (Fig.1). Various values of this common 


/ 
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1200 =1400 §=1600 1800 2000 


Fig. 1. Signal of which the components undergo a com- 
mon phase shift. 


phase shift thus alternately present themselves. If such 
a common phase shift produced an audible change, we 
would hear beats (in this case slow variations of 
timbre). Nothing of the kind occurs, however, and we 
must conclude that such a uniform phase shift is not 
detectable. We have studied this phenomenon for many 
frequency ranges, and the rule seems to hold always. 
The exceptions could all be traced down to artefacts giv- 
ing rise to non-uniform phase shifts. 

Although no complete evidence is available, it is 
expected that the rule holds also for transient signals 
(filtered and modulated pulses). It remains to be stu- 
died for the transient onset of tones but this presents 
large difficulties both experimentally and methodically. 
The rule can be described in different terms if we refer 
to the waveform of the signal. Because of the modula- 
tion, the signal has acquired an envelope which has the 
waveform of the modulating signal. The carrier can be 
recognized in the “fine structure” of the signal, the 
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fluctuations of the waveform between the limits set by 
the envelope. In Fig. 2 this is shown for various values 
of the phase of the carrier signal with respect to the 
modulating signal. A shift of the phase of the “fine 
structure” with respect to the envelope is thus not 
noticeable. All signals of Fig. 2 thus sound the same. 


Fig. 2. Waveforms of signals with various values of the 
phase of the carrier signal with respect to the 


modulating signal. 


If we now return to problems of sound reproduction, 
we can assume that the rule holds for all signals. One 
can then draw the following conclusions regarding the 
effects of pure phase shifts on transmission. Refer- 
ring to Fig. 3, the phase-versus-frequency characteristic 
of a (non-minimum-phase) transducer +, we note that a 
phase shift as shown can approximately be regarded as 
a combination of a linear phase shift (meaning time 
delay) and a uniform phase shift (fine structure delay) . 
Both types pass unnoticed, hence a phase shift as shown 
in the figure cannot be detected. This holds as long as 
the approximation by a straight line is valid. Curvature 
of the phase characteristic will eventually lead to a 
detectable phase change. In the modulated signals 
described above, such a curvature can be simulated by 
shifting the phases of the modulating signal compounds 


1 This term refers to properties of linear networks. 
It can be shown that networks having no zeros of trans- 
mission in the right half of the complex-frequency 
plane, have the smallest phase shift of all networks 
with the same amplitude characteristic. In addition the 
minimum-phase characteristic can be computed from 
the amplitude-versus-frequency characteristic. For such 
networks the phase characteristic obviously cannot be 
changed without changing the amplitude characteristic. 
To study the audibility of phase changes only, one 
should use non-minimum-phase networks and _trans- 
ducers. 
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Fig. 3. Phase-versus-frequency characteristic of an elec- 
troacoustic transducer. 
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Fig. 4. Non-uniform phase shifts of the components of a 
signal. 


with respect to each other (Fig. 4)?. It is found that 
such changes are detectable only if the component fre- 
quencies lie within approximately 20 percent of each 


2 Similar phase changes can be brought about by 
changing the central component (carrier) with respect 
to the side bands. See Maruss, R. C. and Mirzer, R. L., 
Phase effects in monaural perception, J. Acoust. Soc. 
Amer. 19 [1947], 780. 
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other. The critical band mechanism seems to be re- 
sponsible for this. 

Concluding, then, it is.to be expected that a phase 
characteristic showing no appreciable curvature within 
the critical bandwidth of the ear, does not lead to 
audible distortion. This is clearly true for steady signals 
and has partially been proven for transient signals. It 
remains to be studied for the transient onset of a tone. 
At the moment we cannot safely state what the thresh- 
old of detectability is in terms of the second derivative 
d?y/df? per critical band. The number of components 
will almost certainly influence this threshold. 

With the evidence at hand it seems highly improbable 
that phase distortions, caused by irregularities of the 
response in the ultrasonic region of amplifiers, would 
ever be audible. If high-frequency resonances are de- 
tectable at all, this should be caused by distortion or 
instability due to the excitation of these resonances. The 
effect should disappear when the incoming signal is 
truly confined to the audible frequency band and 
distortion reduced to a minimum. 

(Received 27th October 1960.) 
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Sixth International Congress of Audiology 


The sixth International Congress of Audiology will 
be held in Leyden, The Netherlands, from September 
5—8, 1962. 


President: Prof. Dr. H. A. E. van Disnoxrcx; 
Secretary: Drs. A. Spoor. 


In the program are planned three round-table talks, 
with associated and free papers. 


The subjects of the round table talks are: 


“Frequency analysis of the normal and pathological ear” 
Moderator: Prof. Dr. G. von BExEsy. 


“Central deafness in children” 


Moderator: Prof. Dr. J. M. Taro. 


“Psychogenic deafness and simulation” 
Moderator: Prof. Dr. H. A. E. van Disnorck. 


Official languages of the Congress are: English, 
French, German and Spanish. Working language will 
be: English and French. 

The address of the secretariat is: 

Ear-Nose-Throat Department, 
Academisch Ziekenhuis, 
Leiden (The Netherlands) 


